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1 . 0 SUMMARY 


The objective of this program was to develop a thermal load transfer 
module that would be capable of transferring thermal data from a three dimen- 
sional analytical model of a combustor liner to a corresponding three dimen- 
sional stress analysis model. The initial goal of the transfer module was the 
ability to handle different mesh densities for the heat transfer analysis and 
the stress analysis and the capability to utilize both finite difference and 
finite element heat transfer codes. The extended goals of the transfer module 
were to develop a general purpose tool that would effectively allow for the 
coupling of heat transfer analysis and stress analysis while maintaining the 
flexibility of independent analysis codes. 

The software package developed under this contract is called Transfer 
Analysis Code to Interface Thermal/Structural Problems or TRANCITS. The 
TRANCITS code satisfies all of the objectives for transferring thermal data 
between heat transfer and structural models of combustor liners and, in addi- 
tion, can be used as a generic thermal translator between heat transfer and 
stress models of any components regardless of the geometry. TRANCITS can 
accurately and efficiently convert the temperature distributions predicted by 
the heat transfer programs to those required by the stress codes. It can be 
used for both linear and non-linear structural codes and can produce nodal 
temperatures, elemental centroid temperatures or elemental gauss point temper- 
atures. The thermal output of both the MARC and SINDA heat transfer codes can 
be interfaced directly with TRANCITS and it will automatically produce stress 
model temperatures formatted for NASTRAN and MARC. In addition to these codes, 
any thermal program and structural program can be interfaced by utilizing the 
neutral input and output forms supported by TRANCITS. 

In summary, the TRANCITS code can be used to interface temperature data 
between thermal and structural analytical models. The use of this transfer 
module allows the heat transfer analyst to select the thermal mesh density and 
thermal analysis code that is best suited to solve the thermal problem and 
gives the same freedoms to the stress analyst without the efficiency penalties 
associated with common meshes and the accuracy penalties associated with the 
manual transfer of thermal data. 
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2.0 INTRODUCTION 


The overall objective of this program was to develop a thermal data trans- 
fer computer program module for the Burner Liner Thermal-Structural Load Model- 
ing Program. This was accomplished by (1) reviewing existing methodologies 
for thermal data transfer and selecting three heat transfer codes for appli- 
cation in this program, (2) evaluating the selected codes to develop criteria 
for developing a computer program module to transfer thermal data from the 
heat transfer codes to selected stress analysis codes, (3) developing the 
automated thermal load transfer module, and (4) verifying and documenting the 
module . 

In aircraft turbine engine hot section components, combustor liners, and 
turbine blades and vanes, cyclic thermal stresses are the most important damage 
mechanism. Consequently, accurate and reliable prediction of thermal loads is 
essential to improving durability. To achieve this goal, a considerable effort 
over the past 20 years has been devoted to the acquisition of engine temper- 
ature test data, as well as the development of accurate, reliable, and effi- 
cient computer codes for the prediction of steady-state and transient temper- 
atures and for the calculation of elastic and inelastic cyclic stresses and 
strains in hot section components. There is a need for continued development 
of these codes, because the availability of more accurate analysis techniques 
for complex configurations has enabled engine designers to use more sophisti- 
cated designs to achieve higher cycle efficiency and reduce weight. 

It has become apparent in recent years that there is a serious problem of 
interfacing the output temperatures and temperature gradients from the heat 
transfer codes with the input to the stress analysis codes. In part, this is 
a penalty for success. When computers were slower and computer memory smaller, 
the size of problems that could be analyzed, in terms of heat transfer nodes 
and stress analysis finite elements, was limited. Manual transfer of the 
output temperatures from the heat transfer codes to the stress analysis input, 
with manual interpolation to accommodate the mismatch in the nodal meshes of 
the two programs, was not unduly burdensome. With the growth in computer capa- 
city and speed and the development of input preprocessors and output postpro- 
cessors, the analysis of components using hundreds and even thousands of nodes 
in the heat transfer and stress models has become economical and routine. This 
has exacerbated the interfacing problem to where the engineering effort 
required is comparable to that required for the remainder of the analysis. 
Furthermore, a considerable amount of approximation has been introduced in an 
effort to accelerate the process. This tends to introduce errors into the 
temperature data which negates the improved accuracy in the temperature distri- 
bution achieved through use of a finer mesh. 

One solution to these interfacing problems is to force the heat transfer 
model and the stress model to be performed using the same analysis code and 
the exact same mesh densities. Although this technique may be acceptable for 
some applications, it has several disadvantages. In many models the areas of 
high stress gradient may be influenced by thermal loads but are predominantly 
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controlled by geometrical discontinuities such as fillets and holes. A struc- 
tural mesh that was fine enough to represent the gradients due to these dis- 
continuities would be in general much denser than necessary for the thermal 
problem. This "overkill" might be tolerated for 2D analysis but many times 
the additional cost associated with these extra nodes in 3D analysis cannot be 
afforded. It is sometimes argued that because of the reduced degrees of free- 
dom for a thermal problem as compared to a structural problem (typically 1 dof 
for thermal and 2 to 6 dofs for structural) that this additional cost is mini- 
mized. This is true for steady-state analyses, but for transient problems, 
the number of thermal runs far exceeds the number of structural runs usually 
performed. 

Another example of the disadvantage of common meshes is the current trend 
in structural analysis toward adaptive mesh refinement. This technique will 
automatically refine the structural mesh in areas of high strain gradients. 
Unless the high strain gradient were caused by high thermal gradients there is 
certainly no need to refine the thermal mesh and redo the thermal analysis 
just because the structural mesh was made denser. All of these factors point 
toward the need for different thermal and structural mesh densities. 

A further drawback of common heat transfer and stress codes and meshes is 
the widespread use of the finite difference technique to solve the thermal pro- 
blem as opposed to the use of finite element codes for the structural problem. 
Again this suggests the need for an external process to transfer the thermal 
data to an independent structural analysis. 

Therefore, the primary objective to this project is to develop, implement 
and demonstrate an interfacing module that will address the problems associated 
with the transfer of temperature data from thermal analyses to structural 
analyses . 

This interfacing module must be able to handle different thermal and 
structural meshes and must support both finite difference and finite element 
heat transfer codes. 

Conceptually, the overall program structure is shown in Figure 1. 

Direct inputs are temperature data, and the geometries associated with the 
thermal data and with the stress model. The particular heat transfer and 
stress analysis programs used affect the process inside the transfer module. 
Final output is the temperature data at the required locations in the stress 
model . 
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3 . 0 DESCRIPTION OF TASKS 


This project was addressed by completing four tasks. These tasks include 
a review of the current status of thermal load modeling, an assessment of the 
thermal load modeling procedures, development of a thermal load transfer module 
and verification and documentation of the transfer module. 

The four tasks are described below. 


3.1 TASK I - CURRENT STATUS OF THERMAL LOAD MODELING 

General Electric conducted a review of the current methods used to trans- 
fer thermal data from analytical sources to structural analysis programs in 
order to compute thermal stress and strain. The review included both in-house 
capabilities and papers in the open literature. The feasibility of the avail- 
able methods were assessed and incorporated into the thermal load transfer 
module when applicable. In addition to reviewing the existing transfer tech- 
niques, General Electric also conducted a review of existing numerical temper- 
ature prediction codes including both in-house programs and those available in 
the open literature. This review was conducted with particular attention paid 
to the existing thermal load transfer methodologies and also included an exami- 
nation of the temperature output format for each code reviewed. Based upon the 
results of this review, three computer codes were identified from among those 
available to the public to be used in the verification of the thermal load 
transfer module which was developed in Task III of this contract. 

A review was also made of structural analysis computer programs as the 
first step in the process of selecting those codes to be used in this program. 
Both in-house computer programs and those programs defined in the open litera- 
ture were reviewed. The objective of this review was to produce an evaluation 
of each code’s input-output formulation, particularly with respect to tempera- 
ture information. For the greatest public utility, primary emphasis was placed 
on those codes which are both in the open literature and actively supported by 
the commercial computer houses. 


3.2 TASK II - ASSESSMENT OF THERMAL LOAD MODELING 


In Task II, the Aircraft Engine Business Group outlined the thermal data 
transfer procedures for each of the methodologies approved in Task I and 
identified, discussed, and assessed their limitations and deficiencies in 
order to clearly define and evaluate the problem areas associated with thermal 
load modeling. Particular emphasis was placed on the questions of compati- 
bility, ease of implementation, and thermal data format in the discussion of 
the current deficiencies and limitations involved in the transfer of thermal 
data from heat transfer analysis to the various structural analysis codes. 
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Once the limitations and deficiencies were identified in this task, a 
consistent set of evaluation criteria was developed. These criteria included 
factors that address the accuracy and the efficiency of the transfer procedure. 


3.3 TASK III - DEVELOPMENT OF THERMAL LOAD TRANSFER MODULE 


Based on the results of Tasks I and II, General Electric developed a 
thermal data transfer module that efficiently and accurately interfaced 
thermal data into structural analysis codes. The thermal data may come from 
both finite difference and finite element heat transfer analysis programs. 

The transfer module is capable of interfacing 3-D steady-state and transient 
heat transfer data into linear and nonlinear structural analysis programs. 

State-of-the-art features in this module allow different mesh distribu- 
tion for the heat transfer model, and the structural analysis model and slight 
differences in part geometry due to differences in component tolerances used 
in the models. Three-dimensional translation and rotation routines were 
incorporated into the module to account for possible mismatches between the 
origin of the structural model and the thermal model coordinate systems. 

As these techniques were developed, they were encoded into a modular 
software package that efficiently manages the thermal data. The temperature 
data from the selected heat transfer programs is read directly by the inter- 
facing module and automatically converted into a usable form for the selected 
structural analysis codes. The output of the interfacing module is completely 
compatible with the requirements for all of the structural codes selected. In 
addition, the flexible format of the output allows it to be easily implemented 
into other structural codes. 

The techniques developed for the transfer module were designed to be 
flexible and have as much built-in generic capability as possible. The module 
provides options for compatible input with different existing heat transfer 
and structural computer codes. In conclusion, the transfer module provides an 
accurate streamline method of determining burner liner temperatures and gradi- 
ents at specified locations compatible with the structural analysis input 
requirements . 


3.4 TASK IV - VERIFICATION AND DOCUMENTATION OF AUTOMATED THERMAL TRANSFER 
LOAD MODULE 


This Task is concerned with validating the data transfer module which was 
developed in Task III. Several test cases were run utilizing various combina- 
tions of heat transfer and structural analysis computer codes with temperature 
data flowing between the programs by means of the data transfer module. One 
of the test cases involved a 3D model of a combustor liner. Both finite 
difference and finite element heat transfer codes were used to produce the 
thermal results and these thermal loads were then interfaced to structural 
analysis codes using the transfer module. Cases were run using similar and 
dissimilar meshes for the heat transfer and stress analysis. Stress points 
that lie slightly outside the heat transfer model were interfaced using the 
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module. The accuracy of the transfer was evaluated by comparing the tempera- 
tures produced by the transfer module to those predicted by the actual heat 
transfer codes used. 

These demonstration test cases were documented as verification of the 
transfer code. 

The coding of the transfer module was made as hardware independent as 
possible. All computer system dependent routines were isolated into easy 
to modify subroutines. The final version of the transfer module is executable 
on NASA Lewis's IBM system. 

As a separate document a User's Manual for the data transfer module was 
provided. This manual includes a precise set of operating instructions for 
each of the software packages developed. Included is a description of all 
variables needed to execute the codes. 

As part of this documentation, a complete set of Fortran listings for 
each subroutine developed for the data transfer module along with a computer 
tape of this coding was delivered to NASA Lewis Research Center. 
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4.0 TASK I 


Current status of thermal load modelling and a review of heat transfer 
and structural codes. 


4.1 THERMAL ANALYSIS PROGRAMS 

A review was conducted of the in-house and open literature concerning 
numerical temperature prediction codes. The primary purpose of this task was 
to select three codes to be used to verify the thermal load transfer module. 

The other criterion for the codes is that they must be available in the NASA 
Lewis computer system. 

A comprehensive summary of the capabilities of thirty-eight different 
codes is contained in the excellent survey article by A.R. Noor entitled "Sur- 
vey of Computer Programs for Heat Transfer Analysis". (Reference 1) We studied 
this summary closely, paying particular attention to any already existing 
transfer modules between the thermal and structural analysis modules, as well 
as to the output form for the temperature results. 

As to the first question, none of the programs surveyed contained an 
automatic thermal interface. This statement is also applicable to all our 
in-house heat transfer programs as well. Most major finite element programs 
do, however, have the capability to pass thermal data directly to an identical 
mesh which can then be used for structural analysis. Such a capability may be 
extremely useful for certain applications, but for many applications it is 
very restrictive. It imposes mutual restrictions on the heat transfer analyst 
and on the structural analyst which neither of them really need to observe. 

It is also true that this one-to-one correspondence requires that the heat 
transfer program use a finite element technique. This fact is clearly not a 
generally desirable arrangement. 

The second area we focused on was how we could access the temperatures 
produced by the heat transfer code. Most of the major codes surveyed also 
carry a file of output temperatures, which are suitable for post-processing. 

The form of this file is, of course, program dependent, but we believe that it 
can become the basic input to the thermal loads interface module. A thermal 
output file is not absolutely necessary since all the thermal codes print out 
the temperature at the heat transfer elements in some readable form. This 
readable form, known as "hard-copy", can be scanned and the required thermal 
data can be obtained. 

Since all the programs considered possess both the special output file as 
well as the standard "hard copy" output, either form of the output could be 
utilized. One of the unknown questions is: Does the special output file 

contain all the information necessary for the thermal load interface? Fre- 
quently finite difference based programs, for example, do not carry a complete 
geometry definition as a part of the output. The complete geometry is essen- 
tial to the interface module. It can, in most cases, be obtained from the 
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input to the heat transfer program. Questions such as this one cannot be 
answered with confidence until the individual programs are actually tried. 

In general, all programs were found to be roughly equivalent as far as 
producing temperature data usable by a structural program. Many different 
features are available in these codes, for example, different treatments of 
radiation, better plotting capabilities, etc., but the main focus of this 
study is on transferring temperature data. On this basis, there are nine 
major programs which could be used: 

ANSYS 

HEATING 

MARC 

NASTRAN 

SINDA 

THTD 

ADINAT 

SAHARA 

PAFEC 

Of these codes, only four are currently available on NASA Lewis computers and 
these are HEATING, MARC, NASTRAN, and SINDA. Further, the transfer module is 
demonstrated with the interface of both finite difference and finite element 
codes, namely the SINDA program for the finite difference program and MARC as 
the finite element program. Also demonstrated is the interface to NASTRAN, 
because of its wide popularity. For completeness, the portions of Noor's 
survey article which deal with these programs and others are listed in Appen- 
dix A. 

Other codes considered, but not available on NASA Lewis f s computer, were 
ANSYS, ADINAT, SAHARA and PAFEC. ANSYS and ADINAT are both widely used finite 
element codes. Both of these codes have a wide range of elements that can be 
used for heat transfer analysis, and ANSYS is frequently used for gas turbine 
design problems. It would be desirable to have the transfer module deal with 
these heat transfer codes directly. The SAHARA program developed at Sandia 
laboratories is another finite difference program that could be used to vali- 
date the transfer module. The survey indicates that SAHARA, like SINDA, has an 
output file that can be used for postprocessing. This output file would be 
very useful for the transfer module. Finally the PAFEC program was considered 
because it uses boundary integral techniques to perform the heat transfer 
analysis. Here again, the ability to transfer thermal data from such a pro- 
gram is beyond the requirements of this contract but nevertheless will be very 
desirable. It is felt, however, that because of the modular construction of 
the transfer module, future implementation of these codes or any other common 
code could be easily accomplished. 


4.2 STRUCTURAL ANALYSIS CODES 

A survey of stress analysis programs was also conducted. By virtue of 
the huge number of programs available (it is certainly in the hundreds), our 
survey cannot be considered comprehensive, but our conclusions are consistent 
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with the scope of this study. In particular, no structures program which was 
considered carries its own thermal interface transfer module. Temperature 
data are transferred by one of two methods. 

(a) Temperatures are obtained from an exactly equivalent heat transfer 
analysis, utilizing the identical mesh, 

or 

(b) temperatures are interpolated and/or extrapolated from a different 
mesh to discrete points: nodes, gauss points, element centroids, 
etc. 

Method (a) requires no interface module whatsoever, but it is not consid- 
ered to be desirable for a general case. 

Method (b) is ideally suited for a thermal interface module in that it is 
quite general. The module which was developed in Task III produces the temper- 
ature at any point in space which is consistent with a specific stress analy- 
sis. Whether this point is an actual node or gauss point makes little differ- 
ence. The actual details of which points are required for temperature input 
to a structural model are not critical. The module will accept any selection. 

Our survey of commercially available programs relied heavily on a recent 
review with R. Zirin of the General Electric Gas Turbine Division covering 
ANSYS, EASE 2 , MARC, NASTRAN, STARDYNE, and SUPERB. 

All of these finite element structural codes have a large library of 
different elements. The location within these elements to which the thermal 
data must be input varies. Some require nodal temperatures, some need elemen- 
tal centroid temperatures, others need gauss point data. They all, however, 
need temperatures at some discrete point in the structural model. Once again, 
the transfer module will produce temperatures at any location. 

For the purposes of this contract, the transfer module is demonstrated on 
MARC and NASTRAN, the codes available at NASA Lewis. Again, for completeness, 
the information on MARC and NASTRAN is given in Appendix B. 

The thermal load interface module developed under this study is flexible 
enough to drive any structural program which needs temperatures for input. 
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5.0 TASK II 


Assessment of Thermal Load Modeling 

The objective of Task II was to evaluate the thermal transferring capa- 
bilities currently available in existing heat transfer codes. In addition, an 
attempt was made to identify any freestanding thermal data transfer modules 
that might be in use by industry. Once these capabilities were evaluated, a 
set of criteria was compiled that included these factors and several additional 
criteria that were important in the development of a successful transfer module. 

As stated in Section 4.1 no general capability for transferring 3D thermal 
data to structural codes was identified in the thermal codes surveyed. The 
most common mechanism for accomplishing this interfacing was through the use 
of similar mesh densities. One exception to this was the THTD program that did 
produce a thermal output file that could input temperatures directly into 
selected stress codes. The actual procedure of assigning the stress node tem- 
peratures, however, involved a manual correlation between the stress point and 
the heat transfer element chosen to be "closest" to the stress node. Not only 
was this time consuming, but it also did not provide for any interpolation 
between heat transfer elements. The proposed code would consider this short- 
coming in its development. 

The investigation of external thermal load transfer modules identified 
three codes being used by various companies. These codes were 2D TITAN devel- 
oped by General Electric (AEBG) , MERLIN developed by Sandia Lab and LOMAP 
developed by Lockheed. The evaluation of each of these codes follows: 


5.1 2D TITAN EVALUATION 


The 2D TITAN transfer code was developed by GE to interface the thermal 
results of a heat transfer code called THTD to selected structural analysis 
codes. The THTD code is a finite difference formulation. The transfer code 
will interface temperatures into any 2D mesh including plane stress, plane 
strain, axisymmetric and shells of revolution. It has some heat transfer win- 
dowing capability both spatial and temporal. TITAN can interface stress points 
that lie slightly outside of the heat transfer model. It also has built in 2D 
coordinate transformations to align the thermal and structural models. The 
program has some internal accuracy checks used to flag areas of the thermal 
model where the interfacing could yield inaccurate results. Another related 
feature is the use of interpolating functions that are consistent with the 
assumption associated with the THTD heat transfer element. The actual inter- 
polation is accomplished using the 2D planar isoparametric shape functions. 

The following limitations were identified: 

a. No 3D capability. 

b. Only interfaces with one heat transfer code. 
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c. No time interpolations between thermal solutions, 

d. No provision for scaling the thermal results based on variations in 
thermal boundary conditions. 


5.2 MERLIN EVALUATION 

The Sandia code (MERLIN) is very similar to the in-house GE program called 
TITAN. Investigation into MERLIN revealed the following limitations: 

a. Restricted to 2D problems. 

b. Interfaces only with finite element heat transfer codes. 

c. Does not interface with finite difference heat transfer codes. 

d. Only interfaces into a specific set of finite element stress codes. 

e. Meshes (both heat transfer and stress) must be generated using a 
specific mesh generator. 

f. Does not allow for coordinate transformation of meshes. 

g. Does not account for stress points " slightly” outside heat transfer 
mesh. 

h. Limits the number of elements in the heat transfer model. 

i. Limits the number of nodes in stress model. 

j . Uses the same temperature mapping function regardless of the type of 
heat transfer element. 

The interfacing technique used in MERLIN is a three-step process. Step one 
is the mesh search. This search locates which heat transfer element contains 
each stress point. The second step is called the element search. This process 
defines the local coordinates of the stress point with respect to the coord- 
inates of the heat transfer element. The last step computes the temperature 
of each stress point based on the weighting coefficients derived from the 
local coordinates. This entire process is almost identical to the methods used 
in TITAN. 


5.3 LOMAP EVALUATION 


The Lockheed program LOMAP has very limited capabilities as a general 
transfer module. It will transfer temperatures from a 3D heat transfer mesh to 
a 3D stress mesh, however, the method used is very simplistic. The weighting 
coefficients used to interpolate temperatures are Proportional to the temper- 
atures of the four closest heat transfer nodes. For a general transfer program 
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; 


the four closest nodes may not be the correct nodes to use. In fact, the four 
geometrically closest heat transfer nodes may not even be in the same thermal 
region as the stress point. This technique places severe restriction on the 
types of meshes that can be interfaced properly. 


5.4 EVALUATION OF TRANSFERRING CRITERIA 

Based on the study of existing thermal transfer modules, GE internal 
needs, and GE previous experience with TITAN, three levels of criteria for the 
program development were proposed. 

Level I contains the general criteria which must be satisfied for a usable 
product. 

Level II contains specific criteria which must be satisfied to meet the 
requirements associated with gas turbine design problems. This list stems 
mainly from GE internal experience. 

Level III contains criteria which are desirable but not necessary. In most 
cases, items in this class can be achieved through a multi-step Process. Total 
automation might be desirable, but we do not believe this effort is warranted 
at this time. 


Level I: General Criteria For A Thermal Transfer Module 


IA. Independent Heat Transfer and Stress Geometry Meshes 

This criteria lies at the heart of our effort. Useful thermal transfer 
modules must address this feature. Automatically included in this is the 
ability to transfer from finite difference heat transfer mesh to a finite 
element stress analysis mesh. 

IB. Accurate Transfer of Data 

Simplistic approaches such as averaging the closest nodes do not always 
yield accurate results, and the utility of the transfer program is then ques- 
tionable. This criteria was met by using all available temperature information 
to do the interpolation and by using different mappings to correspond with 
different heat transfer elements. 


IC. User Friendly 

Programs which do not meet this criteria tend to be used incorrectly or 
as a last resort. The plan was to construct the thermal transfer module to 
encourage the analyst to use it. Any errors encountered by the module are 
flagged and reported with a diagnostic message. 
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ID. Computationally Efficient 


The program was coded to achieve an efficient flow of data. GE’s past 
experience with TITAN has led to several improvements over the original 
efforts, and these have resulted in similar gains in 3D transfer problems. 
This criteria covers both searches to find the proper heat transfer element 
for a stress node as well as the single element inverse mapping functions. 


IE. Flexible 

The thermal transfer module was constructed such that future modifica- 
tions, or even different applications, for instance, pressure or boundary con- 
dition transfer, could be accomplished without a full rewrite. This criteria 
stems from past experience in having to improve or draw upon old techniques 
which could almost, but not quite, perform the required task. The transfer 
module not only transfers temperatures in a state-of-the-art manner, but it 
will also provide a vehicle for other 3D interpolation based problems. 


Level II - Specific Criteria Required For Gas Turbine Design Problem 

IIA. Coordinate Transforms 

Coordinate transformation that allow the heat transfer model to be aligned 
with the structural model. 


IIB. "Out-of-Box” Provision 

Provision to account for stress nodes that lie just outside the heat 
transfer model due to slight differences in the dimensions used in the heat 
transfer and stress analysis models, as a result of using different tolerances 
on the actual component dimensions. 


IIC. Windowing 

Capability to "window in" on a smaller portion of the heat transfer model. 


IID. Selected Time Steps 

Ability to select temperature distributions at specific time steps form a 
large transient thermal analysis. 

Level III - Desirable But Not Essential Features 

IIIA. Automatic Handling of Temperature Discontinuities 

In the module these will be treated in a two-step manner. Total auto- 
mation is possible, but has not been implemented. 
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IIIB. Scaling of Temperatures Based on Variation In Engine Power 
Level Settings 

Scaling of the interpolated temperatures based on changes in the heat 
transfer boundary condition, such as engine power settings, various inlet con- 
ditions or different cooling flows requires specific knowledge of the heat 
transfer analysis and will not be done inside the transfer module, but could, 
if desired, be applied by another program to the original results from the 
transfer module. 


IIIC. "Altered” Stress Geometry 

Many times the stress analyst wants to alter the part geometry to reduce 
his stresses, but the deviations will not, in the judgement of the heat trans- 
fer analyst, affect the temperatures. In the past, ”ad hoc" procedures to 
transfer temperatures to the new stress geometry have been used. This approach 
is not optimum and has not been included in the transfer module. 

Table I compares each of the external transfer codes evaluated against 
some of the criteria chosen. This table illustrates areas where significant 
effort was placed to develop an accurate and efficient thermal transfer module. 
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Table 1 


Criteria 

MERLIN 

Codes 

LOMAP 

GE 

TITAN 

TRANCITS 

Proposed 

Transfer 

Module 

Independent Meshes 

Yes 

Yes 

Yes 

Yes 

2-D Capability 

Yes 

Yes 

Yes 

Yes 

3-D Capability 
Thermal Mapping 
a Function of Heat 
Transfer Elements 

No 

Yes 

No 

Yes 

Coordinate 

Transformations 

No 

No 

Yes 

Yes 

‘Out-of-Box 1 

No 

No 

Yes 

Yes 

Provision 

Isoparametric 

Yes 

No 

Yes 

Yes 

Mapping 

Geometrical/Temporal 

No 

No 

Yes 

Yes 

Windowing Capability 
For Time Interpolation 

No 

Yes 

No 

No 




6.0 DEVELOPMENT AND DOCUMENTATION OF THE TRANSFER MODULE 


The objectives in this task were to develop an efficient and accurate 
software package that implemented all of the features discussed in Task I and 
Task II. The development was divided into several subtasks, each of which 
addressed important program features that were critical to a successful code. 

The subtasks were: 

1. General Flow of Code 

2. User Input 

3. General Heat Transfer Input 

4. Specific Heat Transfer Code Input 

5. General Stress Code Input 

6. General Output 

7. Specific Stress Code Output 

8. Finite Difference Heat Transfer Code Consideration 

9. Finite Element Heat Transfer Code Consideration 

10. General 3D Search Routines 

11. Weighting Coefficient Routines 

12. Exterior Surfacing Routines 


6.1 GENERAL FLOW OF CODE 

Several factors were considered in developing the general flow of the 
code. These factors included: 

a. Program modularity 

b. Program efficiency 

c. Input/output requirements 

d. Internal data flow 

e . Flexibility/portability 

f. User friendliness 


The basic program flow is illustrated by Figure 6-1. 

Each of the boxes in Figure 6-1 represent program modules that have been 
coded to minimize the system dependent features of the transfer program. In 
addition to this, these modules communicate with each other through an internal 
file structure that is hardware independent, easy to follow and easy to modify. 
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Figure 6-1. Overall Modular Flow of Code 








These features not only make the code very portable but also allow the same 
basic structure to be used in the future to transfer other types of loading. 
Additional heat transfer elements can easily be added to this scheme. 


6.2 INPUT AND OUTPUT 


The input and output information can be grouped into three basic cate- 
gories. They are geometrical, thermal and user control data. Figure 6-2 
illustrates how these data interact with the transfer module. 

The user control data are input to the transfer module with a namelist 
input file. This file contains all of the variables that can be input by the 
user to control the interfacing procedure and is described in detail in the 
TRANCITS Users manual (Reference 2). Options such as geometry and temporal 
windowing, coordinate transformations, time step tolerancing and information 
about the type of input and type of output desired are entered through this 
file and help provide for the user-friendliness of the code. 

The other two categories of input are associated with the geometrical data 
for both the heat transfer and stress models and the temperature data that is 
the output from a thermal code. The transfer module has been structured to 
accept output of the SINDA finite difference code and the thermal output of the 
MARC finite element program. In order to satisfy these requirements and at the 
same time produce a module that can be used with a variety of heat transfer 
codes, a neutral heat transfer input form was developed. This input form will 
support both finite difference and finite element thermal results and is the 
primary input to TRANCITS. Compatibility with SINDA and MARC is accomplished 
by internally converting the output of these codes to this neutral form. 

The heat transfer input file consists of six partitions of data. These 
data are described in detail in the users manual. It stores all the sizing, 
identification, qeometrical and thermal data associated with the heat transfer 
analysis. It is important to note that the geometrical data is an essential 
part of this input. The geometry (nodal coordinates and element connectivity) 
is the only link between the heat transfer and stress models. Without geomet- 
rical data an accurate interface cannot be accomplished. The thermal data in 
this file varies depending on the type of heat transfer analysis. Nodal tem- 
peratures are stored if the thermal code is a finite element type and elemental 
temperatures are stored if the thermal code is a finite difference formulation. 

The automatic interface with SINDA and MARC makes use of output files 
created by these programs. The MARC program produces an output file known as 
the MARC cost file. This output contains all of the data required by the 
transfer module to perform the interface. The SINDA code also has the capa- 
bility to produce an output file. However, since SINDA is a finite difference 
program, there is no geometry stored on the output file. Indeed, the SINDA 
code does not require geometry as input. Only volumes, areas and distances are 
needed to predict the thermal distributions. Therefore, the raw SINDA output 
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Figure 6-2. TRANCXTS Schematic. 














is not adequate for the transfer procedure. This file must be augmented by the 
SINDA geometry file. At first glance this seems to be a severe limitation of 
using TRANCITS in conjunction with any finite difference code. This is not the 
case because in practice large finite difference models are not created man- 
ually. Pre-processors automatically generate the volume and area data needed 
for the finite difference analysis. The input to these preprocessors is grid 
point coordinates and element connectivity. This input is exactly what is 
needed for the geometry portion of the TRANCITS input. 

The input associated with the stress model is much simpler than that of 
the heat transfer model. To interface nodal temperature for the stress model 
all that is required is a file containing the name and rectangular cartesian 
coordinates of the stress points to be interfaced. This information is con- 
tained in the stress model coordinate file. If elemental centroid or elemen- 
tal gauss point temperatures are requested, such as required by MARC, an 
additional file containing the stress element connectivity is also required. 

The primary output of the transfer module is a neutral temperature output 
file. This file consists of the name of the stress node and its corresponding 
temperature for every transient time point requested. This file can easily be 
input into simple formatting routines and the stress point temperatures can be 
configured into any form required. In addition to this neutral output file, 
the program has the capability to automatically format the temperatures 
compatible with the NASTRAN or the MARC analysis codes. 


6.3 FINITE DIFFERENCE CONSIDERATION 

There are two primary differences between using a finite difference code 
versus a finite element code in conjunction with the transfer module. The 
first is the potential lack of geometry associated with the finite difference 
codes. This has been discussed in the previous section. The second is the 
elemental locations at which the thermal results are predicted. The basic out- 
put for a finite element heat transfer code is the temperature at the nodal 
points. This is ideally suited for the interfacing procedure, since these tem- 
peratures along with the assumptions associated with the thermal field inside 
the element can be used to accurately interpolate data within the element 
domain. On the other hand, finite difference codes produce temperatures at the 
element centroids and perhaps at the center of the element faces. Temperatures 
at these locations are much more difficult to use effectively and pose a pro- 
blem. The approach taken in the transfer module is to convert these elemental 
temperatures into accurate temperatures at the elemental vertices. Once this 
has been accomplished, the interfacing technique proceeds in a manner similar 
to the one used for finite element heat transfer codes. This conversion of 
the elemental temperatures to nodal temperatures is performed in two steps. 

The first step is to map on an element by element basis the centroidal and 
face center temperatures to the element vertices. This mapping must be done in 
a fashion that is consistent with the assumptions governing the thermal distri- 
bution within the element. For a 3D linear finite difference volume, the 
assumption is made that the temperatures can vary linearly in all 3 directions 
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between the face centers and the centroid of the element. Given the tempera- 
tures at the 6 face centers and the centroid, the temperature, at the eight 
vertices of the element can be computed making use of the 3D tetrahedral iso- 
parametric shape functions. The equations that must be solved are of the 
form: 

T = [1, x, y, z] [c]' 1 T 

where 

T = temperature at vertice 


y = coordinates of vertice 
z 

T = known temperatures at the centroid and the three adjacent face 
centers . 

where 

1 x i y i z i -1 

[C]" 1 = 1 X T Y. r Z T 

J J J 

1 \ y k z k 

1 x T y t z t 

JL JU JL 

where 

X IJKL 

^IJKL = coordinate of the tetrahedral 
Z IJKL 

- 1 

Note that [C] is only a function of geometry and therefore must only be com- 
puted once regardless of the number of transient time steps. This inverse is 
computed in closed form in the transfer module and is therefore efficient and 
accurate. Figure 6-3 illustrates the finite difference volume and the corre- 
sponding tetrahedral. 

Note that during this initial stage the actual temperatures of the ver- 
tices are not directly computed, only the coefficients that relate the temper- 
atures at the vertices to those at the known locations are computed. 

The second step in this process is to compute a unique temperature at each 
vertice. This is done by weighting all of the vertice coefficients by a scale 
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factor that is inversely proportional to the distance between the centroid of 
each element adjacent to the vertice and the vertice itself. This technique 
makes use. of all of the surrounding temperature information to eventually 
produce a unique temperature at the vertice of each heat transfer element. 

Note that if the finite difference results are only at the centroid, the trans- 
fer module assumes the entire element was at that temperature. The weighting 
process however will generate different vertice temperatures since the temper- 
ature of the vertice depend on the geometry and the number of elements assoc- 
iated with that vertice. 


6.4 3D SEARCH AND WEIGHTING COEFFICIENT ROUTINES 


Two of the criteria that must be met by any transfer module are that the 
interfacing procedure be efficient and accurate. In order to satisfy the 
efficiency requirements, techniques were developed that rapidly determine 
which heat transfer element contains the stress points of interest. The 
approach taken was to invoke a "multi-step filtering" process. This process 
was developed and implemented such that most of the element candidates were 
eliminated using a very simple algorithm. The elements that passed through 
the initial "coarse" filter were then subjected to a more sophisticated test 
("fine filter") that determined if the stress point was contained within the 
bounds of the element. This search technique, like the vertice coefficient 
routine is only a function of geometry and therefore must only be performed 
once regardless of the number of transient solutions interfaced. The simple 
coarse filter requires that the minimum and maximum values of the x, y, z 
coordinates be stored for each heat transfer element. This procedure is out- 
lined below. 

The stress point is outside of the heat transfer element if any of the 
following conditions are met: 


X 

s 

X 


s 

Y 

s 

Y 

s 

Z 

s 

Z 


s 


< X . 

mm 

> X 

max 

< Y . 

mm 

> X 

max 

< Z . 

mm 

> Z 

max 


where X , Y , Z are the coordinates of the stress point. For most models 
this technique drastically reduces the number of heat transfer elements that 
could possibly contain the stress point. Note that this procedure is valid 
for any order 3D heat transfer element as long as the maximum and minimum 
values are based on all of the nodes associated with the element. 
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The technique employed as the "fine filter" actually performs two tasks. 
First it determines if the stress point lies outside of the element, in which 
case the next element is considered. Second, if it is determined that the 
stress point is contained within the element, it automatically returns coef- 
ficients that relate the known temperatures of the vertices of the element to 
the desired temperature at the stress point location. This "fine filter" or 
weighting coefficient routine is based on the inversion of the isoparametric 
shape functions. The choice of these interpolating functions was based on many 
factors. These functions are widely used for many finite element heat transfer 
element formulations. They can be formed to represent a variety of temperature 
fields and therefore can be used to describe the variation within a finite 
difference volume. They are mathematically easy to manipulate and understand. 
The inversion process involves solving for the so-called local isoparametric 
coordinates (r, s, t) based on the global coordinates of the vertices of the 
element and coordinates of the stress point. Regardless of the order of the 
element the following equations can be written: 


n 

X = I a. X, 
i=l 1 1 


n 

Y = I a. Y. 
. - 1 i 
i=l 


n 

Z = X a. Z. 

1 1 


n 

T = 1 a. T. 

i i 


where X, Y, Z 


T 


X. 

i 


Y. 

i 


Z. 

l 


N 

a. 


represents the geometry 
represents the normal field 

are the coordinates of the vertices of the heat 
transfer element 

is the number of vertices 

are the isoparametric shape functions 


6 . 4.2 


The problem is to solve for local coordinates (r, s, t) that correspond to the 
global coordinates of the stress point (Xp, Y^, Z^). The relationships between 
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the shape functions and the local coordinates are well known once the order of 
the element has been fixed. Substitution of these relations into the above 
summation along with the equating of X to X Y to Y and Z to Z yields a set 
of non-linear equations with r, s and t as Unknowns? The local^coordinates 
become functions of each other, global vertice coordinates and the coordinates 
of the stress point* 

That is: 


r = f(r , 
P P 

V 

t , 
p 

V 

V 

V 

X , 
p 

V 

V 


s = f(r , 
o p 

V 

t , 
p 

V 

V 

z i’ 

X , 
p 

V 

V 

6.4.3 

t = f (r , 
P P 

V 

t , 

p 

V 

V 

z i’ 

X , 

p 

V 

V 



These equations can now be solved numerically to yield values of r , s , 
t . Once the values of the local coordinates are known, a simple test Bete?- 
mSnes if the stress point is contained in the heat transfer element. The 
stress point is not contained in the element if any of the following conditions 
are met. 


-1 > r > 1 
P 

-1 > s > 1 

P 

-1 >t >1 
P 


If the point is found to be contained with local coordinates can be substituted 
into the expressions for the shape functions and the thermal weighting coef- 
ficient can then be computed. 

This is an ideal method for this interfacing procedure, it finds the cor- 
rect heat transfer element and at the same time produces accurate weighting 
coefficients. Figure 6-4 illustrates this mapping technique. 

For this study the 3D linear isoparametric formulation was chosen as the 
interpolating functions for the heat transfer elements. The set of non-linear 
equations described in Equation 6.4.3 were solved by using a fixed point ite- 
ration technique in conjunction with an M educated initial value M procedure. 
These initial values are computed as: 
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r . 
P* 


X - X „ 
p min 

X - x 
max min 


Y - Y . 

_ p mm 

~ Y - Y . 
max mm 


6.4.5 


t . 
pi 


Z - Z . 
p mm 

Z - z 
max min 


These initial values correspond to the correct values of r , s , t if the 
heat transfer element is a rectangular parallel piped. The ini?ial^guelses 
are used in conjunction with the righthand side of the non-linear equation to 
produce new values of the local coordinates. This iteration continues until 
convergence is achieved. The initial implementation of this technique showed 
it to be very efficient and fairly insensitive to element distortion. Several 
highly distorted 8-node 3D elements were tested. Figure 6-5 shows two of 
these elements. Convergence was achieved for all points attempted using these 
elements. Once this technique was installed into the transfer module and 
actual problems were interfaced, it was discovered that although the iterative 
procedure was not sensitive to the shape of the element, it was very sensitive 
to the orientation of the element with respect to the global axis system. It 
was found that for some configurations, the technique became unstable and 
either diverged or converged to an undesirable root of the equations. 

Logic was added to the transfer module to avoid these stability problems. 
First, a coordinate transformation was applied to the heat transfer element to 
align the direction defined by I and J joints of the element with the global X 
axis. The iterative technique is then used to solve for the local coordinates 
of the stress point. If it fails to converge to the correct root, the program 
reorients the element such that a new vertex of the element is now the I joint. 
The alignment transformation is then performed and the iterative solution is 
repeated. If it still fails to converge a new vertex is chosen for the I 
joint and the process is again repeated. If convergence is not achieved after 
all the vertices have been used, the stress point is not contained in this 
heat transfer element. 

A 400 element airfoil model, shown in Figure 6-6 was run using this ver- 
sion of the transfer module. The stress points used were the centroids of the 
elements and 8 additional points per element that were very near the corners 
of the element. There were a total of 3735 stress points. After the first 
transformation 3728 of the stress points converged. The other 7 stress points 
converged after the additional reorientation of the element was completed. 
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Figure 6-6. 400- 3D Element Air Foil Model Used As Test Case. 
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These results show that the initial transformation is adequate for most of the 
stress points and that the additional computation required for the vertex 
reorientation is only necessary for a very small percentage of stress points. 

As stated earlier the existing transfer module inverts the 3D linear iso- 
parametric shape function. The technique used, however, can be applied to any 
order element with only small impact on the coding of the transfer module. 


6.5 EXTERIOR POINT SURFACING ROUTINE 

In order for the transfer module to be a useful tool, it must deal with 
stress points that lie slightly outside of the heat transfer module. This 
situation occurs for many reasons in practice, the most common being part 
tolerances used to create the heat transfer model and the stress model. It is 
not the intention of the transfer module to produce temperatures for stress 
Points that are in areas outside of the domain of the thermal analysis and 
therefore no extrapolation of temperatures are performed. On the other hand, 
if a stress point exists that physically represents the surface of the part 
and that point because of slight tolerance mismatches is not contained in the 
heat transfer element on the surface, the transfer module should produce a 
temperature that is representative of the surface temperatures. This has been 
accomplished in TRANCITS by implementing an exterior surfacing routine. This 
technique also makes use of the local coordinates (r, s, t) . If the stress 
point lies outside the element, one of the local coordinates will be greater 
than one. For example, 

r = . 6 
s = -1.2 
t = -.8 

This point is not contained in a heat transfer element since 
s 1 

A point that lies on the surface near the stress point would have local 
coordinate values of 

r = . 6 
s = -1 
t = .8 

The distance between (r, s, t) and (r , s , t ) can easily be computed in phy- 
sical space and stored. The heat transfer element with the smallest distance 
can then be used as the closest element and the temperature of the surface 
point can be produced. This allows the transfer module to predict an approxi- 
mate temperature for the stress points that lie outside the heat transfer 
model. The program flags all stress points that exhibit this behavior and 
also prints out the distance between the actual coordinates of the stress 
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point and the point on the surface for which the temperatures were computed. 
Users should pay close attention to these points and judge if the distance 
outside the model is acceptable for the analysis being performed. Figure 6~7 
illustrates this surfacing technique along with the options of built in trans- 
formation and isoparametric mapping techniques. 
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Figure 6-7. Enhanced Program Features. 





7.0 VERIFICATION AND DOCUMENTATION OF THE TRANSFER MODULE 


Two objectives were met by the verification and documentation phase of 
this contract. The verification subtask provided the necessary test cases to 
ensure that all of the options implemented were producing accurate and effi- 
cient results. Indeed many efficiency improvements were identified and imple- 
mented in this phase. The purpose of the documentation subtasks was to pro- 
vide information required to use the transfer module as well as the documents 
needed to understand and possibly alter the coding for specific problems. 


7.1 VERIFICATION OF THE TRANSFER MODULE 

Verification of the transfer code was performed in parallel with the 
development effort. Each of the modules in the code were tested to ensure 
that the input and output of those modules were consistent with all file struc- 
tures. All user options were executed and thermal results from both SINDA and 
MARC were directly input and processed by the transfer code. Nodal tempera- 
tures for NAS TRAN were produced as well as elemental centroid temperatures for 
MARC . 


Formal verification of the mapped thermal results were conducted using 
three different 3D models. The primary objective of these test cases was to 
ensure that the thermal data produced by the heat transfer codes was not 
degraded in the mapping to the stress analysis models. The verification pro- 
cedure involved comparing the results of the heat transfer code to the inter- 
polated temperatures produced by TRANCITS. Comparisons were made using two 
different methods. The first method plotted the heat transfer results and the 
TRANCITS results against the model geometry for various regions of the 3D 
models. These plots were used to verify that the thermal gradients predicted 
by the heat transfer analysis were indeed mapped into the stress models. The 
second comparison method utilized the PATRAN Reference 3 postprocessing pro- 
gram. Colored isotherms were created for the stress model using the inter- 
polated temperatures and then for the heat transfer model using the original 
heat transfer results. These two sets of isotherms were then compared for a 
qualitative check on the thermal gradients. 

The first test case was a simple model of a rectangular prism. The heat 
transfer mesh was composed of a 2x2x5 grid, the stress mesh was much finer and 

had a mesh density of 4x4x10. These models are shown in Figure 7-1. The 

boundary conditions applied to the heat transfer model give rise to a linear 
thermal gradient along the length of the prism. The transfer module was used 
to map these temperatures into the finer stress model and Figure 7-2 shows a 

plot of the heat transfer code temperatures and the TRANCITS temperatures 

versus the distance along the prism. The figure shows that for this simple 
case the interpolated values agree exactly with those predicted by the heat 
transfer code. Figure 7-3 and 7-4 show the isotherms for the heat transfer 
model and stress model respectively. Visual comparison of these isotherms 
again indicate the same thermal gradients in both models. 
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Figure 7-1* Prism Test Case 
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The model used for the second test case is shown in Figure 7-5. It is 
finite difference model of the tip of a turbine blade composed of approxi- 
mately 450 heat transfer elements. For this case the stress model mesh den- 
sity was identical to the heat transfer model. The transfer module was used 
to map the heat transfer temperature at the centroid of the elements and the 
face centers of the element sides to the vertices of the element. The bound- 
ary condition applied to the heat transfer model represent realistic condi- 
tions for a turbine blade and give rise to gradients in all directions. Fig- 
ures 7-6 thru 7-11 show plots of the original heat transfer temperatures and the 
transferred corner temperatures against a local distance parameter. These plots 
are for 3 spanwise sections through the airfoil (Tip, Middle, Root Sections) 
for both the concave and convex sides of the blade. The local distance param- 
eter is measured from the leading edge on the concave side and the trailing 
edge on the convex side. The first two curves on each plot represent the 
temperatures for a layer of heat transfer elements just above and just below 
the corner node temperatures. As the plots show, the corner temperature dis- 
tribution not only has the same trends as the original heat transfer tempera- 
ture distribution but the value of the temperature at each corner node falls 
in between the temperature above and below the corner as computed by the heat 
transfer code. 

These comparisons provide excellent verification of the accuracy of the 
transferring technique. 

Models representing the final test case are shown in Figure 7-12. These 
models represent a 3D sector of a combustor liner. The mesh density for the 
heat transfer model is finer than the stress model in some locations and 
coarser in others. The boundary conditions used in the heat transfer model 
once again represent typical engine conditions. Figures 7-13 and 7-14 show 
isotherms for both heat transfer and stress grid respectively. These isotherms 
also verify that the gradients are essentially the same for both models. 


7.2 DOCUMENTATION OF THE TRANSFER 


The documentation of this code was divided into two segments. The first 
was the TRANCITS User's Manual. This manual describes in detail all of the 
input and output file contents for the options implemented in this version of 
the transfer module. It shows examples of these files and discusses how the 
code should be executed. It also has a complete description of all the user 
input variables and what defaults if any are used if these inputs are omitted. 
The manual specifies the current limitations of the code in regard to maximum 
size of problems that can be interfaced. The manual should provide all of the 
necessary information required to use the transfer modules. 

The second involved the documentation of the coding of the transfer mod- 
ule. Listings and prologs of each subroutine used and a flow diagram of the 
entire code have been included in this report in Appendices C and D. These 
prologs are part of the source code and they give a brief description of the 
purpose of the subroutine. They also describe all of the calling arguments 
and whether these arguments are inputs or outputs to the routine. All files 
used in the routines are specified as well as any routines that are called 
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Figure 7-5. Airfoil Test Case Model. 
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Figure 7-7. Temperature Verificatio 
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Figure 7-9. Temperature Verif icati 
Large Model (Airfoil) . 












Figure 7-10. Temperature Verification Plots Concave Side Root Section 
Large Model (Airfoil) . 
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from the subroutine. Commons that are used in the routine are also noted. If 
there are any special limitations or comments, these are described in the 
prolog. 

The flow diagram gives a concise picture of the program flow and provides 
a convenient reference for which subroutines are a part of each module. This 
and the program documentation should enable any user to understand and perhaps 
adapt this transfer module to their individual needs. 
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8.0 CONCLUSIONS 


The TRANCITS software package accurately and efficiently transfers thermal 
data from a 3D heat transfer study to a 3D stress analysis mesh of different 
elemental density. The heat transfer analysis can be performed using finite 
element or finite difference techniques. TRANCITS currently supports the 
output of the MARC and SINDA heat transfer codes directly and will format the 
thermal data compatible with the input requirements of NASTRAN and MARC. In 
addition, the neutral heat transfer input file and neutral temperature output 
file can be used to interface the results of any heat transfer code to any 
structural code. The transfer module has been written in a modular fashion 
and is easy to implement and modify. The technology developed to transfer 
thermal data forms an excellent foundation for the transfer of other engi- 
neering parameters such as pressures, loads and displacements. 
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APPENDIX A 


Survey of Computer Programs for Heat Transfer 
(A.K. Noor) 


The following tables and discussions were extracted from A.K. Noor's paper 
"Survey of Computer Programs for Heat Transfer Analysis." 

These tables illustrate the solution capabilities, the heat transfer ele- 
ments supported, the pre and post-processing options, the type of boundary 
conditions supported and any interfacing options for six of the heat transfer 
codes considered. 
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Pari I « Analytical Capabilities of the Progra* 



♦Extracted form "Survey of Computer Frogmen 


ir Heat Analysis” By A.K. Hoor 


























Solids of devolution 
Boundary UnefiU 
Scaler Clements 
Point Contact dements 
Others (see program abstracts) 
Range of Applications and Phenomena 
linear Steady State 
Oonllnear Steady State 
Thermal frequencies and ttode Shapes 
linear Transient Response 
flonllnear Transient Response 
Others (see program abstracts) 
formulation 

a) Fundamental Unknowns 

Temperatures 
Temperatures and flux 
Enthalpy 

Others (see program abstracts) 

b) Elemental Matrices 

Conduction 


ANSYS 



























ert (tie Multipoint Clmslraffils} 
Iheiial Properties 























Latent Heat and Phase Change effects 
Material Added or R e m oved During Analysis 
User Supplied (see prog ram abstracts) 

Others (see program abstracts) 

Initial Conditions 
Homogeneous. 

Varying Throughout the Region 
Initial tnthalpy (for Phase Change) 

User* Supplied 

Others (see program abstracts) 
boundary Conditions and thermal toads 

a) Steady State 

b) time Dependent 

a) Steady 

b) Te»pprit»re Dependent 

c) time Varying 

a) Steady State 

b) time Dependent 


Prescribed temperatures 
thermal flux Input 

Convection from a Surface 
to Its Surroundings 


ANSYS 
































oi 

00 





ANSYS 

o 

ss 

X 

| 

43 

S3 

VONIS 

| 

e 



Convection from Surroundings 
to a Surface 

a) Steady State 

• 

a 

a 

a 

a 



h) lime Dependent 

• 

a 

a 

a 

a 

a 


forced Convedfoti 


• 


a 

a 

a 



Radiation frm a Surface 
to Its Surroundings 

a) Steady State 

• 

a 

a 

a 

a 



b) flue Dependent 

a 

a 

a 

a 

a 

a 


Radiation fro* Surroundings 
to a Surface 

i.Jl.il^. at.. ao. - . . _ 

a) Steady State • 

a 

a 

a 

a 

a 



b) flue Dependent 

• 

a 

a 

a 

a 

a 


Muiation Dfuren wirrw ups 

a 

a 


a 

a 



Radiation Betneen 
n Surfaces with 

a) User* Supplied 
flew factors 

a 

a 


a 

a 

a 


1) Internal ty 
Calculated 
flea factors 




a 




Prescribed field Hot# 

a 



a 

a 



Boundary layer Convection 



a 

a 




Votwetrlc Heat Generation 

a) On Clement level 

• 


a 



a 


b) 0n Node level 


a 

a 

a 

a 



Sap (Thenaat Resistance) 




_a__._. 

a 

- a . 

a 


Boundary Condi t lons/loads Added or leiiived Durfny Analysis 

• 


a 

a 

a 



Olliers (see program abstracts) 




































a 

o 



Iciferilwt Field ted Transmitted Directly trow Neat 
transfer Nodules to t hemal Stress Nodules 

enclosure Radiation with View Factor Calculation 

Heat Input /Out put at Cons Ira Ined Boundaries 

Cyclic Symmetry 

Substructure 

a) Repeated Use of Identical Substructures 

b) Nixing linear and Nonlinear Substructures 

c) Hiring Substructures with Different 
lypes of Nonttnearltles 

Restart Capability 

Others (see program abstracts} 

a 

Program Operational On 


COC 

I BN 

UNIVAC 

Honeywell 

felefunten 

AMDAHL 

litmus 























































k) Aetoaetlc or Seial-Autoaetlc Genera tor^for^ 

0ne>0laenslona1 tl wwtl 

trtanpnlar Elements 

Quadrilateral EleMents 

Body or Shell of Revolotl on tleaents 

Three- Diarm tonal Solid llwwH 

Tno*Dlaenslonat Sh ell Elements 
Curvilinear Finite Differ ence Grids 
Others I see prograa ab stracts) ^ 

c) Data Checking Facilities 

tine Frlnter 

Hotter _______ 

Interactive Graphics 

Others (see prograa abstracts) 

d) Flots and Gra phics Display of Hodel ^ 

Coaplete Analysts R egion 

Fart of Analysis Region 

•llow-Mp* Option 


o 

CO 


AN SYS 






Hidden Lines or Surfaces 
Orthographic Views 
Perspective and Isometric flews 
Section View on Arbitrary Plane 
Others (see program abstracts) 
e) Other Tael I Hies 
Digitizer Input 

Automatic Renumbering of Nodes, Clements or equations 

Table lookup of Oata 

Others (see program abstracts) 

Results Output Fonu 
a) Tabular Output 
Fixed Set 

User Defined Set and Sequences 
Maximum and Minimum Quantities 
Average and Maxima for Slocks of Modes 
temperature or Flux Exceedances 
Olliers (see program abstracts) 


























o 

u% 

£ 

5 

o 

53 

JL 

MARC 

£3 

53 H 

VflNIS 

1 

H 



b} file Output far User Poil-rmMilwi and Plotting 

9 

9 

9 

9 

o 



c) Plots 








Isotherm Plot* (Contour*) of Tim^riturw/Fliw 

o 

9 

O 

9 


o 


Surface function* 






9 


Selective Output (e.e.* by i learnt* or Pee Ions) 

o 

O 

O 

O 




Histories (#.§., tlw History) 

o 

O 

o 

' o 

O 

• 


Other* (see program abstract*) 




# 




4. 

Interactive Input and Control 









Parameter Specie!! feat ton (e.g. 9 flux or 11 we Step*) 

o 





• 

j 

Singularity Check 



9 





Crror Correct f on/Pecovery 

o 







lf*er Control of Matrix Decomposition 








Others (see program abstracts) 
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APPENDIX B 


Survey of Structural Analysis Codes 


The following tables were extracted from a survey performed by R. Zirin of 
General Electric Gas Turbine Division. 

These tables illustrate the solution capabilities, the element libraries, the 
types of loadings and the output features of the following structural codes: 
EASE2 , STARDYNE, NASTRAN, ANSYS, MARC, SUPERB. 
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TYPES OF ANALYSIS 

ft 

III 

• 

| ANALYTICAL CAPABILITY 

2 

UNEAR 
STATICS * 

MECHANICAL LOAOS 



TEMPERATURE LOADS 



EULER BUCKLING 



INERTIA RELIEF 

r 


DYNAMICS 

MODE/FREQUENCY 

ID 


FREQUENCY RESPONSE 

r 


TRANSIENT RESPONSE 

ID 


SHOCK SPECTRA 

D 


RANDOM RESPONSE 



NONLINEAR TRANSIENT 



NONLINEAR 

STATICS 

NONLINEAR BUCKLING 

■ 

1 

LARGE DISPLACEMENT 


1 

PLASTICITY 

D 

1 

CREEP 



VISCOELASTICITY 

■ 

1 

LARGE STRAINS 

□ 


HEAT TRANSFER 

STEADY STATE 

■ 


TRANSIENT 

■ 


SUBSTRUCTURES 

(SUPER- 

ELEMENTS) 

STATIC 

D 


OYNAMIC 

D 

1 

CYCLIC SYMMETRY 

n 



FRACTURE MECHANICS 

■ 



FLUIDS 

D 


litefti t iftiefitK 

ELECTRIC CIRCUITS 




OPTIMIZATION 




ACOUSTIC CAVITIES 




FATIGUE DAMAGE 

■1 



















































STRUCTURAL ANALYSIS 
ELEMENT/MATRIX LIBRARY 



WOTCS: 


M Mamdrana and/or plant strain only (no piatt btndtnp) 

S Inciudai suoparamatnc (ormt witn ftwtr nodts 

C Alto >ndvd«t cubic itooaramatr* altmtnt wim two 
m«ds*dc nodat 
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ANSYS 






































STRUCTURAL ANALYSIS 
ELEMENT/MATRIX LIBRARY (continued) 

| .PROGRAM 

ee 

STAROYNE 

5 

1 

Z 

ANSYS 

MARC 

IB 

S 

i 

| . CLIME NT 

III 

AX!* 

SYMMETRIC 

ELEMENTS 

SHELLS 

. CONICAL 

\ 



B 

B 

B 

a 

CURVED 

A 



fl 

B 

B 

B 

TRIANGULAR 

RINGS 

3 NODE 

A 




B 

B 

fl 

6 NODE 

A 



B 

B 


B 

QUAD 

RINGS 

4 NODE 

□ 



fl 

B 

B 

B 

8 NODE 

□ 



B 

B 

B 

B 

SOLIO 

ELEMENTS 

TETRA- 

HEDRON 

4 NODE 




B 

B 

B 

B 

WEDGES 

6 NOOE 




B 

B 

B 

B 

15 NOOE 

<3 



B 


B 

fl 

HEXA- 

HEDRONS 

8 NOOE 



! 

i 

B 

B 

B 

B 

20 NOOE 

a 



B 

B 

B 

B 

PIPE ELEMENTS 

STRAIGHT 




fl 

fl 

g 

B 

ELBOW 

r 



B 

B 

B 

B 

TEE 

•MjM* 



B 

fl 


B 


NOTES; 

S Indudes aubperametric forms with fnwf nodes 
C Also indudes cubic isoparametric element widi too 
midside nodes 
0 Degenerate case 
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STRUCTURAL ANALYSIS 

ELEMENT /MATRIX LIBRARY [continued) 

| PROGRAM | 

w 

m 

< 

m 

tu 

Z 

o 

c 

? 

¥» 

NASTRAN 

ANSVS 

MARC 

SUPERB 

| - ELEMENT 

GENERAL 
STIFFNESS * 
ELEMENTS 

STRING 


t 

9 

r_ 

• 

• 

• 

SCALAR SPRING 




• 




8 * 6 or 12 * 12 MATRIX 



9 

• 

• 


• 

GENERAL MATRIX 



• i 


MASSES 

ELEMENT 

LUMPED (DIAGONAL) 


z 

Z 

| z 

Z 


CONSISTENT 




z 

* 

Z 

NON* 

STRUCTURAL 

SCALAR (OOF) 


9 


• 


NODAL 


9 9 

9 

9 


DISTRIBUTED 


1 

1 

• 




GUYAN REDUCTION 


1 -• 

t 

• • 

1 k 

«, 8 

1 


GENERAL MATRIX 




9 9 

i ! 
! 1 

0 AMP! BIG 

SCALAR 




9 




OASHPOT 

o 



9 

9 



OISCRETE VISCOUS (C) --(Kj * A(M| 


• 


9 | 

* i 

• 


STRUCTURAL ft «* ig)(Kj 




• 

1 

i 



MODAL VISCOUS 


9 


• 


m 


GENERAL MATRIX 



1 • i • i 

| 

OTHER ELEMENTS 

GAP 

HI- 




• j 

• 


FRICTION 

m. 




• 

• 1 


RIGIO 


* 

• 

• 




REBAR SOLIO 








ELASTIC FOUNDATION 





• 1 

• 


CRACK TIP 





9 

# 


•LAMINATED SHELL 

C 9 




9 

• 

*5 

PLOT ONLY * 

i 



• 1 i 

- 


NOTES: 

f S«« mtrtino 
2 G*rmr* tat from tftmity 
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PROGRAM 


HEAT TRANSFER-CONDUCTING ELEMENTS 


1 ELEMENT | 

LINEAR 

/ 

PLANAR 

3 NODE TRIANGLE 

A 

4 NODE QUAD 

n 

e node quad 

□ 

transverse conducting 

SHELL 

d 

AXISYMMETRIC 

TRIANGULAR RING 

A 

4 NODE QUAD RING 

□ 

8 NODE QUAD RING 

□ 

a 

SOLID 

TETRAHEORON 


WEDGE 

V 

8 NODE BRICK 

0 

IS NOOE WEDGE 

$ 

30 NODE 8RICK 

0 

GENERAL MATRIX INPUT 

■1 



< 

z 




NOTES: 

S Contain* aubparamatrie form* with fawar number *of nodat 
P Alto contain* parabolic isoparametric aiamant with ona midtid* nod* 
C Abo contain* cubic isoparametric aiamant with two midsida node* 

0 Dapanar a t a cm 
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8U3<mS 











































COORDINATE SYSTEMS 
AND MATERIAL PROPERTIES 

PROGRAM 

« 

tu 

< 

m 

STAROYNE 

NASTRAN 

ANSYS 

MARC 

SUPERB 

| ~ FEATURE 


■ASIC 

(GLOBAL) 

CARTESIAN 

# 

m 

• 

! # 

• 

• 

• cylindrical 

• 

• 

• 

• 

• 

SPHERICAL 


• 

• 

• 

GENERAL 


1 i 

COORDINATE 

SYSTEMS 

SXEWED 

(LOCAL) 

CARTESIAN 

• 

I_l 

• i • 

• 

• 

CYLINDRICAL 

• 

i • 

• • 

• 

SPHERICAL 



• 

• 

• 

GENERAL 





% 



MIXED 

• | • 

• 

• • | • 

MATERIAL PROPERTIES 

ISOTROPIC 

• 1 • 

• • 

• 

• 

2-0 ORTHOTROPIC 


• 

• • 

1 

• 

3-0 ORTHOTROPIC 



• 

1 

m 

TEMPERATURE DEPENDENT 

m 


• m m 

• 

STRESS DEPENDENT 

. 1 . . 

i 


TIME DEPENDENT 


• ! • i 



NONLINEAR ELASTIC 


! • 


WORK 

HAROENING 

ISOTROPIC 




KINEMATIC 



* 

i • 

i 

• ] 


COMBINED 


1 • 

• 


ORNL 10 CYCLE 



• 

• 


GENERAL 



i 

Z 


NOTES; 

t Pcrfomwd by uur wbrautin* 
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BOUNDARY CONDITIONS 


FEATURE 


FRO CRAM 


X 

I 2 « 

1 SI 




CONCENTRATED 

IDDBDB 



DISTRIBUTED (BEAM) 

DDDDD 




PLATES/SHELLS 

Donna 



FRESSURE 

AXISYMMETRIC ELEMENTS 

nnnnn 




SOLIOS 

nnnnn 


STATIC 

TEMPERATURE 

DDDDD 



ACCELERATION 

nnnnni 

LOAOING 


ROTATIONAL VELOCITY 

nnnnni 



COMBINATION 

DDDDDI 



AXI* 

AXISYMMETRIC SHELLS 

■■nnei 



SYMMETRIC 

AXISYMMETRIC RINGS 

• 



TIME OEPENOENT 

nnnnn 


OYNAMIC 

FREQUENCY DEPENDENT 

MHIM 


PSO RANOOM 

nnnnn 

/ 


SHOCK SPECTRUM 

nnnnn 



SINGLE POINT* 

M— » 

DISPLACEMENT 

CONSTRAINTS 

MULTI POINT* 

nnnnn 



SPECIFIED NONZERO DISPLACEMENT 

nnnnni 



HEAT SOURCE/SINK 

—BHH 

HEAT 

CONVENTION 

nnnnn 

TRANSFER 

RA0IAT10N 

nnnnn 

L_ 


SPECIFIED TEMPERATURE 

nnnnn 


NOTES; *$2ngie point constraint s enforced tare trantiation(s) and/or rotation(s) in coordinates} associated with • node point 
Multi-point constraint • enforced Sneer constraint relationships between translation (s I end/or rotation(s) which may ^be 
associated with different node points 

1 Applies to some elements 3 Displacement components set equal on different nodes 

2 Specialised forms of rigid and interface coupling ♦ Stand alone program 
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SUPERB 








































PRE- AND POST-PROCESSING 



FEATURE 


UNOEFORMEO GEOMETRY 


NODE LABELS 


ELEMENT LABELS 


PROPERTY LABELS 


2*0 SECTIONS 


BOUNOARY CONDITION LA8ELS 


HIOOEN LINES REMOVED 


DEFORMED GEOMETRY 


CONTOURS 20 STRUCTURE 


CONTOURS SOLIO STRUCTURE 


TIME HISTORY 


REQUENCY RESPONSE 


er" 


ARBITRARY. X VS. Y 


NODES 


RESTRAINTS 


LOAOS 


BY LOAD CASES 


BY ELEMENT 


SORTING 

MAX/MIN SUMMARY 

BBMMB 


SELECTEO NOOES AN 0/0 R ELEMENTS 

• 

BBMBB 

IANOWIOTH MINIMIZATION 

bbhheeieq 



NOTES; 


1 Ganarctat (tote In I "dlmamiorT 

2 Ganarita* data in 2 “dimansaera** 

3 data in 3 "dimerwora" 

4 Printer plots 

* Stand dent program 
W Wavefront totutien 
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SUPERB 


















































APPENDIX C 


Subroutine Prolog 

This appendix lists the prologs for all of the subroutines used in the 
transfer module. The prologs represent a standard description of the purpose 
of the routine, the calling arguments and their attributes, the routines that 
are called, any commons that are used and special comments about the sub- 
routine. 


75 





* DIRECTORY OF SUBROUTINES/FUNCTIONS - SORTED ALPHABETICALLY BY ROUTINE NAME * 
********#**************************«************************-************(******** 


FILE NAME 


GE/TRANCITS 
GE/TRANCITS 
GE/TRANCITS 
GE/TRANCITS 
GE/TRANCITS 
GE/TRANCITS 
GE/TRANCITS 
GE/TRANCITS 
GE/TRANCITS 
GE/TRANCITS 
GE/TRANCITS 
GE/TRANCITS 
GE/TRANCITS 
GE/TRANCITS 
GE/TRANCITS 
GE/TRANCITS 
GE/TRANCITS 
GE/TRANCITS 
GE/TRANCITS 
GE/TRANCITS 
GE/TRANCITS 
GE/TRANCITS 
GE/TRANCITS 
GE/TRANCITS 
GE/TRANCITS 
GE/TRANCITS 
GE/TRANCITS 
GE/TRANCITS . 
GE/TRANCITS 
GE/TRANCITS 
GE/TRANCITS 
GE/TRANCITS 
GE/TRANCITS 
GE/TRANCITS 
GE/TRANCITS 
GE/TRANCITS 
GE/TRANCITS 
GE/TRANCITS . 
GE/TRANCITS . 
GE/TRANCITS . 
GE/TRANCITS 
GE/TRANCITS 
GE/TRANCITS 
GE/TRANCITS 
GE/TRANCITS 
GE/TRANCITS . 
GE/TRANCITS 
GE/TRANCITS . 
GE/TRANCITS 
76 GE/TRANCITS . 

GE/TRANCITS . 


SUBROUTINE/ 

PAGE NO 

FUNCTION 


AREAD 

17 

CALCRD 

12 

CCNVEC 

22 

CHKPT 

92 

CHNODT 

7 

CINV 

37 

CKWIND 

53 

CORNCC 

21 

CTRAN 

9 

CWCRNT 

31 

DATINT 

75 

EINCOR 

49 

ELTRAN 

98 

EMINMX 

40 

ERRPRT 

84 

ETCORE 

64 

FCRNTP 

34 

FILREW 

48 

FNDELM 

50 

FNDJAC 

94 

FTIME 

41 

GENTMP 

61 

GETCRD 

23 

GETEWT 

77 

GETMPE 

26 

GETPD 

13 

GETWNO 

52 

GNLSUR 

46 

GPCOEF 

101 

GPTEMP 

102 

GTGAU5 

19 

GTNUMT 

78 

HTFIO 

58 

HTFPAS 

67 

HTICON 

69 

I3DSF 

87 

INITCM 

99 

JACBCK 

96 

KNTSPF 

43 

MARCRD 

70 

MARCTO 

100 

MAXTDF 

35 

MKRGF 

54 

MPERW 

27 

NTCORE 

65 

PARD 

16 

PDOFSF 

93 

. PRCGEO 

79 

PRCSOF 

76 

PREPGM 

8 

PRNTIM 

39 


GE/TRANCITS 

PROFAC 

56 

GE/TRANCITS 

PUTWFI 

24 

GE/TRANCITS 

RDCNTL 

59 

GE/TRANCITS 

RDSPF 

47 

GE/TRANCITS 

RDTEMP 

30 

GE/TRANCITS 

RDTIME 

4 

GE/TRANCITS 

REAARY 

89 

GE/TRANCITS 

REACRD 

20 

GE/TRANCITS 

READR 

71 

GE/TRANCITS 

REMELM 

11 

GE/TRANCITS 

RFILCR 

44 

GE/TRANCITS 

RWCCF 

25 

GE/TRANCITS . 

RWELMS 

2 

GE/TRANCITS 

RWNODS 

1 

GE/TRANCITS 

RWSPCF 

97 

GE/TRANCITS 

SCTION 

10 

GE/TRANCITS 

SETFLP 

29 

GE/TRANCITS 

SFDIST 

91 

GE/TRANCITS 

SHAPEF 

15 

GE/TRANCITS 

SHAPFV 

95 

GE/TRANCITS 

SINTIT 

73 

GE/TRANCITS 

SKIP 

57 

GE/TRANCITS 

SKPREC 

72 

GE/TRANCITS 

SORT 

68 

GE/TRANCITS . 

S0RT2 

90 

GE/TRANCITS 

SSURCH 

42 

GE/TRANCITS 

STSTMP 

62 

GE/TRANCITS 

TEMFOR 

85 

GE/TRANCITS 

TFLSET 

60 

GE/TRANCITS . . 

TIMCK 

6 

GE/TRANCITS 

TPNAST 

86 

GE/TRANCITS 

UNPAKF 

33 

GE/TRANCITS 

VOLM 

14 

GE/TRANCITS 

WRCRNT 

66 

GE/TRANCITS 

WRITMP 

36 

GE/TRANCITS 

WRSTMP 

63 

GE/TRANCITS 

WRTDIR 

81 

GE/TRANCITS 

WRTEMP 

83 

GE/TRANCITS 

WRTGEO 

82 

GE/TRANCITS 

WRTNTF 

80 

GE/TRANCITS 

ZERORF 

45 


77 


TITLE * 


DATE =07/15/85 


TIME = 23 42 


PAGE= 1 


2640 SUBROUTINE RWNODS( IREC.NN.XYZ, IRW) 

2670CPR 

2680CPR 

2690CPR ROUTINE-- RWNODS 
2700CPR 

2710CPR DATE/WRITTEN BY— 8/04/82 RU MAFFEO 

2720CPR DATE/REVISED BY — 

2730CPR 

2740CPR FUNCTION/PURPOSE-- 

2750CPR THIS I/O ROUTINE READS AND WRITES THE NODAL INFO 
2760CPR TO THE INTERNAL DATA BASE FILES 
2770CPR THIS FILE IS NOW IFL3 

2780CPR IN THE WRITE MODE IT AUTOMATICALLY ZEROS THIS FILE 
2790CPR FOR NODES NAMES THAT ARE NOT USED 

2800CPR THIS RANDOM FI LE ( I FL3 ) COULD EASILY BE REPLACED BY 

2810CPR AN ARRAY IF ENOUGH CORE WAS AVAI ALABEL 

2820CPR 

2830CPR CALLING SEQUENCE-- CALL RWNODS ( I REC , NN , XYZ , I RW ) 
2840CPR 

2850CPR CALLING ARGUMENTS — 

2860CPR NAME ATTRIBUTES DEFINITION 


2870CPR 

IREC 

(i) 

RECORD TO BE PROCESSED 

2880CPR 

NN 

(I/O) 

NODE NUMBER (NAME) 

2890CPR 

XYZ 

(I/O) (3X1) 

ARRAY OF COORDINATES 

2900CPR 

IRW 

(I) 

READ/WRITE SWITCH 


2910CPR I RW*0 READ 

2920CPR I RW= 1 WRITE 

2930CPR 

2940CPR 

2950CPR 

2960CPR FILES USED — 

2970CPR I FL3 RANDOM FILE TO STORE COORDINATES 

2980CPR 

2990CPR 

3000CPR COMMONS USED — 

3010CPR AFIL FILE CODE COMMON 

3020CPR 

3030CPR 

3040CPR FUNCTIONS/ROUTINES CALLED — 

3050CPR NONE 

3060CPR 

3070CPR LIBRARIES ACCESSED — 

3080CPR NONE 

3090CPR 

3100CPR LOCAL VARIABLES — 

3110CPR N/A 

3120CPR 

3130CPR SPECIAL COMMENTS ABOUT THIS ROUTINE 

3140CPR NONE 

3150CPR 

3160CPR 

3170CPR SPECIAL REMARKS/ INST RUCTIONS — 

3180CPR HOST = H6000 

3190CPR 

3200CPR ********************************** 

3630 END 
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GE/TRANCITS 


RWNODS 


PAGE 5 1 



TITLE* 


DATE *07/15/85 


TIME * 23 42 


PAGE* 2 


3650C** 

3670 

36900 * * 

3700CPR 

3710CPR 

3720CPR 

3730CPR 

3740CPR 

3750CPR 

3760CPR 

3770CPR 

3780CPR 

3790CPR 

3800CPR 

3810CPR 

3820CPR 

3830CPR 

3840CPR 

3850CPR 

3860CPR 

3870CPR 

3880CPR 

3890CPR 

3900CPR 

3910CPR 

3920CPR 

3930CPR 

3940CPR 

3950CPR 

3960CPR 

3970CPR 

3980CPR 

3990CPR 

4000CPR 

4010CPR 

4020CPR 

4030CPR 

4040CPR 

4050CPR 

4060CPR 

40‘70CPR 

4090CPR 

4100CPR 

41 10CPR 

4120CPR 

4130CPR 

4140CPR 

4150CPR 

4160CPR 

4170CPR 

4180CPR 

4190CPR 

4200CPR 

4210CPR 

4220CPR 

4230CPR 

4240CPR 

4250CPR 

4260CPR 

4270CPR 

4280CPR 

4290CPR 

4300CPR 

4310CPR 

4320CPR 

4330CPR 

4340CPR 


SUBROUTINE RWELMS( I REC , A , IRW , ICODE ) 


ROUTINE— RWELMS 

DATE/WRITTEN BY-- 8/04/82 RJ MAFFEO 
DATE/REVISED BY- 

FUNCTION/ PUR POSE— 

THIS ROUTINE IS AN I/O ROUTINE 

ITS PURPOSE IS TO READ/WRITE THE ELEMENT ATTRIBUTES TO THE 
INTERNAL DATA BASE FILE 

THIS DATA BASE FILE IS RANDOM FILE IFL4 
THE ATTRIBUTES ARE 

— ELEMENT NAME 

-- COORDINATES OF ALL 15 POINTS IN THE ELM 
-- ELEMENT CONNECTIVITY OF THE 8 CORNERS 
-- THE PACKED THTD FACE NUMBER 
-- VOLUME OF ELEMENT 

THE ROUTINE CAN PROCESS DIFFERENT AMOUNTS OF THESE 

ATTRIBUTES BASED ON THE VALUE OF ICODE 

THE RANDOM FILE IS NOT INITIALIZED PRESENTLY WITH THIS 

ROUTINE 


CALLING SEQUENCE-- CALL RWELMS ( IREC , A , I RW f ICODE ) 

CALLING ARGUMENTS — 

NAME ATTRIBUTES DEFINITION 

IREC (I) RECORD NUMBER TO BE PROCESSED 

A (I/O) (56X1) ARRAY OF ATTRIBUTES 

IRW (I) READ/WRITE SWITCH 

I RW*0 READ 

IRW* 1 WRITE 

ICODE (I) SWITCH TO SET AMOUNT OF DATA TO PROCESS 

IC0DE*0 — PROCESS ALL DATA 
ICODE* 1 -- PROCESS ELM NAME ONLY 

IC0DE*2 -- PROCESS ELM NAME AND CORNER COORDINATES 
IC0DE*4 -- PROCESS EVERYTHING EXCEPT VOLUME 
IC0DE*5 -- PROCESS ALL DATA 


FILES USED — 

I FL4 -- RANDOM FILE TO STORE ELM ATTRIBUTES 


COMMONS USED — 

AFIL -- COMMON TO DEFINE FILECODES 


FUNCTIONS/ROUTINES CALLED — 
NONE 


LIBRARIES ACCESSED — 

NONE 

LOCAL VARIABLES — 

N/A 

SPECIAL COMMENTS ABOUT THIS ROUTINE 
NONE 
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4350CPR 

4360CPR 

4370CPR SPECIAL REMARKS/ INSTRUCT IONS-- 

4380CPR HOST = HSOOO 

4390CPR 

4400CPR ************************************ 

4850 END 
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4870C = = 

4890 

4920C = 3 

4930CPR 

4940CPR 

4950CPR 

4960CPR 

4 9700 PR 

4980CPR 

4990CPR 

5000CPR 

5010CPR 

5020CPR 

5030CPR 

5040CPR 

5050CPR 

5060CPR 

5070CPR 

5080CPR 

5090CPR 

5100CPR 

51 10CPR 

5120CPR 

5130CPR 

5140CPR 

5150CPR 

5160CPR 

5170CPR 

5180CPR 

5190CPR 

5200CPR 

5210CPR 

5220CPR 

5230CPR 

5240CPR 

5250CPR 

5260CPR 

5270CPR 

5280CPR 

5290CPR 

5300CPR 

5310CPR 

5320CPR 

5330CPR 

5340CPR 

5350CPR 

5360CPR 

5370CPR 

5380CPR 

5390CPR 

5400CPR 

5410CPR 

5420CPR 

5430CPR 

5440CPR 

5450CPR 

5460CPR 

5470CPR 

5480CPR 

5490CPR 

5500CPR 

5510CPR 

5520CPR 

5530CPR 

5540CPR 

5550CPR 

556O0PR 


SUBROUTINE RDTIME 


ROUTINE-- RDTIME 

DATE/WRITTEN BY-- 8/18/82 RJ MAFFEO 

DATE/REVISED BY— 9/13/83 SR MARTIN 

REASON — CONVERT FROM F4BT FILES TO THE NEW 

HEAT XFER OUTPUT FILE FORM 

FUNCTION/PURPOSE — 

THIS ROUTINE READS THE TIME VALUES FROM THE HEAT XFER OUTPUT FILE 
AND DETERMINES IF THEY ARE DESIRED FOR THIS RUN IT 
AUTOMATICALLY ELIMINATES THE NON-STP TIME STEPS IT ALSO 
CALLS THE REQUIRED ROUTINE TO PERFORM INITIAL CHECKS ON 
TIMES AND ELM COUNTS 

THE HEAT TRANSFR OUTPUT FILE IS REWOUND BEFORE RETURNING 
FROM THIS ROUTINE 


CALLING SEQUENCE— CALL RDTIME 

CALLING ARGUMENTS — 

NAME ATTRIBUTES DEFINITION 

FILES USED — 

INTITS — HEAT XFER OUTPUT FILE (TITAN INPUT) 
I ERF -- ERROR REMARK FILE 
I NOT 6 — HARD COPY FILE 


COMMONS USED — 

AFIL — FILE CODES 

SIZE -- PROBLEM SIZE PARAMETERS 

CNTLTM -- TIME/TEMPERATURE CONTROL PARAMETERS 


FUNCTIONS/ROUTINES CALLED — 

TIMCK 

CHNODT 


LIBRARIES ACCESSED — 
NONE 


LOCAL VARIABLES — 

NSTPTM — COUNT OF STP TIMES PROCESSED 

NCNTT -- COUNT OF TOTAL TIMES PROCESSED 

NTMDO — USER SPECIFIED NUMBER OF STP TIMES TO DO 

TIMT — ARRAY OF TIMES PROCESSED THIS ARRAY WILL 

-- BE USED AS A DIRECTORY TO POINT TO THE 
— LOCATION ON THE HEAT XFER OUTPUT FILE FOR EACH TIME 


SPECIAL COMMENTS ABOUT THIS ROUTINE 


SPECIAL REMARKS/INSTRUCTIONS 
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5570CPR HOST = H6000 

5580CPR 

5590CPR *************************** 

6640 END 
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6660C SSSSSSSSCSSSSSSKeSSSSSSSE=SSSKSS5SSSSSSSSSS==SSSSS=SSS=£SSSS= 

6680 SUBROUTINE TIMCK( TIME , ITIME ) 


67 10C SSS3SSSSSSSSSSSSSS ss = = = S = SSS = SSS&XSSeaS£SSS = 3E = SSSSE = ESXSSS = a 


6720CPR 
6730CPR 
6740CPR 
6750CPR 
6760CPR 
6770CPR 
6780CPR 
6790CPR 
6800CPR 
6810CPR 
6820CPR 
6830CPR 
6840CPR 
6850CPR 
6860CPR 
6870CPR 
6880CPR 
6890CPR 
6900CPR 
6910CPR 
6920CPR 
6930CPR 
6940CPR 
6950CPR 
6960CPR 
6970CPR 
6980CPR 
6990CPR 
70OOCPR 
7010CPR 
7020CPR 
7030CPR 
7040CPR 
7050CPR 
7060CPR 
7070CPR 
7080CPR 
7090CPR 
7100CPR 
71 10CPR 
7120CPR 
7130CPR 
7140CPR 
7150CPR 
7160CPR 
7170CPR 
7180CPR 
7190CPR 
7840 


ROUTINE — TIMCK 

DATE/WRITTEN BY— 8/18/82 RJ MAFFEO 
DATE/REVISED BY — 

FUNCTION/PURPOSE — 

THIS ROUTINE CHECKS TO SEE IF THE PRESENT TIME VALUE 
IS ONE Or THE TIMES DESIRED 


CALLING SEQUENCE— CALL TIMCK(TIME , ITIME ) 
CALLING ARGUMENTS— 


NAME 

ATTRIBUTES 

DEFINITION 

TIME 

(I) 

PRESENT TIME VALUE OF TIME 

ITIME 

(0) 

DESIRED SWITCH 

ITIME=0 »> NOT DESIRED 
ITIME* 1 »> DESIRED 

FILES USED- 

NONE 


COMMONS USED- 

- 



CNTLTM — TIME/TEMPERATURE CONTROL PARAMATERS 


FUNCTIONS/ROUTINES CALLED-- NONE 

LIBRARIES ACCESSED — 

NONE 


LOCAL VARIABLES — 

SPECIAL COMMENTS ABOUT THIS ROUTINE 

SPECIAL REMARKS/ INSTRUCTIONS— 

HOST = H6000 

**«***********************K****** 

END 
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7860C SS * 8 S8SSSS3S&SS3S3SSSSSSSBSSSSSSSB 


sssssse 


7880 

79 10C*= 

7920CPR 

7930CPR 

7940CPR 

7950CPR 

7960CPR 

7970CPR 

7980CPR 

7990CPR 

8000CPR 

8010CPR 

8020CPR 

8030CPR 

8040CPR 

8050CPR 

8060CPR 

8070CPR 

8080CPR 

8090CPR 

8100CPR 

81 10CPR 

8120CPR 

8130CPR 

8140CPR 

8150CPR 

8160CPR 

8170CPR 

8180CPR 

8190CPR 

8200CPR 

8210CPR 

8220CPR 

8230CPR 

8240CPR 

8250CPR 

8260CPR 

8270CPR 

8280CPR 

8290CPR 

8300CPR 

8310CPR 

8320CPR 

8330CPR 

8340CPR 

8350CPR 

8360CPR 

8370CPR 

8380CPR 

8390CPR 

8400CPR 

8410CPR 

8420CPR 

8430CPR 

8440CPR 

8450CPR 

8460CPR 

8470CPR 

8480CPR 

9510 


SUBROUTINE CHNODT 


ROUTINE— CHNODT 

DATE/WRITTEN BY — 8/18/82 RJ MAFFEO 

DATE/REVISED BY — 

FUNCTION/PURPOSE — 

THIS ROUTINE MAKES TWO BASIC CHECKS ON THE INPUT 

IT CHECKS TO SEE THAT THE NUMBER OF ELEMENTS IN EACH TIMESTEP 

IS GREATER THAN OR EQUAL TO THE NUMBER OF GEOMETRY ELEMENTS 

IT ALSO CHECKS TO INSURE THAT ALL ELMS ON THE 

NEUTRAL ELM FILE HAVE TEMPERATURES ON THE THT 

OUTPUT FILE 


CALLING SEOUENCE— CALL CHNODT 

CALLING ARGUMENTS — 

NAME ATTRIBUTES DEFINITION 


FILES USED — 

I FLO — THT OUTPUT FILE 
I ERF — ERROR REMARK FILE 

COMMONS USED — 

AFIL 

WORK 

SIZE 

ELMDAT 

CNTLFL 


FUNCTIONS/ROUTINES CALLED — 

ROMANY 

RWELMS 

LIBRARIES ACCESSED — 

NONE 


LOCAL VARIABLES — 

SPECIAL COMMENTS ABOUT THIS ROUTINE 

SPECIAL REMARKS/INSTRUCTIONS— 

HOST « H6000 

******************************* 

END 
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9530C=== 

9550 

9580C=== 

9590CPR 

9600CPR 

9610CPR 

9620CPR 

9630CPR 

9640CPR 

965CCPR 

9660CPR 

9670CPR 

9680CPR 

9690CPR 

9700CPR 

9710CPR 

9720CPR 

9730CPR 

9740CPR 

9750CPR 

9760CPR 

9770CPR 

9780CPR 

9790CPR 

9800CPR 

9810CPR 

9820CPR 

9830CPR 

9840CPR 

9850CPR 

9860CPR 

9870CPR 

9880CPR 

9890CPR 

9900CPR 

9910CPR 

9920CPR 

9930CPR 

9940CPR 

9950CPR 

9960CPR 

9970CPR 

9980CPR 

9990CPR 

10000CPR 

10010CPR 

10020CPR 

10030CPR 

10040CPR 

10050CPR 

10060CPR 

10070CPR 

10080CPR 

10090CPR 

10100CPR 

10110CPR 

10120CPR 

10130CPR 

10750 


SUBROUTINE PREPGM 


ROUTINE-- PREPGM 

DATE/WRITTEN BY— 8/18/82 RJ MAFFEO 
DATE/REVISED BY — 

FUNCTION/PURPOSE — 

THIS ROUTINE CALLS ALL ROUTINES NEEDED TO 
TRANSFORM, SECTION (WINDOW) AND/OR 
ELIMINATE UNWANTED ELMS FROM THE PROBLEM 


CALLING SEQUENCE— CALL PREPGM 

CALLING ARGUMENTS-- NONE 

NAME ATTRIBUTES DEFINITION 

NONE NONE NONE 

FILES USED-- NONE 

COMMONS USED- 

SIZE 

ELMDAT 

CNTLGM 


FUNCTIONS/ROUTINES CALLED — 

RWELMS 

CTRAN 

SECTION 

REMELM 


LIBRARIES ACCESSED — 
NONE 


LOCAL VARIABLES — 

JSLC -- WINDOW SWITCH 

JSLC* 1 >>> ELM IS IN WINDOW 

USLC»-1 »> ELM 1$ NOT IN WINDOW 

SPECIAL COMMENTS ABOUT THIS ROUTINE 


SPECIAL REMARKS/ INSTRUCT IONS — 
HOST * H6000 

END 


GE/TRANCITS 


PREPGM 


PAGE* 8 


85 



TITLE' DATE =07/15/85 TIME = 23 42 PAGE' 9 

1 0770C 8S88S8S8S88SS8SS33S8eSSS:SSSS83SSBSSB888SSa3S8S8SS88E=S888SSa 

10790 SUBROUTINE CTRAN( XO , YO , ZO , C ) 

10830CPR 

10840CPR 

10850CPR ROUTINE-- CTRAN 

10860CPR 

10870CPR DATE/WRITTEN BY-- 8/18/82 RJ MAFFEO 

10880CPR DATE/REVISED BY-- 

10890CPR 

10900CPR FUNCTION/PURPOSE — 

10910CPR 

10920CPR THIS ROUTINE TRANSFORMS THE COORDINATES OF THE EIGHT 

10930CPR CORNER POINTS OF EACH ELM 

10940CPR TRANSFORMS ARE LIMITED TO TRANSLATIONS 

10950CPR 

10960CPR 

10970CPR CALLING SEQUENCE-- CALL CTRAN ( XO , YO , ZO , C ) 

10980CPR 

10990CPR CALLING ARGUMENTS-- 


1 1000CPR 

NAME 

ATTRIBUTES 

DEFINITION 

1 1010CPR 

XO 

(I) 

X OFFSET 

1 1020CPR 

YO 

(I) 

Y OFFSET 

1 1030CPR 

ZO 

(I) 

Z OFFSET 

1 1040CPR 

C 

(I/O) ARRAY 

COORDINATE ARRAY 


1 1050CPR 
1 1060CPR 
1 1070CPR 

1 1080CPR FILES USED-- NONE 

1 1090CPR 

1 1 100CPR COMMONS USED-- NONE 

1 1 1 10CPR 
1 1 120CPR 

11130CPR FUNCTIONS/ROUTINES CALLED-- NONE 

1 1 140CPR 

1 1 150CPR LIBRARIES ACCESSED-- 

1 1 160CPR NONE 

1 1 170CPR 
1 1 180CPR 

1 1 190CPR LOCAL VARI ABLES-- 

1 1200CPR 

11210CPR SPECIAL COMMENTS ABOUT THIS ROUTINE 

1 1220CPR 
1 1230CPR 

1 1240CPR SPECIAL REMARKS/ INSTRUCTIONS-- 

1 1250CPR HOST * H6000 

1 1260CPR 

1 1 270CPR ********************************* 

11440 END 
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1 1460C*********************************************************** 

11480 SUBROUTINE SCTION(C, XMN, YMN, ZMN, XMX , YMX , 2MX , INDT ) 

1151 OC >>5SSS>3SSisss&aaesmscs8::srsE:esecE>£zE3S3ss: = siEE:;EisisE! 

1 1520CPR 
1 1530CPR 

1 1540CPR ROUTINE-- 5CTI0N 

1 1550CPR 

1 1 560CPR DATE/WRITTEN BY-- 8/18/82 RJ MAFFEO 

1 1570CPR DATE/REVISED BY — 

1 1 580CPR 

1 1590CPR FUNCTION/PURPOSE — 

1 1600CPR 

11610CPR THIS ROUTINE CHECKS TO SEE IF ANY OF THE CORNER 

1 1620CPR COORDINATES ARE OUTSIDE OF THE WINDOW SPECIFIED 

1 1630CPR BY THE MIN AND MAX VALUES OF X,Y,Z 
1 1640CPR NOTE IF ANY POINT OF AN ELM IS OUTSIDE THE WINDOW 
11650CPR THE ELM IS CONSIDERED OUT OF THE WINDOW. 

1 1660CPR 
1 1670CPR 
1 1680CPR 

1 1690CPR CALLING ARGUMENTS — 


1 1700CPR 

NAME 

ATTRIBUTES 

DEFINITION 



11710CPR 
1 1720CPR 

C 

(I) 

COORDINATE ARRAY 



1 1730CPR 

XMN > 





1 1740CPR 

YMN » 

(I) 

MIN VALUE OF X,Y,Z 

FOR 

WINDOW 

1 1750CPR 

ZMN > 





1 1760CPR 

XMX > 





1 1770CPR 

YMX » 

(I) 

MAX VALUE OF X.Y.Z 

FOR 

WINDOW 

1 1 780CPR 
1 1790CPR 

ZMX > 





1 1800CPR 

INOT 

(0) 

IN/OUT SWITCH 




11810CPR I NOT * 1 » ELM IS IN THE WINDOW 

1 1820CPR I NOT*- 1 >> ELM IS OUT OF WINDOW 

1 1830CPR 
1 1840CPR 
1 1850CPR 

11860CPR FILES USED-- NONE 

1 1870CPR 

1 1880CPR COMMONS USED— NONE 
1 1890CPR 
1 1900CPR 

11910CPR FUNCTIONS/ROUTINES CALLED-- NONE 

1 1920CPR 

1 1930CPR LIBRARIES ACCESSED— 

1 1940CPR NONE 

1 1950CPR 
1 1960CPR 

1 1970CPR LOCAL VARIABLES— 

1 1980CPR 

1 1990CPR SPECIAL COMMENTS ABOUT THIS ROUTINE 

12000CPR 
1 2010CPR 

12020CPR SPECIAL REMARKS/INSTRUCTIONS— 

1 2030CPR HOST ■ H6000 

12040CPR 

12050CPR ********************************** 

12360 END 


GE/TRANCITS 


SCTION 


PAGE* 10 


87 


TITLE* 


DATE *07/15/85 


TIME * 23 42 


PAGE* 1 1 


12380C******** — 


=S*SSS3=ESC=SS3e==3=SSC8=SS=S=3= 


12400 SUBROUTINE REMELM 

1 24 20C 

12440CPR 


1 2450CPR* SSSBS8SSXSSSS33SSSSSSZSSBSSSSSS&5SS&S3BS33SSSBSBSSXS 


12460CPR 

12470CPR 

12480CPR 

12490CPR 

12500CPR 

12510CPR 

12520CPR 

12530CPR 

12540CPR 

12550CPR 

12560CPR 

1 2570CPR 

12580CPR 

12590CPR 

12600CPR 

12610CPR 

12620CPR 

12630CPR 

12640CPR 

12650CPR 

1 2660CPR 

12670CPR 

12680CPR 

1 2690CPR 

1 2700CPR 

12710CPR 

1 2720CPR 

12730CPR 

1 2740CPR 

12750CPR 

12760CPR 

1 2770CPR 

1 2780CPR 

12790CPR 

1 2800CPR 

128 10CPR 

12820CPR 

12830CPR 

12840CPR 

12850CPR 

12860CPR 

12870CPR 

12880CPR 

12890CPR* 

13470 


ROUTINE— REMELM 

DATE/WRITTEN BY— 09/13/82 RJ MAFFEQ 

DATE/REVISED BY- 

FUNCTION/ PURPOSE — 

THIS SUB ELIMINATES UNDESIRED ELMS FROM THE RANDOM 
ELM GEOM FILE 


CALLING ARGUMENTS — 

NAME ATTRIBUTES DEFINITION 


COMMONS USED — 

ELMDAT 

CNTLGM 

FUNCTIONS/ROUTINES CALLED 
RWELMS 

FILES USED -- 


LIBRARIES ACCESSED — 
NONE 


LOCAL VARIABLES — 

SPECIAL COMMENTS ABOUT THIS ROUTINE 

SPECIAL REMARKS/INSTRUCTIONS — 

HOST * H6000 


«*i3is*i««tsEC#Ee3r:a2nnaii*rsit:cirEi2:»s:ta3Ecr 

END 
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13490C===== 

13510 

135300===** 

13540CPR 

13550CPR 

13560CPR 

13570CPR 

1 3580CPR 

13590CPR 

13600CPR 

13610CPR 

13620CPR 

13630CPR 

13640CPR 

1 3650CPR 

13660CPR 

13670CPR 

13680CPR 

13690CPR 

13700CPR 

13710CPR 

13720CPR 

13730CPR 

13740CPR 

13750CPR 

13760CPR 

13770CPR 

13780CPR 

13790CPR 

13800CPR 

13810CPR 

13820CPR 

13830CPR 

13840CPR 

1 3850CPR 

13860CPR 

13870CPR 

13880CPR 

13890CPR 

13900CPR 

13910CPR 

13920CPR 

1 3930CPR 

13940CPR 

13950CPR 

13960CPR 

13970CPR 

13980CPR 

13990CPR 

14000CPR 

14010CPR 

14020CPR 

14030CPR 

14040CPR 

14050CPR 

14060CPR 

14070CPR 

14080CPR 

14090CPR 

14100CPR 

14110CPR 

14120CPR 

14 130CPR 

15130 


■XSS=SXS=XSCSS=2C 


SUBROUTINE CALCRD 


ROUTINE-- CALCRD 

DATE/WRITTEN BY-- 8/08/82 RJ MAFFEO 
DATE/REVISED BY — 

FUNCTION/PURPOSE-- 

THE PURPOSE OF THIS ROUTINE IS TO CALL ALL 
REQUIRED ROUTINES TO COMPUTE THE COORDINATES OF THE 
CENTROID OF THE ELMS AND THE COORDINATES OF THE 
FACE CENTERS IT THEN CALLS THE ELM I/O ROUTINE TO STORE 
THIS INFO INTO THE INTERNAL DATA BASE 


CALLING SEQUENCE— CALL CALCRD 

CALLING ARGUMENTS-- NONE 

NAME ATTRIBUTES DEFINITION 

NONE NONE NONE 

FILES USED— 

IERF — ERROR REMARK FILE 

COMMONS USED— 

AFIL -- FILE CODE COMMON 

ELMDAT -- ELM ATTRIBUTE COMMON 

SIZE — PROBLEM SIZE COMMON 

WORK — SCRATCH WORK SPACE COMMON 

CRNDEF -- LOCAL COMMON TO COORD CALC ROUTINES 

FUNCTIONS/ROUTINES CALLED— 

GETPD 

RWELMS 

VOLM 

AREAD 

REACRD 

LIBRARIES ACCESSED— 

NONE 

LOCAL VARIABLES — 

SPECIAL COMMENTS ABOUT THIS ROUTINE 

MANY OF THE ROUTINES CALLED BY CALCRD ARE ALSO USED 
IN SIGMA THTD DECK GENERATOR THEY CAN COMPUTE MUCH MORE THAN 
JUST COORD I NATES ( SUCH AS VOLUMES AREAS AND DELTA X). 


SPECIAL REMARKS/INSTRUCTIONS— 
HOST * H6000 

******************************** 

END 
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151 50C SSSSS=5SSSSSKSSeS35SSSSSSSSe=SSSSSSXSSSSSSSSSSSS=3Z= 

15170 SUBROUTINE GETPD( P3D , P2D ) 

151 90C =SSSSae==BSSasaS=B==SSSSS3SC=SSS3aSS3SS3SB3SSSSSSS3=== 


15210CPR 

15220CPR== 

15230CPR 

15240CPR 

15250CPR 

15260CPR 

15270CPR 

1 5280CPR 

15290CPR 

15300CPR 

15310CPR 

15320CPR 

15330CPR 

15340CPR 

15350CPR 

1 5360CPR 

15370CPR 

15380CPR 

15390CPR 

15400CPR 

15410CPR 

15420CPR 

1 5430CPR 

15440CPR 

15450CPR 

15460CPR 

1 5470CPR 

1 5480CPR 

1 5490CPR 

1 5500CPR 

15510CPR 

15520CPR 

1 5530CPR 

1 5540CPR 

15550CPR 

15560CPR 

1 5570CPR 

15580CPR 

15590CPR 

15600CPR 

15610CPR 

15620CPR 

15630CPR 

15640CPR 

15650CPR 

15660CPR 

15670CPR*=i 

16020 


S33SSS=33333SSSSSS3333S3SSSS£333333SaS=3=B33SSS333S 


ROUTINE— GETPD 

DATE/WRITTEN BY-- 09/13/82 RJ MAFFEO 

DATE/REVISED BY- 

FUNCTION/ PURPOSE — 

THIS SUB COMPUTES AND STORES ALL OF THE PARTIAL DERIVATIVES 
OF THE SHAPE FUNCTIONS . BOTH 3D AND 2D 


CALLING ARGUMENTS — 

NAME ATTRIBUTES DEFINITION 

P3D (0) 3D PARTIAL DERIVATIVES 

P2D (0) 2D PARTIAL DERIVATIVES 

COMMONS USED — 

CRNDEF — SEE CALCRD 

FUNCTIONS/ROUTINES CALLED — 

PARD 

GTGAUS 

FILES USED — 


LIBRARIES ACCESSED — 
NONE 


LOCAL VARIABLES — 

SPECIAL COMMENTS ABOUT THIS ROUTINE 

SPECIAL REMARKS/INSTRUCTIONS — 

HOST = H6000 


END 
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16060 SUBROUTINE VOLM(X , Y , Z , GP , P3D , VOL , XC , YC , ZC , Q , I VER ) 

1 6080C SSS *SSSSXaa:EtSSK8 = SSKSCS = EXSSSS83KSS = SS3=CSZS = SSSS = SSSSSS3 = Z = = 

16 100CPR 

161 1 OCPR ss =sx e v KS s KS ss KK as esksss s sesx3sassss=8s===se assesses* 

16 120CPR 
16130CPR 

16140CPR ROUTINE— VOLM 

16 150CPR 

16 160CPR DATE/WRITTEN BY— 09/13/82 RJ MAFFEO 

16170CPR DATE/REVISED BY — 

16 180CPR 

16 190CPR FUNCTION/PURPOSE— 

16200CPR 

16210CPR THIS SUB DOES NUMERICAL INTEGRATION FOR VOL.XC.YC.ZC 

16220CPR IT COMPUTES DETJ AND THE JACOBIAN 

16230CPR 

16240CPR 

16250CPR 

16260CPR CALLING ARGUMENTS — 


16270CPR 

NAME 

ATTRIBUTES 

DEFINITION 

16280CPR 

16290CPR 

X 

(I) 

ARRAY OF X COORDS 

16300CPR 

Y 

(X) 

/ 4 Y 

16310CPR 

Z 

(I) 

' ' Z 

16320CPR 

GP 

(X) 

ARRAY OF GAUSS POINTS 

16330CPR 

P3D 

(I) 

3D PARTIAL DERIVATIVES 

16340CPR 

VOL 

(0) 

VOLUME 

16350CPR 

XC 

(0) 

X COORD OF CENTROID 

16360CPR 

YC 

(0) 

Y 

16370CPR 

ZC 

(0) 

Z 

16380CPR 

0 

(0) 

JACOBIAN 

16390CPR 

IVER 

(0) 

IVER* 1 > ZERO VOLUME ELM 

16400CPR 

16410CPR 

COMMONS USED 

_ _ 


16420CPR 

16430CPR 

16440CPR FUNCTIONS/ROUTINES CALLED — 
16450CPR 

16460CPR SHAPEF 

16470CPR 

16480CPR FILES USED — 

16490CPR 

16500CPR 

16510CPR LIBRARIES ACCESSED — 

16520CPR NONE 

16530CPR 

16540CPR 

16550CPR LOCAL VARIABLES — 

16560CPR 

16570CPR SPECIAL COMMENTS ABOUT THIS ROUTINE 

16580CPR 

16590CPR 

16600CPR SPECIAL REMARKS/INSTRUCTIONS — 


16610CPR HOST * H6000 

16620CPR 

16630CPR 

1 6640CPR® * ***** B ** 11 ************»*««Bre*»«BEC***S*E*3S*lStEB 

17150 END 
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1 7 170C* atBSSSSBSSSSaB * ss = s = csssBss = 33 = sssz = ssse = z3Z3Ssaa = es = ass = E2 
17190 SUBROUTINE SHAPEF( X , V , 2 , TETA f YNETA t ZETA , XP f YP , ZP ) 

172 10C s = 3S = :2 = == = = - 5 *= SSSB = = = = stnnnl: > : = = 3 3*3i = = ss33 = 3 ===;;s:r;: = : 

17230CPR 

17240CPP sssssssss * 3 * sssssssxssssie83E 83BSBSstesssessssssEE5sss 

17250CPR 

17260CPR 

17270CPR ROUTINE — SHAPEF 

17280CPR 

17290CPR DATE/WRITTEN BY-- 09/14/82 RJ MAFFEO 

17300CPR DATE/REVISED BY — 

17310CPR 

17320CPR FUNCTION/PURPQSE — 

17330CPR 

17340CPR THIS SUB COMPUTES THE VALUES OF X,Y,Z AT THE GUASS POINTS 

17350CPR BASED ON LINEAR 3D SHAPE FUNCTIONS 

17360CPR 

17370CPR 

17380CPR 

17390CPR 

17400CPR CALLING ARGUMENTS — 

17410CPR NAME ATTRIBUTES DEFINITION 

17420CPR 

17430CPR X (I) > 

17440CPR Y (I) » ARRAYS OF COORDINATES 

17450CPR Z (I) > 

17460CPR TETA (I) > 

17470CPR YNETA (I) » LOCAL COORDINATES 

1 7480CPR ZETA (I) > 

1 7490CPR XP (0) > 

17500CPR YP (0) » DESIRED GLOBAL COORDINATES 

17510CPR ZP (0) > 

1 7520CPR 

1 7530CPR COMMONS USED -- 
1 7540CPR 

17550CPR CRNDEF 

17560CPR 

17570CPR FUNCTIONS/ROUTINES CALLED -- 

1 7580CPR 

17590CPR 

17600CPR FILES USED — 

17610CPR 

17620CPR 

17630CPR LIBRARIES ACCESSED — 

17640CPR NONE 

17650CPR 

17660CPR 

17670CPR LOCAL VARIABLES — 

17680CPR 

17690CPR SPECIAL COMMENTS ABOUT THIS ROUTINE 

17700CPR 

17710CPR 

17720CPR SPECIAL REMARKS/INSTRUCTIONS — 

17730CPR HOST * H6000 

17740CPR 

17750CPR 

17760CPR“* r ** 

17950 END 
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1 7970C s ** B = B **** = »****S33*****a* = s:s = a = **«a = a* = ac3*a* = sa**ac* = 

17990 SUBROUTINE PARD(TETA f YNETA t ZETA , IGP , IND. PD) 

1801 OC 3E = = = = = *:i = csccs!z = 3 = = = = = c = t = ctE = = = s = = = zEa = = c = = = = zai= = 3 = 

18030CPR 

1 8040CPR = s = = = = * = “= : -** = = = = = as*ss = = = = ss = s = = = x = s = = *c = s = = = = s3mne = 

1 8050CPR 
1 8060CPR 

1 8070CPR ROUTINE-- PARD 

1 8080CPR 

18090CPR DATE/WRITTEN BY— 09/14/82 RJ MAFFEQ 

18 100CPR DATE/REVISED BY — 

18110CPR 

18 120CPR FUNCTION/PURPOSE — 

18130CPR 

18140CPR THIS SUB COMPUTES THE DERIVATIVES OF THE 
18150CPR SHAPE FUNCTIONS WRT THE LOCAL SYSTEM 
18160CPR BOTH FOR V0LUMES(3D) AND AREAS ( 2D ) 

18170CPR 
18 180CPR 
1 8 190CPR 

1 8200CPR CALLING ARGUMENTS — 


18210CPR 

NAME 

ATTRIBUTES 

DEFINITION 

18220CPR 




18230CPR 

TETA 

(I) > 


18240CPR 

YNETA 

(I) » 

LOCAL COORDINATES 

1 8250CPR 

ZETA 

(I) > 


18260CPR 

IGP 

(I) 

GAUSS POINT COUNTER 

18270CPR 

IND 

(I) 

2D/3D INDICATOR 

18280CPR 



IND*0 » 3D 

18290CPR 



I ND* 1 THRU 6 » 2D IND*FACE NUMBER 

1 8300CPR 

PD 

(0) 

ARRAY OF DERIVATIVES 

18310CPR 




18320CPR 

COMMONS USED 

-- 



18330CPR 

1 8340CPR CRNDEF 

1 8350CPR 

18360CPR FUNCTIONS/ROUTINES CALLED — 

1 8370CPR 
1 8380CPR 

1 8390CPR FILES USED -- 

18400CPR 

18410CPR 

18420CPR LIBRARIES ACCESSED -- 
18430CPR NONE 
18440CPR 
1 8450CPR 

18460CPR LOCAL VARIABLES — 

18470CPR 

18480CPR SPECIAL COMMENTS ABOUT THIS ROUTINE 

18490CPR 

18500CPR 

18510CPR SPECIAL REMARKS/INSTRUCTIONS — 

1 8520CPR HOST = H6000 

1 8530CPR 

18540CPR 

1 8550CPR* ** = a*= = e: = :itit = is: = z = iEi:Bris*XEt!B£rEESBc = = = a:i*3E£3 
18950 END 
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1 8970C=== SS3BSS8sa3SESSS8SSSSSStSSai = rSSSSS3t3SSSSSSSSaSSS 

18990 SUBROUTINE AREAD( X , Y , 2 , P2D , 0 . ARA , DX2 , 

1 9020C sssBasaaaaassasaaxsasassBssaass&aassaasssEzaasssasaa 

19040CPR 

1 9050CPR* = 3 = = ^= = ssssss:sss = = sssssssssrsssx3:ssssssssssssss:sssss 

19060CPR 

19070CPR 

19080CPR ROUTINE-- AREAD 

19090CPR 

19 100CPR DATE/WRITTEN BY-- 09/14/82 RJ MAFFEO 

19110CPR DATE/REVISED BY — 

19120CPR 

19130CPR FUNCTION/PURPOSE-- 
19 140CPR 

1 9150CPR THIS SUB COMPUTES THE FACE AREAS AND DELTA X 
19 160CPR FOR 3D THT ELMS 

19 170CPR NOTE IT COMPUTES DELTA X TWO WAYS 

19 180CPR IDX* 1 DX=DIST FROM ELM CENTROID TO FACE CENTROID 

19190CPR IDX=2 . DX=PERPENDICULAR DIST FROM ELM CENT TO FACE 
19200CPR IDX*3 SWITCH TO JUST CALC CENTROID OF FACE 5 

19210CPR 

19220CPR IT IS BEING USED HERE TO COMPUTE THE COORDINATES 

19230CPR OF THE FACE CENTERS 

19240CPR 

19250CPR 

19260CPR 

19270CPR 


192BOCPR 

CALLING 

ARGUMENTS 

— 


19290CPR 
1 9300CPR 

NAME 

ATTRIBUTES 

DEFINITION 

19310CPR 

X 

a) 

> 


19320CPR 

> 

(i) 

>> 

ARRAY OF CORNER COORDS 

19330CPR 

Z 

(i) 

> 


1 9340CPR 

P2D 

(i) 


DERIVATIVES OF SHAPE FUNCTIONS 

19350CPR 

0 

(0) 


JACOBIAN 

19360CPR 

ARA 

(0) 


ARRAY OF FACE AREAS 

19370CPR 

DX2 

(0) 


ARRAY OF FACE DELTA XS 

19380CPR 

IDX 

(I) 


DELTA X COMPUTATION INDICATOR 

19390CPR 

XC 

(I) 

> 


1 9400CPR 

YC 

(I) 

>> 

COORDS OF VOLUME CENTROID 

19410CPR 

ZC 

(I) 

> 


19420CPR 

XA 

(0) 

> 


19430CPR 

YA 

(0) 

>> 

ARRAYS OF COORDS FOR FACE CENTERS 

19440CPR 

ZA 

(0) 

> 



19450CPR 

19460CPR COMMONS USED -- 
19470CPR 

19480CPR CRNDEF 

19490CPR 

19500CPR FUNCTIONS/ROUTINES CALLED — 
19510CPR 

19520CPR GTGAUS 

19530CPR SHAPEF 

19540CPR 

19550CPR FILES USED -- 
19560CPR 
1 9570CPR 

19580CPR LIBRARIES ACCESSED — 

19590CPR NONE 

19600CPR 

19610CPR 

19620CPR LOCAL VARIABLES — 

19630CPR 

19640CPR SPECIAL COMMENTS ABOUT THIS ROUTINE 

19650CPR 

19660CPR 

19670CPR SPECIAL REMARKS/INSTRUCTIONS — 
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19680CPR HOST 

19690CPR 

19700CPR 

197 10CPR==****===== 
20780 END 


H6000 

sscesssasssssssssssssssssssssssssssssrss 
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208000==== SSSSSSSSSSSSSIBSaBSSSSBIBSSBBBSBBSXSBSSlSSS 

20820 SUBROUTINE GTGAUS( IND, 1 1 . 12 , 13 , IP ) 

20840C 88888888888 ***B**=**3*888*8SBB8B88BBBBBBBB8BBBS 

20860CPR 

20870CPR=*= S8 =*SSBSSSSSSBI8SSSX88BBBSSXSX8B8SS3SX8SBBSSS3SXBB 

20880CPR 

20890CPR 

20900CPR ROUTINE— GTGAUS 

20910CPR 

20920CPR DATE/WRITTEN BY— 09/14/82 RU MAFFEO 

20930CPR DATE/REVISED BY — 

20940CPR 

20950CPR FUNCTION/PURPOSE — 

20960CPR 

20970CPR THIS SUB GETS APPROPRIATE GAUSS POINTS 

20980CPR FOR EACH FACE OF THE 3D BRICK 

20990CPR 

21000CPR 

21010CPR 

21020CPR 

21030CPR CALLING ARGUMENTS — 


21040CPR 

21050CPR 

NAME 

ATTRIBUTES 

DEFINITION 

21060CPR 

IND 

(I) 

FACE INDICATOR 

21070CPR 

11 

(0) > 


2 1080CPR 

12 

(0) > 


2 1090CPR 

13 

(0) » 

GAUSS POINT COUNTERS 

21 lOOCPR 

IP 

(0) > 



211 10CPR 

21 120CPR COMMONS USED — 

21130CPR 

21140CPR 

21 150CPR FUNCTIONS/ROUTINES CALLED — 

21 160CPR 
21 170CPR 

21180CPR FILES USED — 

21 190CPR 
21200CPR 

21210CPR LIBRARIES ACCESSED -- 

2 1220CPR NONE 

21230CPR 

21240CPR 

21250CPR LOCAL VARIABLES — 

21260CPR 

21270CPR SPECIAL COMMENTS ABOUT THIS ROUTINE 

2 1280CPR 

21290CPR 

21300CPR SPECIAL REMARKS/INSTRUCTIONS -- 

21310CPR HOST = H6000 

21320CPR 

21330CPR 

2 1 340CPR BBBSSBBBSBBBBBBSBBBBBBBBBaBBBSSESBEBBBBBBBBESSBBrBBr 

21450 END 
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2147 0C = 

21490 SUBROUTINE REACRD(CRD,XC,YC,ZC,XA,YA,2A) 

2151 OC = = = t: = = a = = 3c = = :s = = = = = rx = = = C 3-E = = - = t = z = = ri 3 :; - EE!: __ = = _ = = 

2 1 530CPR 

2 1 540CPR* * = c:r = = 3:a: = :: = EC£::zaioi = r:22s = ;crisa3 = :2tisrr = :Ba; 

2 1 550CPR 
21560CPR 

21570CPR ROUTINE-- REACRD 

2 1 580CPR 

2 1 590CPR DATE/WRITTEN BY-- 09/15/82 Rd MAFFEO 

21600CPR DATE/REVISED BY-- 

21610CPR 

2 1620CPR FUNCT I ON/PURPOSE -- 

2 1630CPR 

2 1640CPR THIS SUB REARRANGES THE FACE CENTER COORDINATE ARRAYS 

2 1650CPR AND THE VOLUME CENTROIDS INTO THE CRD ARRAY 

2 1660CPR 
2 1670CPR 
2 1680CPR 
2 1 690CPR 


2 1 700CPR 

CALLING 

ARGUMENTS 

-- 


21710CPR 

NAME 

ATTRIBUTES 

DEFINITION 

2 1 720CPR 





2 1 730CPR 

CRD 

(0) 


BRICK COORD ARRAY 

2 1740CPR 

XC 

(I) 

> 


2 1 750CPR 

YC 

(I) 

>> 

VOLUME CENTROIDS 

21760CPR 

ZC 

(I) 

> 


21770CPR 

XA 

(I) 

> 


21780CPR 

YA 

(I) 

>> 

FACE CENTER ARRAYS 

2 1 790CPR 

ZA 

(I) 

> 


2 1800CPR 





21810CPR 

COMMONS 

USED -- 




2 1820CPR 
21830CPR 

2 1840CPR FUNCTIONS/ROUTINES CALLED -- 

2 1850CPR 

21860CPR 

2 1870CPR FILES USED -- 

2 1880CPR 

21890CPR 

2 1900CPR LIBRARIES ACCESSED -- 
21910CPR NONE 
2 1 920CPR 
2 1930CPR 

2 1940CPR LOCAL VARIABLES -- 

2 1950CPR 

2 1960CPR SPECIAL COMMENTS ABOUT THIS ROUTINE 

2 1970CPR 
2 1980CPR 

2 1990CPR SPECIAL REMARKS/INSTRUCTIONS — 

22000CPR HOST * H6000 

22010CPR 

22020CPR 

22030CPR*** ******* ************=«=»*B=«**===**s=s=a***««c=*== 
22240 END 
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22260C== 

22280 

223000 = = 

22320CPR 

22330CPR* 

22340CPR 

22350CPR 

22360CPR 

22370CPR 

22380CPR 

22390CPR 

22400CPR 

22410CPR 

22420CPR 

22430CPR 

22440CPR 

22450CPR 

22460CPR 

22470CPR 

22480CPR 

22490CPR 

22500CPR 

22510CPR 

22520CPR 

22530CPR 

22540CPR 

22550CPR 

22560CPR 

22570CPR 

22580CPR 

22590CPR 

22600CPR 

22610CPR 

22620CPR 

22630CPR 

22640CPR 

22650CPR 

22660CPR 

22670CPR 

22680CPR 

22690CPR 

22700CPR 

22710CPR 

22720CPR 

22730CPR 

22740CPR 

22750CPR 

22760CPR 

22770CPR 

22780CPR 

22790CPR 

22800CPR 

22810CPR 

22820CPR 

22830CPR 

22840CPR 

22850CPR 

22860CPR 

22870CPR 

22880CPR* 

24890 


SUBROUTINE CORNCC 


ROUTINE-- CORNCC 

DATE/WRITTEN BY-- 11/10/82 RJ MAFFEO 

DATE/REVISED BY — 

FUNCTION/PURPOSE — 

THIS SUB CALLS ALL ROUTINES NEEDED TO COMPUTE 
AND STORE THE INFORMATION ASSOCIATED WITH THE 
CORNER NODE WEIGHTING COEFFICIENTS 


CALLING ARGUMENTS — 

NAME ATTRIBUTES DEFINITION 

NONE 

COMMONS USED — 

AFIL 

ELMDAT 

SIZE 

WORK 

CNTLGM 

FUNCTIONS/ROUTINES CALLED — 

RWELMS 

SORT 

CCNVEC 

GETCRD 

CINV 

PUTWFI 

RWCCF 

FILES USED -- 
I ERF S 

LIBRARIES ACCESSED -- 
NONE 


LOCAL VARIABLES — 

SPECIAL COMMENTS ABOUT THIS ROUTINE 

SPECIAL REMARKS/INSTRUCTIONS — 

HOST * H6000 


3SSBXSXBSXXS8SS8X8XS88&3&XSXSS8XXXX888SSX8XBXCSXB 

END 
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249 IOC******* — *** — ****************************** — * — ******** 

24930 SUBROUTINE CCNVEC ( NOD , IRW , NLOC . NUMA , MCON, MXNCON, NSIZW ) 

24950C* * *****33833333338338333383338333338333333888338833383 

24970CPR 

24990CPR 

25000CPR 

25010CPR ROUTINE— CCNVEC 

25020CPR 

25030CPR DATE/WRITTEN BY— 11/10/82 RU MAFFEO 

25040CPR DATE/REVISED BY — 

25050CPR 

25060CPR FUNCTION/PURPOSE-- 
25070CPR 

25080CPR THIS SUB WRITES TO (CREATES) AND READS FROM THE 
25090CPR NODE -CONNECTION VECTOR (MCON) 

25100CPR 

25110CPR 

25 120CPR CALLING ARGUMENTS — 


25130CPR 
25 140CPR 

NAME 

ATTRIBUTES 

DEFINITION 

25150CPR 

NOD 

(I) 

NODE NUMBER 

25 160CPR 

IRW 

(I) 

READ/WRITE SWITCH(0-READ . 1-WRITE) 

25170CPR 

NLOC 

(I/O) 

CURRENT CONNECTION TO NODE 

25 180CPR 

NUMA 

(0) 

MAX 0 CONNECTIONS TO NODE 

25190CPR 

MCON 

(I/O) 

PACKED CONNECTION ARRAY 

25200CPR 

MXNCON 

(0) 

MAX 0 OF CONNECTIONS TO ANY NODE 

25210CPR 

25220CPR 

NSIZW 

(I) 

SIXE OF WORK ARRAY (MCON HERE) 

25230CPR 

25240CPR 

COMMONS USED 

AFIL 



25250CPR 

25260CPR 

25270CPR FUNCTIONS/ROUTINES CALLED — 

25280CPR 

25290CPR 

25300CPR FILES USED — 

25310CPR 

25320CPR IERF S 

25330CPR 

25340CPR LIBRARIES ACCESSED — 

25350CPR NONE 

25360CPR 

25370CPR 

25380CPR LOCAL VARIABLES — 

25390CPR 

25400CPR SPECIAL COMMENTS ABOUT THIS ROUTINE 
25410CPR THE NODE-CONNECTION VECTOR (MCON) IS KEYED 
25420CPR TO THE NODE NAME AND IS COMPOSED OF THE MAX 

25430CPR CONNECTION TO THAT NODE AND THE CURRENT CONNECTION TO 

25440CPR THAT NODE IN THE FORM — 

25450CPR 

25460CPR MCON l NOD ) * 1000* MAX CONNECTION + CURRENT CONNECTION 

25470CPR 

25480CPR THE INITIAL CURRENT COUNT IS EQUAL TO THE MAX CONNECTION. 
25490CPR AS THE CONNECTIONS ARE PROCESSED THE CURRENT COUNT IS 
25500CPR REDUCED BY ONE 

25510CPR 

25520CPR SPECIAL REMARKS/INSTRUCTIONS — 

25530CPR HOST « H600O 

25540CPR 

25550CPR 

25560CPR®* ■«>>Ba*BaKBB8>8SCimCBISBB«8B*K«IISSIBCSB8S8ISSSSSX8 

26040 END 
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26060C sssssssssassssssssssessssBssssssssssssssssssssssss 

26080 SUBROUTINE GETCRD ( NE , CRD . BRD , XCN, YCN, ZCN, JREC ) 

26 1 OOC SSS = SSS8S88SSSSSSBSS88SSSSStB3t8SSS8SasSSS3S3XS3SSS 

26 120CPR 

26 130CPR ==ssssssssssss=ss3ssssssscssssssssssssssss * ssssss **s = 
26 140CPR 
26 150CPR 

26 160CPR ROUTINE— GETCRD 

26 170CPR 

26 180CPR DATE/WRITTEN BY— 11/11/82 RJ MAFFEO 

26 190CPR DATE/REVISED BY — 

26200CPR 

26210CPR FUNCTION/PURPOSE — 

26220CPR 

26230CPR THIS SUB EXTRACTS THE APPROPRIATE SUBSET OF 

26240CPR COORDINATES FOR EACH CORNER NODE FORM THE 

26250CPR FULL SET OF COORDS OF EACH ELM 

26260CPR 

26270CPR 

26280CPR 

26290CPR 

26300CPR 


26310CPR 

CALLING 

ARGUMENTS 

— 


26320CPR 

26330CPR 

NAME 

ATTRIBUTES 

DEFINITION 

26340CPR 

NE 

c) 


ELMEMENT NAME 

26350CPR 

CRD 

(i) 


FULL SET OF COORDS 

26360CPR 

BRD 

(0) 


SUBSET OF COORDS 

26370CPR 

XCN 

(0) 

>» 


26380CPR 

YCN 

(0) 

>>>>>>>>> 

COORDS OF CORNER 

26390CPR 

ZCN 

(0) 

>>> 


26400CPR 

26410CPR 

JREC 

(I) 


CORNER LABEL 


26420CPR COMMONS USED — 

26430CPR 

26440CPR AFIL 

26450CPR 

26460CPR FUNCTIONS/ROUTINES CALLED — 

26470CPR 

264QOCPR 

26490CPR FILES USED — 

26500CPR 

26510CPR I ERF S 

26520CPR 

26530CPR LIBRARIES ACCESSED -- 

26540CPR NONE 

26550CPR 

26560CPR 

26570CPR LOCAL VARIABLES — 

26580CPR 

26590CPR SPECIAL COMMENTS ABOUT THIS ROUTINE 
26600CPR 

26610CPR THE ARRARY IDX IS USED TO DEFINE THE RELATIONSHIP 
26620CPR BETWEEN THE CORNER POINT AND THE APPROPRIATE 3 FACES 
26630CPR FOR EXAMPLE CORNER POINT 1 IS CONNECTED TO FACES 
26640CPR 1.4,5 

26650CPR THEREFORE IDX( 1 . 1 )*1 . IDX( 1 ,2) *4 , IDX ( 1 , 3 ) =5 

26660CPR 

26670CPR 

26680CPR 

26690CPR SPECIAL REMARKS/INSTRUCTIONS — 

26700CPR HOST 3 H6000 

26710CPR 

26720CPR 

26730CPR® sie*8ssss388s:xE:::s8iBscG::3:ss3sx:e£3C(8a8S8ss8iE 
27340 END 
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27360C SI8B **ss8SB* 

27380 SUBROUTINE PUTWFI (NOD , NE , NUMA , JREC . XSF , NPOS , ACWF ) 

27400C B * = » = S « a »««« 

27420CPR 

27430CPR- 5=s2ss:c:i;rs:r8iss3aasiir Eee3:=33= - I i =c=IS -. Bt3 _. =B 

27440CPR 

27450CPR 

27460CPR ROUTINE— PUTWFI 

27470CPR 

27480CPR DATE/WRITTEN BY— 11/11/82 R«J MAFFEO 

27490CPR DATE/REVISED BY — 

27500CPR 

27510CPR FUNCTION/PURPOSE — 

27520CPR 

27530CPR THIS SUB POSITIONS THE WEIGHTING COEFF 

27540CPR INFO INTO CORRECT LOCATIONS OF THE ACWF ARRAY 

27550CPR 

27560CPR 

27570CPR 

27580CPR 

27590CPR 

27600CPR CALLING ARGUMENTS — 


27610CPR 

NAME 

ATTRIBUTES 

DEFINITION 

27620CPR 




27630CPR 

NOD 

(I) 

CORNER NODE NUMBER 

27640CPR 

NE 

(I) 

ELM NUMBER 

27650CPR 

NUMA 

(I) 

MAX NUMBER OF CONNECTIONS 

27660CPR 

JREC 

(I) 

PACKED CORNER - FACE 0 

27670CPR 

XSF 

(I) 

WEIGHTING COEFFS 

27680CPR 

NPOS 

(I) 

STARTING POSITION FOR INFO 

27690CPR 

ACWF 

(0) 

ARRAY CONTAINING INFO 

27700CPR 




27710CPR 

COMMONS USED 

— 


27720CPR 




27730CPR 

AFIL 




27740CPR 

277BOCPR FUNCTIONS/ROUTINES CALLED — 

27760CPR 

27770CPR 

27780CPR FILES USED — 

27790CPR 

27800CPR 

27810CPR LIBRARIES ACCESSED — 

27820CPR NONE 

27830CPR 

27840CPR 

27850CPR LOCAL VARIABLES — 

27860CPR 

27870CPR SPECIAL COMMENTS ABOUT THIS ROUTINE 

27880CPR 

27890CPR 

27900CPR SPECIAL REMARKS/INSTRUCTIONS — 

27910CPR HOST * H6000 

27920CPR 

27930CPR 



28260 END 
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282800=*= 

28300 

283200=== 

28340CPR 

28350CPR* 

28360CPR 

28370CPR 

28380CPR 

28390CPR 

28400CPR 

28410CPR 

28420CPR 

28430CPR 

28440CPR 

28450CPR 

28460CPR 

28470CPR 

28480CPR 

28490CPR 

28500CPR 

28510CPR 

28520CPR 

2S530CPR 

28540CPR 

28550CPR 

28560CPR 

28570CPR. 

28580CPR 

28590CPR 

28600CPR 

28610CPR 

28620CPR 

28630CPR 

28640CPR 

28650CPR 

28660CPR 

28670CPR 

28680CPR 

28690CPR 

28700CPR 

28710CPR 

28720CPR 

28730CPR 

28740CPR 

28750CPR 

28760CPR 

28770CPR 

28780CPR 

28790CPR 

28800CPR 

28810CPR 

28820CPR 

28830CPR 

28840CPR* 3 

29120 


SUBROUTINE RWCCF ( NOD . ACWF , NPOS . I RW ) 


ROUTINE-- RWCCF 

DATE/WRITTEN BY-- 11/11/82 RJ MAFFEO 

DATE/REVISED BY-- 

FUNCTION/PURPOSE-- 

THIS ROUTINE READS AND WRITES THE CORNER NODE 
COEFF INFO TO AND FROM THE DATA BASE 


CALLING ARGUMENTS — 


NAME 

ATTRIBUTES 

DEFINITION 

NOD 

(I/O) 

NODE NAME (RECORD 0 ) 

ACWF 

(I/O) 

ARRAY OF INFO 

NPOS 

(I) 

0 OF TERMS TO BE READ/WRITTEN 

IRW 

(I) 

READ/WRITE SWITCH(0-READ , 1-WRITE) 

COMMONS USED 

-- 


AFIL 




FUNCTIONS/ROUTINES CALLED -- 


FILES USED -- 
IRCCF R 

LIBRARIES ACCESSED -- 
NONE 


LOCAL VARIABLES -- 

SPECIAL COMMENTS ABOUT THIS ROUTINE 

SPECIAL REMARKS/INSTRUCTIONS — 

HOST * H6000 


= = * = — sss = ssss3xs3se33 = ccess = sscss = = s = ssxsa = szs = = s 

END 
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29140C*** 8 

29160 SUBROUTINE GETMPE (MP , NWDMPE ) 



29200CPR 

29210CPR K3 SSBSSBSesSSS8&&te8SSSSS88«tSrS3B8S8SSSSSCB&teX8SSSS 

29220CPR 

29230CPR 

29240CPR ROUTINE-- GETMPE 

29250CPR 

29260CPR DATE/WRITTEN BY— 09/19/83 RJ MAFFEO 

29270CPR DATE/REVISED BY — 

29280CPR 

29290CPR FUNCTION/ PURPOSE — 

29300CPR 

29310CPR THIS ROUTINE READS THE FIRST TIME STEP TEMPERATURE 

29320CPR SOLUTION AND EXTRACTS THE ELM NAME AND RELATIVE POSITION 

29330CPR OF THIS ELM IN THE OUTPUT LIST. IT THEN CALLS MPERW TO 

29340CPR STORE THIS INFO INTO THE MP ARRAY 

29350CPR 

29360CPR 

29370CPR 

29380CPR 

29390CPR 

29400CPR CALLING ARGUMENTS — 

29410CPR NAME ATTRIBUTES DEFINITION 

29420CPR 

29430CPR MP (I/O) ELM MATRIX POSITION ARRAY 

29440CPR NWDMPE (I) tf OF WORDS IN MP ARRAY 

29450CPR 

29460CPR COMMONS USED -- 
29470CPR 

29480CPR AFIL 

29490CPR SIZE 

29500CPR CNTLFL 

29510CPR 

29520CPR FUNCTIONS/ROUTINES CALLED — 

29530CPR 

29540CPR TFLSET 

29550CPR HTFIO 

29560CPR MPERW 

29570CPR 

29580CPR FILES USED -- 

29590CPR 

29600CPR 

29610CPR LIBRARIES ACCESSED — 

29620CPR NONE 

29630CPR 

29640CPR 

29650CPR LOCAL VARIABLES — 

29660CPR 

29670CPR SPECIAL COMMENTS ABOUT THIS ROUTINE 

29680CPR 

29690CPR 

29700CPR SPECIAL REMARKS/INSTRUCTIONS — 

29710CPR HOST « H6000 * 

29720CPR 

29730CPR 

29740CPR-»« 

30240 END 
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30260C s:s 

30280 SUBROUTINE MPERW( MP , NE , NPOS , IRW) 

30300C == * ==3=I =3==*=*»===3s:3=3=3;crscz2n;ra=cx=;3==:r=22n!:r== 

30320CPR 

30330CPR s = = a = 3 = = = - = *»*** = = = = = * = = a = = 3* = = =* = s = = = ass = ea = = = = = a = = = s: 

30340CPR 

30350CPR 

30360CPR ROUTINE-- MPERW 

30370CPR 

30380CPR DATE/WRITTEN BY-- 10/08/83 RJ MAFFEO 

30390CPR DATE/REVISED BY-- 

30400CPR 

304 10CPR FUNCT I ON/PURPOSE -- 

30420CPR 

30430CPR THIS ROUTINE PACKS AND UNPACKS THE INFO NEEDED TO RELATE 

30440CPR THE ELEMENT NAME TO THE RELATIVE POSITION OF THE ELEMENT 

30450CPR IN THE OUTPUT LIST 

30460CPR 

30470CPR 

30480CPR 

30490CPR 

30500CPR CALLING ARGUMENTS — 


30510CPR 

NAME 

ATTRIBUTES 

DEFINITION 

30520CPR 

30530CPR 

MP 

(I/O) 

MATRIX POSITION ARRAY 

30540CPR 

NE 

(I) 

ELM NAME 

30550CPR 

NPOS 

(I/O) 

RELATIVE POSITION IN OUTPUT 

30560CPR 

IRW 

(I) 

SWITCH (0--RE AD , 1--WRITE) 

30570CPR 

30580CPR 

COMMONS USED 

-- NONE 


30590CPR 

30600CPR 

30610CPR FUNCTIONS/ROUTINES CALLED -- 
30620CPR MOD (FORTRAN SUPPLIED) 

30630CPR 

30640CPR FILES USED -- 
30650CPR 

30660CPR 06 S 

30670CPR 

30680CPR LIBRARIES ACCESSED — 

30690CPR NONE 

30700CPR 

30710CPR 

30720CPR LOCAL VARIABLES -- 
30730CPR 

30740CPR SPECIAL COMMENTS ABOUT THIS ROUTINE 
30750CPR 

30760CPR THIS ROUTINE USES DECIMAL PACKING TO PACK 2 POSITIONS 

30770CPR PER WORD OF THE MP ARRAY IT MAKES USE OF THE FORTRAN 

30780CPR SUPPLIED 'MOD ' FUNCTION 
30790CPR THE POSITIONS ARE PACKED AS FOLLOWS — 

30800CPR POSITION 1 »> LEFT HAND SIDE OF WORD 

30810CPR POSITION 2 »> RIGHT HAND SIDE OF WORD 

30820CPR THE MAX POSITION THAT CAN BE PACKED IS NPOS-9999 

30830CPR THE ADVANTAGE OF THIS SCHEME IS THAT THE DIMENSION OF 

30840CPR MP NEED ONLY BE 4000 TO STORE ELMS UP TO 8000 
30850CPR 

30860CPR FOR EXAMPLE 

30870CPR IF THE POSITION OF ELM 1 IS 987 

30880CPR AND THE POSITION OF ELM 2 IS 79 

30890CPR THEN THE VALUE STORED IN MP(1) IS 09870079 

30900CPR NOTE THAT MP(1) CONTAINS THE POSITIONS OF BOTH 

30910CPR ELMS 1 AND 2 

30920CPR 

30930CPR SPECIAL REMARKS/INSTRUCTIONS — 

30940CPR HOST « H6000 

30950CPR 
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30960CPR 

3097 OOP R =es = s * = = = =a = = *— = — — = = 

31620 END 
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3 1 640C ======= B=S==3======a==K S CC=*X===E5=2X======KSSC==S= 

31660 SUBROUTINE SETFLP ( NNEW , NOLD ) 

3 1680C*====*====*===*=========*===*===*==========*===3s:== 

31700CPR 

317 10CPR = = ** = S * = = = = = = = = = = = = * = = = = = = = = = = = = = * = = = = = = = = « = = = = = = = = = = = 
31720CPR 
3 1 730CPR 

31740CPR ROUTINE— SETFLP 

3 1 750CPR 

3 1 760CPR DATE/WRITTEN BY-- 09/19/83 RJ MAFFEO 

31770CPR DATE/REVISED BY — 

3 1 780CPR 

3 1 790CPR FUNCTION/PURPOSE — 

3 1 800CPR 

31810CPR THIS ROUTINE SETS THE H T INPUT FILE TO THE CORRECT 

3 1820CPR POSITION TO READ THE SOLUTION DESIRED 

3 1830CPR 

3 1840CPR 

3 1 850CPR 

31860CPR 

31870CPR 

3 1 880CPR CALLING ARGUMENTS — 

3 1 890CPR NAME ATTRIBUTES DEFINITION 

3 1900CPR 

31910CPR NNEW (I) FILE DESIRED 

3 1 920CPR NOLD (I) CURRENT FILE COUNT 

31930CPR OF THT OUTPUT FILE 

3 1940CPR 

3 1950CPR COMMONS USED — 

3 1960CPR 

31970CPR AFIL 

31980CPR 

3 1 990CPR FUNCTIONS/ROUTINES CALLED — 

32000CPR 

32010CPR TFLSET 

32020CPR HTFPAS 

32030CPR 

32040CPR FILES USED — 

32050CPR 

32060CPR 

32070CPR LIBRARIES ACCESSED — 

32080CPR NONE 
32090CPR 
32 100CPR 

32110CPR LOCAL VARIABLES — 

32 120CPR 

32130CPR SPECIAL COMMENTS ABOUT THIS ROUTINE 

32140CPR 

32150CPR 

32 160CPR SPECIAL REMARKS/INSTRUCTIONS — 

32170CPR HOST * H6000 

32 180CPR 

32190CPR 

3 2 2 00 C P R * * * = * S = = = « = = = = = = = * B S X 2 S S S = B S E E 3 = S S S S = = = = = 3 3 C S S S = E = S X 

32550 END 


106 


GE/TRANCITS 


SETFLP 


PAGE= 29 



TITLE® 


DATE =07/15/85 


TIME = 23 42 


PAGE=30 


325700=== 

32590 

326100== 

32630CPR 

32640CPR 

32650CPR 

32660CPR 

32670CPR 

32680CPR 

32690CPR 

32700CPR 

32710CPR 

32720CPR 

32730CPR 

32740CPR 

32750CPR 

32760CPR 

32770CPR 

32780CPR 

32790CPR 

32800CPR 

32810CPR 

32820CPR 

32830CPR 

32840CPR 

32850CPR 

32860CPR 

32870CPR 

32880CPR 

32890CPR 

32900CPR 

32910CPR 

32920CPR 

32930CPR 

32940CPR 

32950CPR 

32960CPR 

32970CPR 

32980CPR 

32990CPR 

33000CPR 

33010CPR 

33020CPR 

33030CPR 

33040CPR 

33050CPR 

33060CPR 

33070CPR 

33080CPR 

33090CPR 

33 100CPR 

33110CPR 

33120CPR 

33130CPR 

33140CPR 

33150CPR* 

33670 


SUBROUTINE RDTEMP( TEMP , NFLVAL ) 




ROUTINE- 




RDTEMP 


RJ MAFFEO 


DATE/WRITTEN BY— 09/19/83 
DATE/REVISED BY — 

FUNCTION/PURPOSE — 

THIS ROUTINE READS THE TEMPERATURE INFO FOR THE DESIRED 
TIME STEP FROM THE THT OUTPUT FILE AND STORES IT IN 
THE TEMP ARRAY 


CALLING ARGUMENTS — 

NAME ATTRIBUTES 


TEMP 

NFLVAL 


( 0 ) 

(I) 


DEFINITION 

TEMPER STORAGE ARRAY 
CURRENT FILE COUNT 


COMMONS USED — 

AFIL 

SIZE 

CNTLFL 

FUNCTIONS/ROUTINES CALLED — 

HTFIO 

FILES USED — 

IERF S 

LIBRARIES ACCESSED -- 
NONE 

LOCAL VARIABLES — 

SPECIAL COMMENTS ABOUT THIS ROUTINE 


SPECIAL REMARKS/INSTRUCTIONS — 
HOST « H6000 


t*m*B«**Bna*«E = si:t = = «cnKCEa*s = r*saci 


END 
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33690C* 3 = = 

33710 SUBROUTINE CWCRNT(MP , TEMP , TIMVAL . NNODMX , NNODMN , TCRIT , MXNCON) 

337 30C = = = = = e = = = = = = 

33750CPR 

33760CPR=== =============»==============================:==== == 

33770CPR 

33780CPR 

33790CPR ROUTINE-- CWCRNT 

33800CPR 

33810CPR DATE/WRITTEN BY— 09/19/83 RJ MAFFEO 

33820CPR DATE/REVISED BY — 

33830CPR 

33840CPR FUNCTION/PURPOSE — 

33850CPR 

33860CPR THIS ROUTINE COMPUTES A UNIQUE TEMPERATURE FOR EACH 

33870CPR CORNER OF THE HEAT TRANSFER ELM BASED ON THE INDIVIDUAL 

33880CPR TEMPERATURES FOR EACH ELM CONNECTION 

33890CPR 

33900CPR 

33910CPR 

33920CPR 

33930CPR 

33940CPR CALLING ARGUMENTS — 


33950CPR 

33960CPR 

NAME 

ATTRIBUTES 

DEFINITION 

33970CPR 

MP 

(I) 

ELM MATRIX POSITION ARRAY 

33980CPR 

TEMP 

(I) 

ELM TEMPERATURE ARRAY 

33990CPR 

TIMVAL 

(I) 

CURRENT VALUE OF TIME 

34000CPR 

NNODMX 

(I) 

MAX NODE NAME 

34010CPR 

NNODMN 

(I) 

MIN NODE NAME 

34020CPR 

34030CPR 

34040CPR 

TCRIT 

(I) 

TEMPERATURE USED IN MAXTDF TO 
DETERMINE IF CORNER TEMPERATURE 
DIFFERENCE SHOULD BE OUTPUT 

34050CPR 

34060CPR 

34070CPR 

MXNCON 

(I) 

MAX NUMBER CONNECTIONS TO 
ANY NODE 

34080CPR 

34090CPR 

34100CPR 

COMMONS USED 
AFIL 




341 10CPR 

34120CPR FUNCTIONS/ROUTINES CALLED — 

34130CPR 

34140CPR RWCCF 

34 150CPR UNPAKF 

34160CPR MPERW 

34170CPR FCRNTP 

34 180CPR MAXTDF 

34 190CPP WRITMP 

34200CPR FILREW 

34210CPR 

34220CPR FILES USED -- 

34230CPR 

34240CPR 

34250CPR LIBRARIES ACCESSED — 

34260CPR NONE 

34270CPR 

34280CPR 

34290CPR LOCAL VARIABLES -- 
34300CPR 

34310CPR SPECIAL COMMENTS ABOUT THIS ROUTINE 

34320CPR 

34330CPR 

34340CPR SPECIAL REMARKS/INSTRUCTIONS — 

34350CPR HOST « H6000 

34360CPR 

34370CPR 

34380CPR a3SIIes » BC::BS8S = B'csst:B&s:x3::sss::ssse:::sBB3ss>8c 
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35450 


END 
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354 70C = = *’* : = = = == = = = = = =* = = = = = = *= = = = * = = = * = = = * = = = = = = = a= = == = = == = 

35490 SUBROUTINE UNPAKF ( NPACK , NAME , I F ACE ) 

35530CPR 

35540CPR* = = = = = = = = = = = = = s = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = 

35550CPR 

35560CPR 

35570CPR ROUTINE-- UNPAKF 

35580CPR 

35590CPR DATE/WRITTEN BY-- M/29/82 Rd MAFFEO 

35600CPR DATE/REVISED BY-- 

35610CPR 

35620CPR FUNCTION/PURPOSE-- 
35630CPR 

35640CPR THIS ROUTINE UNPACKS THE CORNER-FACE NUMBER INTO AN CORNER 
35650CPR NUMBER ( NAME ) AND THE FACE ARRAY ( I FACE ) 

35660CPR 

35670CPR 

35680CPR 

35690CPR 

35700CPR 

35710CPR CALLING ARGUMENTS — 

35720CPR NAME ATTRIBUTES DEFINITION 

35730CPR 

35740CPR NPACK (I) PACKED CORNER-FACE # 

35750CPR NAME (0) CORNER NAME 

35760CPR I FACE (0) FACE INFO ARRAY 

35770CPR 

35780CPR COMMONS USED -- 
35790CPR 

35800CPR AFIL 

35810CPR 

35820CPR FUNCTIONS/ROUTINES CALLED -- 

35830CPR 

35840CPR 

35850CPR FILES USED -- 
35860CPR 

35870CPR IERF S 

35880CPR 

35890CPR LIBRARIES ACCESSED -- 

35900CPR NONE 

35910CPR 

35920CPR 

35930CPR LOCAL VARIABLES — 

35940CPR 

35950CPR SPECIAL COMMENTS ABOUT THIS ROUTINE 

35960CPR 

35970CPR NOTE PACK CORNER-FACE # = CORNER NAME * 1000000+1 FACPC 
35980CPR WHERE IFACPC IS A SINGLE INTEGER THAT CONTAINS ALL 

35990CPR ALL THE FACE NUMBERS 

36000CPR FOR EXAMPLE 

36010CPR IF THE FACE NUMBERS WERE 6,5, 4, 3, 2,1 

36020CPR THEN I FACPC-65432 1 

36030CPR 

36040CPR 

36050CPR 

36060CPR SPECIAL REMARKS/INSTRUCTIONS — 

36070CPR HOST * H6000 

36080CPR 

36090CPR 

36 100CPR saaxtlKSKXBSZ38aK *E*sa»SRi»B«es«t»«scssiB»aBiBiii 
36480 END 
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36500C = = — = — — — — — • — — — — — = — = = = — = = — — sssssaBzsssxxxssszBBzsssaBSB 

36520 SUBROUTINE FCRNTP(TEMP ,MPOS , IFACE .NCORN , CTP , COF ) 

36540C = *= Z3Za«Z*B3«*ZZB*3 3 3 3 3 3 

36560CPR 

36570CPR* SSZ3333SBSS3S3ES = £S3BEBSSSSS&EB3S = e3B8S3X333SXE33SXZ = 

36580CPR 

36590CPR 

36600CPR ROUTINE— FCRNTP 

36610CPR 

36620CPR DATE/WRITTEN BY-- 11/29/82 RJ MAFFEO 

36630CPR DATE/REVISED BY — 

3S640CPR 

36650CPR FUNCTION/PURPOSE — 

36660CPR 

36670CPR THIS ROUTINE EXTRACTS THE 7 TEMPERATURES ASSC WITH A 

36680CPR CONNECTING ELM IT THEN FINDS THE 4 TEMPS ASSC WITH 

36690CPR THE CORNER OF INTEREST AND DOES THE MULTIPLICATION 

36700CPR TO COMPUTE THE TEMPERATURE FOR THIS CORNER BASED ON 

36710CPR THIS ELM CONNECTION 

36720CPR 

36730CPR 

36740CPR 

36750CPR 

36760CPR 

36770CPR CALLING ARGUMENTS — 


36780CPR 

NAME 

ATTRIBUTES 

DEFINITION 

36790CPR 




36800CPR 

TEMP 

(I) 

ELM TEMPER ARRAY 

36810CPR 

MPOS 

(I) 

MATRIX POS OF THE 

36820CPR 

IFACE 

(I) 

ELM FACE ARRAY 

36830CPR 

NCORN 

(I) 

CORNER NUMBER 

36840CPR 

CTP 

(0) 

CORNER TEMPERATURE 

36850CPR 

COF 

(I) 

CORNER COEFF ARRAY 

36860CPR 




36870CPR 

COMMONS USED 

— 


36880CPR 




36890CPR 

AFIL 




36900CPR 

36910CPR FUNCTIONS/ROUTINES CALLED — 

36920CPR 

36930CPR 

36940CPR FILES USED — 

36950CPR 

36960CPR 

36970CPR LIBRARIES ACCESSED — 

36980CPR NONE 

36990CPR 

37000CPR 

37010CPR LOCAL VARIABLES — 

37020CPR 

37030CPR SPECIAL COMMENTS ABOUT THIS ROUTINE 

37040CPR 

37050CPR 

37060CPR SPECIAL REMARKS/INSTRUCTIONS — 

37070CPR HOST * H6000 

37080CPR 

37090CPR 

37 lOOCPR------------------------- ------- ------ 

37970 END 


GE/TRANCITS 


FCRNTP 


PAGE* 34 


111 


TITLE' DATE =07/15/85 TIME = 23 42 PAGE=35 

38010 SUBROUTINE MAXTDF ( NEL , CRNT , NCON , NDONE . TIMVAL , TCRIT , IREC ) 

38050CPR 

3 8 06 OC P R s * = — = = = = = = = * = = = = = = = = a = = = = = = = = = = K = = = = s = = = = 

38070CPR 

38080CPR 

38090CPR ROUTINE-- MAXTDF 

38 100CPR 

38110CPR DATE/WRITTEN BV-- 11/30/82 RJ MAFFEO 

38 120CPR DATE/REVISED BY-- 

38130CPR 

38 140CPR FUNCTION/PURPOSE — 

38150CPR 

38 160CPR THIS ROUTINE COMPUTES THE MAX DIFFERENCE BETWEEN 

38170CPR ALL THE POSSIBLE CORNER TEMPERATURES AT A CORNER 

38 180CPR AND FLAGS IT IF IT IS GREATER THAN TCRIT 

38 190CPR 

38200CPR 

38210CPR 

38220CPR 

38230CPR 

38240CPR CALLING ARGUMENTS — 


38250CPR 

NAME 

ATTRIBUTES 

DEFINITION 

38260CPR 




38270CPR 

NEL 

(I) 

ARRAY OF ELMS CONNECTED TO THIS NODE 

38280CPR 

CRNT 

(I) 

ARRAY OF TEMPS FOR THIS NODE 

38290CPR 

NCON 

(I) 

NUMBER OF ELM CONNECTED TO THIS NODE 

38300CPR 

NDONE 

(I) 

PAGE HEADER COUNTER 

38310CPR 



( *0 WRITE HEADER, NOT = 0 DONT WRITE) 

38320CPR 

TIMVAL 

(I) 

CURRENT VALUE OF TIME 

38330CPR 

TCRIT 

(I) 

CRITICAL TEMP FOR PRINTOUT 

38340CPR 



IF MAX DIFF < TCRIT -- DONT WRITE 

38350CPR 



IF MAX DIFF > TCRIT -- WRITE 

38360CPR 

IREC 

(I) 

RECORD NUMBER ( NODE NUMBER )0F CORNER 

38370CPR 




38380CPR 

COMMONS USED 

— 


38390CPR 




38400CPR 

AFIL 




38410CPR 

38420CPR FUNCTIONS/ROUTINES CALLED -- 
38430CPR 

38440CPR SORT 

38450CPR 

38460CPR FILES USED — 

38470CPR 

38480CPR IN0T6 S 

38490CPR 

38500CPR LIBRARIES ACCESSED -- 

38510CPR NONE 

38520CPR 

38530CPR 

38540CPR LOCAL VARIABLES — 

38550CPR 

38560CPR SPECIAL COMMENTS ABOUT THIS ROUTINE 
38570CPR 


38580CPR 

38590CPR SPECIAL REMARKS/INSTRUCTIONS — 

38600CPR HOST * H6000 

38610CPR 

38620CPR 

38630CPR** «*«***nsEi**iEiia»iisr**aEEiiEiB£!!B = *ai8 = i::3rE = t; 

39160 END 
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391 80C ******* XHKK *KK«KKKSS8ttK3mK*KSXSSCXBKSaCS = X3S = SX2X = = X = 

39200 SUBROUTINE WRITMP ( TBAR , TIMVAL , NOD , NDONE ) 

39220C SS * SSS *8"Sea»8SS8SIB«IS*BBBS8S3ZS8SS«aSS3&CSSSfSSB8 

39240CPR 

39250CPR* * * K = BE = = = -* = - = ** = * = = = »s = = ***Bssc3S!s= = s = as*ss*B = = = 2:as 

39260CPR 

39270CPR 

39280CPR ROUTINE-- WRITMP 

39290CPR 

39300CPR DATE/WRITTEN BY— 11/30/82 RJMAFFEO 

39310CPR DATE/REVISED BY-- 

39320CPR 

39330CPR FUNCTION/PURPOSE-- 
39340CPR 

39350CPR THIS ROUTINE WRITES THE VALUES OF THE TEMPERATURES AT THE 

39360CPR STRESS POINTS TO HARD COPY AND TO THE OUTPUT FILE 

39370CPR 

39380CPR 

39390CPR 

39400CPR 

39410CPR 

39420CPR CALLING ARGUMENTS — 


39430CPR 

NAME 

ATTRIBUTES 

DEFINITION 

39440CPR 

39450CPR 

TBAR 

(I) 

TEMPERATURE AT THE STRESS POINT 

39460CPR 

TIMVAL 

(I) 

CURRENT VALUE OF TIME 

39470CPR 

NOD 

(I) 

STRESS POINT NAME ( NODE #) 

39480CPR 

NDONE 

(I) 

SWITCH 

39490CPR 

39500CPR 

39510CPR 

39520CPR 

39530CPR 

39540CPR 

39550CPR 

COMMONS USED 
AFIL 

— 

■0 — WRITER PAGE HEADER 

NOT 0 -- DO NOT WRITE HEADER 

»-1 -- FLUSH BUFFERS 


39560CPR 

39570CPR FUNCTIONS/ROUTINES CALLED — 

39580CPR 

39590CPR 

39600CPR FILES USED -- 
39610CPR 

39620CPR IN0T6 S 

39630CPR IOUTF S 

39640CPR 

39650CPR LIBRARIES ACCESSED -- 

39660CPR NONE 

39670CPR 

39680CPR 

39690CPR LOCAL VARIABLES — 

39700CPR 

39710CPR SPECIAL COMMENTS ABOUT THIS ROUTINE 

39720CPR 

39730CPR 

39740CPR SPECIAL REMARKS/INSTRUCTIONS — 

39750CPR HOST = H6000 

39760CPR 

39770CPR 

39780CPR» — 

40580 END 
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40600C 3E3SSBBSSSStESSSSSSS3SSSESBSS = SSSSS3SSaEBSSSSSSSSSSS = 

40620 SUBROUTINE CINV(CD , XSP , YSP , ZSP , XSF ,NE . JCN) 

40640C :sss3s3s3ssssssssssssassssssssssssasss«s=sssssssasrs 

40660CPR 

40670CPR 5 sissssssssssssssesetsssesssssssssssssessssssssss 

40680CPR 

40690CPR 

40700CPR ROUTINE-- CINV 
40710CPR 

40720CPR DATE/WRITTEN BY-- 9/10/82 RJ MAFFEO 

40730CPR DATE/REVISED BY-- 01/12/83 RJ MAFFEO 
40740CPR 

40750CPR FUNCTION/PURPOSE -- 

40760CPR THIS SUB FIND THE VALUES OF THE SHAPE FUNCTIONS FOR 
40770CPR ANY POINT NEAR A TETRAHEDRAL ELEMENT THE SHAPE FUNCTIONS 

40780CPR ARE BASED ON THE SIMPLEX TETRAHEDRAl ELEMENT 

40790CPR IT FIRST FINDS THE INVERSE OF THE COORDINATE MATRIX (-C-) 

40800CPR IT THEN THENS MULTIPLIES- THIS INVERSE INTO THE COORD 

40810CPR VECTOR (1 X Y Z) TO FIND THE SHAPE FUNCTIONS 

4C820CPR THIS ROUTINE CAN BE USED TO COMPUTE THE WEIGHTING 

40830CPR COEFF FOR THE CORNERS OF A THTD ELM BASED ON THE TEMPERATURES 

40840CPR AT THE CENTER OF THE FACES AND THE CENTROID 

40850CPR 

40860CPR CALLING ARGUMENTS-” 


40870CPR 

NAME 

ATTRIBUTES 

DEFINITION 

40880CPR 

CD 

m 

12X1 ARRAY 

COORDINATE ARRAY 

40890CPR 




X , Y , Z COORDS OF CORNERS OF TETRA 

40900CPR 

XSP 

(i) 



40910CPR 

YSP 

(i) 



40920CPR 

ZSP 

(XI 


COORDS OF THE POINT FOR WHICH 

40930CPR 




THE SHAPE FUNCTIONS ARE DESIRED 

40940CPR 

XSF 

(0) 

5X1 ARRAY 

VECTOR OF SHAPE FUNCTIONS 

40950CPR 

NE 

(I) 


ELMEMENT NAME 

40960CPR 

JCN 

(I) 


CORNER NUMBER 


40970CPR 

40980CPR 

40990CPR 

41000CPR FILES USED-- 
4 1010CPR AFIL 

4 1020CPR 

4 1030CPR COMMONS USED — 

41040CPR NONE 

4 1050CPR 
41060CPR 

4 1070CPR ROUTINES CALLED — 

4 1080CPR NONE 

4 1090CPR 

41100CPR LIBRARIES ACCESSED-- 
41110CPR NONE 

41 120CPR 

4 1 130CPR LOCAL VARI ABLES-- 

41140CPR I NX 18X1 THIS VECTOR IS USED TO SELECT THE 

41150CPR CORRECT INDICES OF THE C MATRIX FOR 

4 1 160CPR EACH ELEMENT OF THE C INVERSE MATRIX 

41170CPR Cl 4X4 INVERSE OF THE COORD MATRIX(-C-) 

4 1 180CPR C 5X4 2D MATRIX CONTAINING CORRDS OF TETRA 

41 190CPR 
4 1200CPR 

4 1210CPR SPECIAL COMMENTS ABOUT THIS ROUTINE — 

4 1220CPR 

4 1230CPR NOTE THAT THE COORD MATRIX(-C-) IS AUGMENTED WITH 

4 1240CPR A ROW OF ONES THIS ROW OF ONES IS USED IN THE COMPUTATION 

41250CPR OF THE INVERSE 

4 1260CPR 

4 1270CPR NOTE XSF(5)«1/DIST 

4 1280CPR DIST*DISTANCE FROM CENTROID TO CORNER 

4 1290CPR 
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4 1300CPR NOTE A LOCAL COORO SYSTEM IS ESTABLISHED FOR EACH ELM 
41310CPR THIS SYSTEM IS FORMED AS 0(I)-QC 

41320CPR WHERE 0(1) IS THE COORDS OF THE VERTICES 

41330CPR AND OC IS THE COORDS OF THE CENTROID OF THE ELM 

4 1340CPR 
41350CPR 
4 1360CPR 

4 1370CPR SPECIAL REMARKS/INSTRUCTIONS— 

41380CPR HOST * H6000 

41390CPR 

4 1400CPR* *************** * **************************************************** 
42720 END 
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42740CC* * 

42760 SUBROUTINE PRNTIM( ISTOP , 1 1 , TB , TE ) 

42930 END 
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42950C sss:ies:sss:&& & 

4297 0^ SUBROUTINE EMINMX ( XC . YC , 2C , CADD . XYZV , NEMNMX , NE t NNPE ) 

42990C* SSS «* = S = SSS 

43010CPR 

43020CPR 58 

43030CPR 

43040CPR 

43050CPR ROUTINE— EMINMX 

43060CPR 

43070CPR DATE/WRITTEN BY— 04/05/83 RJ MAFFEO 

43080CPR DATE/REVISED BY — 

43090CPR 

43100CPR FUNCTION/PURPOSE— 

43110CPR 

43 120CPR THIS ROUTINE FINDS THE MIN AND MAX XYZ COORDS 
43130CPR OF THE HEAT TRANSFER ELEMENT AND STORES THIS INFO 
43 140CPR INTO THE XYZV ARRAY 

43150CPR IT ALSO KEEPS TRACK OF THE ELMS WITH THE SMALLEST 
43 160CPR AND LARGEST XYZ COORDS AND STORES THEM IN THE 

43 1 70CPR NEMNMX ARRAY 

43 180CPR 
43190CPR 

43200CPR CALLING ARGUMENTS — 

43210CPR NAME ATTRIBUTES DEFINITION 

43220CPR 

43230CPR XC (I) » 

43240CPR YC (I) »» COORD ARRAYS FOR HT ELM VERTICES 

43250CPR ZC (I) » 

43260CPR CADD (I) CONSTANT COORD ADDER 

43290CPR NE (I) CURRENT ELM NAME 

43300CPR NNPE (I) * OF NODES PER ELM 

43310CPR 

43320CPR COMMONS USED — 

43330CPR 

43340CPR 

43350CPR FUNCTIONS/ROUTINES CALLED — 

43360CPR 

43370CPR SORT 

43380CPR 

43390CPR FILES USED — 

43400CPR 

43410CPR 

43420CPR LIBRARIES ACCESSED — 

43430CPR NONE 

43440CPR 

43450CPR 

43460CPR LOCAL VARIABLES — 

43470CPR 

43480CPR SPECIAL COMMENTS ABOUT THIS ROUTINE 

43490CPR 

43500CPR 

43510CPR SPECIAL REMARKS/INSTRUCTIONS — 

43520CPR HOST • H6000 

43530CPR 

43540CPR 

43550CPR»*»»»»»>— 

43910 END 


GE/TRANCITS 


EMINMX 


PAGE* 40 


117 



TITLE* DATE =07/15/85 TIME = 23 42 PAGE=41 

4 39 30C a 3 = = = ■* S s* = = = = = = = X s B s = = = = = = = = = = a = = :=5SS.S a a = = = = = = = = = = s 

43950 SUBROUTINE FTIME ( A 1 ) 

43990CPR 

4 4 OOOC PR = I3B - = - = :: = 34 = = = = = = = = = s3 = :s:smz = = :;n = :r = = = = s£;; = = :r 

44010CPR 

44020CPR 

44030CPR ROUTINE-- FTIME 

44040CPR 

44050CPR DATE/WRITTEN BY-- 04/05/83 RJ MAFFEO 

44060CPR DATE/REVISED BY-- 

44070CPR 

44080CPR FUNCTION/PURPOSE-- 
44090CPR 

44 100CPR THIS ROUTINE CALLS THE APPROPRIATE SYSTEM ROUTINE 

44110CPR TO MEASURE THE PROCESSOR TIME 

44120CPR 

44130CPR 

44140CPR CALLING ARGUMENTS — 

44150CPR NAME ATTRIBUTES DEFINITION 

44 160CPR 

44170CPR A 1 (0) PROCESSOR TIME 

44180CPR 

44190CPR COMMONS USED -- NONE 

44200CPR 

44210CPR 

44220CPR FUNCTIONS/ROUTINES CALLED -- 
44230CPR 

44240CPR PTIME -- HONEY WELL TIMER 

44250CPR 

44260CPR FILES USED -- 

44270CPR 

44280CPR 

44290CPR LIBRARIES ACCESSED -- 

44300CPR NONE 

44310CPR 

44320CPR 

44330CPR LOCAL VARIABLES — 

44340CPR 

44350CPR SPECIAL COMMENTS ABOUT THIS ROUTINE 

44360CPR 

44370CPR 

44380CPR SPECIAL REMARKS/INSTRUCTIONS — 

44390CPR HOST = H6000 

44400CPR 

44410CPR 

44420CPR 8BS3ls = SIS *»:tsE8smiBEn33iB::szziE: = 3st£: = :ttt:: = : 

44470 END 
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44490C = I = ==c = = sa = sae:»*:sE = »3 = = = isi:s3 = 3i= 5 ea3ttr = r= =3 r5se 
44510 SUBROUTINE SSURCH 

44530 C ssa * sssssssa sssss 3 B«**iissisM 5 seasess:«*s 8 »* 5 sss*»s 

44550CPR 

44560CPR a » 8 «. s ,. 3 »,«. e8 , s „ ss , 3SBSS , i;1 . sseBSSS , e , sss „ ssssce 

44570CPR 

44580CPR 

44590CPR ROUTINE— SSURCH 

44600CPR 

44610CPR DATE/WRITTEN BY— 04/05/83 RJ MAFFEO 

44620CPR DATE/REVISED BY — 

44630CPR 

44640CPR FUNCTION/PURPOSE — 

44650CPR 

44660CPR THIS ROUTINE IS THE MINI-EXEC FOR THE STRESS POINT 

44670CPR TO HEAT TRANSFER ELM SEARCH ROUTINES IT CALLS ALL 

44680CPR REQUIRED ROUTINES TO FIND WHICH HT ELM CONTAINS 

44690CPR EACH STRESS POINT AND THEN FINDS THE WEIGHTING COEFFS 

44700CPR FOR THAT STRESS POINT 

44710CPR 

44720CPR 

44730CPR 

44740CPR CALLING ARGUMENTS — 

44750CPR NAME ATTRIBUTES DEFINITION 

44760CPR 

44770CPR 

44780CPR COMMONS USED — 

44790CPR 

44800CPR AFIL 

44810CPR SIZE 

44820CPR SPDAT 

44830CPR 

44840CPR FUNCTIONS/ROUTINES CALLED — 

44850CPR 

44860CPR KNTSPF 

44870CPR RFILCR 

44880CPR ZERORF 

44890CPR GNLSUR 

44900CPR 

44910CPR FILES USED — 

44920CPR 

44930CPR 

44940CPR LIBRARIES ACCESSED — 

44950CPR NONE 

44960CPR 

44970CPR 

44980CPR LOCAL VARIABLES — 

44990CPR 

45000CPR SPECIAL COMMENTS ABOUT THIS ROUTINE 

45010CPR 

45020CPR 

45030CPR SPECIAL REMARKS/INSTRUCTIONS — 

45040CPR HOST « H6000 

45050CPR 
45060CPR 
4 507 OCPR = 

45600 END 
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45620C SSS **- ▼VSSS = SSSSSS = 3SSS3S2tEESSS3SS2SSBS33SSXSSSS = SSSS 

45640 SUBROUTINE KNTSPF (NUMSP , MXNSP ) 

456G0C======================= : ============*===== =========== 

45680CPR 

45690CPR* ss = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = x = * = r = * = = a = = e = = = = = 

45700CPR 

45710CPR 

45720CPR ROUTINE-- KNTSPF 

45730CPR 

45740CPR DATE/WRITTEN BY— 04/05/83 RJ MAFFEO 

45750CPR DATE/REVISED BY — 

45760CPR 

45770CPR FUNCTION/PURPOSE— 

45780CPR 

45790CPR THIS ROUTINE COUNTS THE NUMBER OF STRESS POINTS 

45800CPR TO BE INTERFACED FROM THE STRESS NODE FILE AND 

45810CPR ALSO FINDS THE LARGEST STRESS NODE NAME 

45820CPR 

45830CPR 

45840CPR 

45850CPR CALLING ARGUMENTS — 

45860CPR NAME ATTRIBUTES DEFINITION 

45870CPR 

45880CPR NUMSP (0) * OF STRESS NODES 

45890CPR MXNSP (0) LARGEST STRESS NODE NAME 

45900CPR 

45910CPR COMMONS USED — 

45920CPR 

45930CPR AFIL 

45940CPR 

45950CPR FUNCTIONS/ROUTINES CALLED — 

45960CPR 

45970CPR FILREW 

45980CPR RDSPF 

45990CPR 

46000CPR FILES USED — 

46010CPR 

46020CPR I ERF S 

46030CPR 

46040CPR LIBRARIES ACCESSED — 

46050CPR NONE 

46060CPR 

46070CPR 

46080CPR LOCAL VARIABLES — 

46090CPR 

46 100CPR SPECIAL COMMENTS ABOUT THIS ROUTINE 

46110CPR 

46 120CPR 

46130CPR SPECIAL REMARKS/INSTRUCTIONS — 

46 140CPR HOST « H6000 

46 150CPR 

46160CPR 

461 70CPR**®*« aiiiKrn*nirE***:nn=nn=2Ec:cca*Ec*=32E*3=i=3i3* 

46560 END 
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46580C ==SCBBSBSaB = s * ssss =***:BsrisaBssrsa=ssss==sr===aa=sB= 

46600 SUBROUTINE RFILCR 

46620 C Sa * SS> 8 S 8 SSBSSaa S 3 Sa&S 8 SSSSICSSSSBSBBSB 8 IBSS:StB 3 SSr 

46640CPR 

46G50CPR ===3=a * 3SS==es * K==SE= s=5ass*8sssr sissBsaasassasasssss 

46660CPR 

46670CPR 

46680CPR ROUTINE-- RFILCR 

46690CPR 

46700CPR DATE/WRITTEN BY-- 04/06/83 RJ MAFFEO 

46710CPR DATE/REVISED BY-- 

46720CPR 

46730CPR FUNCT I ON/PURPOSE -- 
46740CPR 

46750CPR THIS ROUTINE SIZES ALL THE RANDOM FILES 

46760CPR IT CALLS THE APPROPRIATE SYSTEM ROUTINE TO DEFINE 

46770CPR THE WIDTH OF THE FILE 

46780CPR 

46790CPR 

46800CPR 

46810CPR CALLING ARGUMENTS — 

46820CPR NAME ATTRIBUTES DEFINITION 

46830CPR 

46840CPR 

46850CPR COMMONS USED — 

46860CPR 

46870CPR 

46880CPR FUNCTIONS/ROUTINES CALLED -- 
46890CPR 

46900CPR RANSIZ -- HONEYWELL ROUTINE TO SET WIDTH OF FILE 

46910CPR DEFINE FILE — IBM ROUTINE TO SET WIDTH 

46920CPR 

46930CPR FILES USED — 

46940CPR 

46950CPR 

46960CPR LIBRARIES ACCESSED -- 

46970CPR NONE 

46980CPR 

46990CPR 

47000CPR LOCAL VARIABLES -- 
47010CPR 

47020CPR SPECIAL COMMENTS ABOUT THIS ROUTINE 
47030CPR ' 

47040CPR 

47050CPR SPECIAL REMARKS/INSTRUCTIONS — 

47060CPR HOST « H6000 

47070CPR 

47080CPR 

47090 CPR 3 B 3 SBtEC 3 *<* 333 >e 3 cssBt 8 EBBsi:a=ssBtese 3 * 

47800 END 
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47820C - = a - ac = s = s3 = = ss SS = S3 = SXXES3 = S = 3 = SS = ZSS = SSSXCSS = = = = 3 = = 
47840 SUBROUTINE ZERORF ( I FC , MXN , NWD , A ) 

4 7 8 60C * e = = = = = = s = s = = = s = = s = is = = = = = s = 3sst: = = = sss s = = = = = = = = = = = s = 

47880CPR 

47890CPR B ***S = *S*=S«**SS*SSSS«**SSS*»S*»=»SSB88BS»»BBSSSSSSSSSBB 

47900CPR 

47910CPR 

47920CPR ROUTINE-- ZERORF 

47930CPR 

47940CPR DATE/WRITTEN BY — 04/07/83 RJ MAFFEQ 

47950CPR DATE/REVISED BY-- 

47960CPR 

47970CPR FUNCT I ON/ PURPOSE -- 
47980CPR 

47990CPR THIS ROUTINE INITIALIZES RANDOM FILES 

48000CPR 

48010CPR 

48020CPR CALLING ARGUMENTS — 


48030CPR 

48040CPR 

NAME 

ATTRIBUTES DEFINITION 

48050CPR 

I FC 

(I) 

FILECODE 

48060CPR 

MXN 

(I) 

MAX RECORD NUMBER 

48070CPR 

NWD 

(I) 

* OF WORDS IN -A- ARRAY 

48080CPR 

48090CPR 

48100CPR 

A 

COMMONS USED 

(I) 

INITIALIZEING ARRAY 


481 10CPR 
48 120CPR 

48130CPR FUNCTIONS/ROUTINES CALLED -- 

48140CPR 

48150CPR 

48 160CPR FILES USED -- 
48170CPR 

48 180CPR I FC R 

48190CPR 

48200CPR LIBRARIES ACCESSED — 

48210CPR NONE 

48220CPR 

48230CPR 

48240CPR LOCAL VARIABLES — 

48250CPR 

48260CPR SPECIAL COMMENTS ABOUT THIS ROUTINE 

48270CPR 

48280CPR 

48290CPR SPECIAL REMARKS/INSTRUCTIONS — 

48300CPR HOST - H6000 

48310CPR 

48320CPR 

48330CPR****®****-** 1 -* = * = *=«* BISBSSSBSSESaSStSIISSSESSSeaaSCS 

48390 END 
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484 10C= 


48430 
484500*== 
48470CPR 
48480CPR* 
48490CPR 
48500CPR 
48510CPR 
48520CPR 
48530CPR 
48540CPR 
48550C D R 
48560CPR 
48570CPR 
48580CPR 
48590CPR 
48600CPR 
48610CPR 
48620CPR 
48630CPR 
48640CPR 
48650CPR 
48660CPR 
48670CPR 
48680CPR 
48690CPR 
48700CPR 
48710CPR 
48720CPR 
48730CPR 
48740CPR 
48750CPR 
48760CPR 
48770CPR 
48780CPR 
48790CPR 
48800CPR 
48810CPR 
48820CPR 
48830CPR 
48840CPR 
48850CPR 
48860CPR 
48870CPR 
48880CPR 
48890CPR 
48900CPR 
48910CPR 
48920CPR 
48930CPR 
48940CPR 
48950CPR 
48960CPR 
48970CPR* 
• 49580 


SUBROUTINE GNLSUR 


ROUTINE— GNLSUR 

DATE/WRITTEN BY-- 04/07/83 RJ MAFFEO 

DATE/REVISED BY- 

FUNCTION/ PURPOSE — 

THIS ROUTINE IS CALL BY THE SEARCH MINI-EXEC (SSURCH) 
IT INITIATES THE SEARCH THRU THE HEAT TRANSFER ELMS 
TO FIND WHICH ELM CONTAINS EACH STRESS POINT 


CALLING ARGUMENTS — 

NAME ATTRIBUTES DEFINITION 


COMMONS USED — 

AFIL 

SIZE 

WORK 

FUNCTIONS/ROUTINES CALLED — 

EINCOR 

RDSPF 

FNDELM 

FILES USED — 


IN0T6 S 

LIBRARIES ACCESSED -- 
NONE 


LOCAL VARIABLES — 

SPECIAL COMMENTS ABOUT THIS ROUTINE 

SPECIAL REMARKS/INSTRUCTIONS — 

HOST * H6000 


*»*«****»***«**»****»*S*** = *SB88B* = *B8BBtSEBBB**eBBa 

END 
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496000************= * =** 





49620 

49640C*** 

49660CPR 

49670CPR* 

49680CPR 

49690CPR 

49700CPR 

49710CPR 

49720CPR 

49730CPR 

49740CPR 

49750CPR 

49760CPR 

49770CPR 

49780CPR 

49790CPR 

49800CPR 

49810CPR 

49820CPR 

49830CPR 

49840CPR 

49850CPR 

49860CPR 

4987CCPR 

49880CPR 

49890CPR 

49900CPR 

49910CPR 

49920CPR 

49930CPR 

49940CPR 

49950CPR 

49960CPR 

49970CPR 

49980CPR 

49990CPR 

50000CPR 

50010CPR 


SUBROUTINE RDSPF (NN. X , Y f Z , IERC ) 


ROUTINE— RDSPF 

DATE/WRITTEN BY— 04/07/83 RU MAFFEO 
DATE/REVISED BY- 

FUNCTION/ PURPOSE— 

THIS ROUTINE READS THE STRESS POINT COORD INFO 
FROM THE STRESS POINT NODE FILE 


CALLING ARGUMENTS — 

NAME ATTRIBUTES DEFINITION 


NN (0) 

X (0) 

Y (0) 

2 ( 0 ) 

IERC (0) 


NODE NAME 
X COORD 
Y COORD 
Z COORD 
ERROR CODES 

0 — NO ERRORS 

1 — END OF FILE 

2 — READ ERROR 


COMMONS USED -- 


AFIL 


FUNCTIONS/ROUTINES CALLED 


FILES USED — 


50020CPR 


50030CPR 

50040CPR 

50050CPR 

5006 OCPR 

50070CPR 

50080CPR 

50090CPR 

50100CPR 

501 10CPR 

50120CPR 

50130CPR 

50140CPR 

50150CPR 

50160CPR 

50170CPR 

50180CPR* 

50440 


ISPNF S 

LIBRARIES ACCESSED — 
NONE 


LOCAL VARIABLES — 

SPECIAL COMMENTS ABOUT THIS ROUTINE 

SPECIAL REMARKS/INSTRUCTIONS — 

HOST * H6000 


END 
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50^60C== =3S32E2SSB=SX====22==SKS====S=SS=S==X2SC==SS==SCSS 

50480 SUBROUTINE FILREW(IFC) 

50500C ** = **s: = ** = = * = = = = = = = s* = = = = = = s = = = = = = ss=ec = *: = = = = = = 2 = 3 = = 2 
50520CPR 

50530CPR** = 33SC = * = * S * = * aIS ** = ** = = ** = = * = = = « = = « = = = * = = *2 = = = = B222 = = 

50540CPR 

50550CPR 

50560CPR ROUTINE-- FILREW 

50570CPR 

50580CPR DATE/WRITTEN BY-- 04/07/83 RJ MAFFEO 

50590CPR DATE/REVISED BY — 

50600CPR 

50610CPR FUNCTION/PURPOSE-- 
50620CPR 

50630CPR THIS ROUTINE REWINDS FILES 

50640CPR 

50650CPR 

50660CPR CALLING ARGUMENTS — 


50670CPR 

50680CPR 

NAME 

ATTRIBUTES 

DEFINITION 

50690CPR 

50700CPR 

IFC 

(I) 

FILE CODE 

50710CPR 

COMMONS USED 

-- 



50720CPR 

50730CPR 

50740CPR FUNCTIONS/ROUTINES CALLED -- 
50750CPR 

50760CPR FILES USED -- 

50770CPR 

50780CPR 

50790CPR LIBRARIES ACCESSED -- 

50800CPR NONE 

50810CPR 

50820CPR 

50830CPR LOCAL VARIABLES — 

50840CPR 

50850CPR SPECIAL COMMENTS ABOUT THIS ROUTINE 

50860CPR 

50870CPR 

50880CPR SPECIAL REMARKS/INSTRUCTIONS — 

50890CPR HOST * H6000 

50900CPR 

50910CPR 

50920CPR* = 

50960 END 
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509800* 
51000 
510200=* 
5 1040CPR 
5 1050CPR 
5 1060CPR 
5 1070CPR 
5 1080CPR 
5 1090CPR 
51 100CPR 
511 10CPR 
51 120CPR 
51 130CPR 
51 140CPR 
51 160CPR 
51170CPR 
51 180CPR 
51 190CPR 
5 1200CPR 
512 IOC PR 
5 1220CPR 
5 1 230CPR 
5 1 240CPR 
5 1 250CPR 
5 1 260CPR 
51270CPR 
5 1280CPR 
5 1 290CPR 
5 1 300CPR 
51310CPR 
5 1 320CPR 
5 1 330CPR 
5 1 340CPR 
5 1 350CPR 
5 1 360CPR 
51370CPR 
5 1 380CPR 
5 1 390CPR 
51400CPR 
51410CPR 
5 1420CPR 
5 1 430CPR 
5 1440CPR 
5 1450CPR 
5 1460CPR 
5 1470CPR 
5 1480CPR 
5 1490CPR 
5 1 500CPR 
51510CPR 
5 1 520CPR 
5 1 530CPR 
51540CPR 
5 1 550CPR 
5 1 560CPR 
5 1 570CPR 
51580CPR 
51590CPR- 
52000 


SUBROUTINE EINCOR( ENX , NELMN, NELMX ) 




ssasBs:: 


ROUTINE-- EINCOR 

DATE/WRITTEN BY-- 04/07/83 RJ MAFFEO 

DATE/REVISED BY-- 

FUNCTIGN/PURPQSE-- 

COORDS INTO CORE FROM THE RANDOM ELM FILE 

THE DATA STORED IN CORE IS — 

ELM NUMBER 

XMIN 

YMIN 

ZMIN 

XMAX 

YMAX 

ZMAX 


CALLING ARGUMENTS — 

NAME ATTRIBUTES DEFINITION 

ENX (0) ARRAY CONTAINING WINDOW INFO 


NELMN 

(i) 

MIN 

ELM 

NAME 

NELMX 

(i) 

MAX 

ELM 

NAME 

COMMONS USED 

— 





AFIL 

ELMDAT 

FUNCTIONS/ROUTINES CALLED -- 
RWELMS 

FILES USED -- 


LIBRARIES ACCESSED — 
NONE 


LOCAL VARIABLES — 

SPECIAL COMMENTS ABOUT THIS ROUTINE 

SPECIAL REMARKS/INSTRUCTIONS — 

HOST = H6000 

= SS ** = = = ***SS****ZSSBSSSSS2*K = KSS8JB«BS8SKXSBKBB = S 

END 
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5 2020C 

52040 SUBROUTINE FNDELM(N$P . XP . YP , ZP , I FOND ) 

52060C==*==«*=====**====*===*=======««=*==*============= 

52080CPR 

52090CPR»»=**»»**=========*==»*»===-===«====*»=**«= ========== 

52100CPR 
521 10CPR 

52120CPR ROUTINE — FNDELM 

52130CPR 

52140CPR DATE/WRITTEN BY-- 04/07/83 RU MAFFEO 

52 150CPR DATE/REVISED BY — 

52 160CPR 

52170CPR FUNCTION/PURPOSE— 

52 180CPR 

52190CPR THIS ROUTINE DOES THE ACTUAL SEARCH THRU THE HEAT 

52200CPR TRANSFER ELMS TO FIND WHICH ELM CONTAINS THE STRESS 

52210CPR POINT 

52220CPR IT CALLS THE REQUIRED ROUTINES TO GET THE ELM 

52230CPR COORDINATES, CHECKS TO SEE IF THE STRESS POINT 

52240CPR IS WITHIN THE ELM AND IF IT IS COMPUTES 

52250CPR AND WRITES THE WEIGHTING COEFFS FOR THIS STRESS 

52260CPR POINT WITH RESPECT TO THE VERTICES OF THE 

52270CPR HEAT TRANSFER ELEMENT 

52280CPR 

52290CPR 

52300CPR 

52310CPR CALLING ARGUMENTS — 


52320CPR 

52330CPR 

NAME 

ATTRIBUTES 

DEFINITION 

52340CPR 

NSP 

(I) 

STRESS NODE NAME 

52350CPR 

XP 

(I) » 


52360CPR 

YP 

(I) »» 

STRESS POINT COORDS 

52370CPR 

ZP 

(I) » 


52380CPR 

I FOND 

(0) 

INDICATOR 


52390CPR 0 — DID NOT FIND CONTAINING ELM 

52400CPR 1 -- DID FIND CONTAINING ELM 

52410CPR 

52420CPR COMMONS USED — 

52430CPR 

52440CPR AFIL 

52450CPR ELMDAT 

52460CPR SIZE 

52470CPR SPDAT 

52480CPR 

52490CPR FUNCTIONS/ROUTINES CALLED -- 
52500CPR 

52510CPR CKWIND 

52520CPR GETWND 

52530CPR I3DSF 

52540CPR RWSPCF 

52550CPR 

52560CPR FILES USED — 

52570CPR 

52580CPR 

52590CPR LIBRARIES ACCESSED -- 

52600CPR NONE 

52610CPR 

52620CPR 

52630CPR LOCAL VARIABLES — 

52640CPR 

52650CPR SPECIAL COMMENTS ABOUT THIS ROUTINE 

52660CPR 

52670CPR 

52680CPR SPECIAL REMARKS/INSTRUCTIONS — 

52690CPR HOST « H6000 

52700CPR 

52710CPR 
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54420 END 
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544 40C * B * = s = aB *** ass * iSBE ss2=ssssa**sBSs8srB*rr*3*ESSSsssss 
54460 SUBROUTINE GETWND( IDONE. IWHR) 

54480C SttBBS * :3;3SS =SS3SSBSSSrBSSSS*tSSBSSSSS 

54500CPR 

545 10CPR 5 = a = = = B *»sssi* = »aetxss3szss«sassBtatB«csBss = *ss 

54520CPR 

54530CPR 

54540CPR ROUTINE-- GETWND 

54550CPR 

54560CPR DATE/WRITTEN BY-- 04/07/83 RJ MAFFEO 

54570CPR DATE/REVISED BY-- 

54580CPR 

54590CPR FUNCTI ON/PURPOSE- - 

54600CPR 

54620CPR OF THE HEAT TRANSFER ELMS AND THE COORDS OF THE 
54630CPR CORNERS OF THE ELM 

54640CPR IT CAN GET THIS INFO FROM FILE OR FROM CORE 

54650CPR 

54660CPR 

54670CPR 

54680CPR CALLING ARGUMENTS — 

54690CPR NAME ATTRIBUTES DEFINITION 

54700CPR 

54710CPR IDONE (I) SEQUENTIAL POINTER TO HT ELM DATA 

54720CPR IWHR (I) INDICATOR 

54730CPR 1 -- GET ONLY WINDOW COORDS 

54740CPR 2 — GET CORNER COORDS 

54750CPR 

54760CPR COMMONS USED — 

54770CPR 

54780CPR AFIL 

54790CPR ELMDAT 

54800CPR SIZE 

54810CPR WORK 

54820CPR 

54830CPR FUNCTIONS/ROUTINES CALLED -- 
54840CPR 

54850CPR RWELMS 

54860CPR 

54870CPR FILES USED — 

54880CPR 

54890CPR 

54900CPR LIBRARIES ACCESSED -- 

54910CPR NONE 

54920CPR 

54930CPR 

54940CPR LOCAL VARIABLES — 

54950CPR 

54960CPR SPECIAL COMMENTS ABOUT THIS ROUTINE 
54970CPR 

54980CPR THE COORD DATA IS PASSED BACK TO THE ROUTINE -FNDELM- 
54990CPR VIA THE COMMON CALLED ELMDAT IN THE ARRAY AED 
55000CPR 
55010CPR 

55020CPR SPECIAL REMARKS/INSTRUCTIONS -- 

55030CPR HOST * H6000 

55040CPR 

55050CPR 

550€OCPft sttlKB *»BI«KSMII«a8S«SSBtl&CX»ZStSiaXSBSSISSBISCC 
55730 END 
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55750C= 


ssssssssssssssscsssssssssrsss 


55770 

55790C««* 

55810CPR 

55820CPR* 

55830CPR 

55840CPR 

55850CPR 

55860CPR 

55870CPR 

55880CPR 

55890CPR 

55900CPR 

55910CPR 

55920CPR 

55930CPR 

55940CPR 

55950CPR 

559S0CPR 

55970CPR 

55980CPR 

55990CPR 

5600OCPR 

56010CPR 

56020CPR 

56030CPR 

56050CPR 

5606 OCPR 

56070CPR 

56080CPR 

56090CPR 

56100CPR 

56110CPR 

56120CPR 

56130CPR 

56140CPR 

56150CPR 

56160CPR 

56170CPR 

56180CPR 

56190CPR 

56200CPR 

5621 OCPR 

56220CPR 

56230CPR 

56240CPR 

56250CPR 

56260CPR 

56270CPR 

56280CPR 

56290CPR 

56300CPR 

56310CPR 

56320CPR-" 

56610 


SUBROUTINE CKWIND( XP , YP . ZP, V , INOT ) 

SSS = = = = a3SSSSSSSSSSaK = BS = ZSSZSSSSB = SS 


:ci3iatscsi2c=:rssi:s;s33£ 


:sss==szsssss 


ROUTINE-- CKWIND 

DATE/WRITTEN BY— 04/07/83 RU MAFFEO 

DATE/REVISED BY— 

FUNCTION/PURPOSE — 

THIS ROUTINE CHECKS TO SEE IF THE COORDS OF THE 

STRESS POINT ARE WITHIN THE WINDOW OF THE HEAT TRANSFER 

ELEMENT 


CALLING ARGUMENTS — 

NAME ATTRIBUTES DEFINITION 

XP » 

YP »» (I) COORDS OF STRESS POINT 

ZP » 

COORDS DF HEAT TRAN ELM 
INOT (0) IN - OUT INDICATOR 

0 — OUT OF WINDOW 

1 -- IN WINDOW 


COMMONS USED -- 


FUNCTIONS/ROUTINES CALLED - 

FILES USED — 

LIBRARIES ACCESSED -- 
NONE 


LOCAL VARIABLES — 

SPECIAL COMMENTS ABOUT THIS ROUTINE 

SPECIAL REMARKS/INSTRUCTIONS — 

HOST « H6000 


==SSS=CSZXZSS=SKSZ«SB 


END 
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56620 SUBROUTINE MKRGF 

5€640C SSSSSSSS8BSSSalSSSSSS *SSSSSSSBSSf SSSSSSSBSSSBXBSSSSSStSS 

56650CPR 

56660CPR 

56670CPR ROUTINE-- MKRGF 
56680CPR 

5S690CPR DATE/WRITTEN BY-- 8/04/82 RJ MAFFEO 

56700CPR DATE/REVISED BY-- 4/04/83 RJ MAFFEO 

56710CPR REASON -- ADDED MIN MAX WINDOW LOGIC 

56720CPR DATE/REVISED BY— 8/30/83 SR MARTIN 

56730CPR REASON — CHANGED I/O TO READ NEW INPUT FILE FORM 

56740CPR 

56750CPR FUNCTION/PURPOSE — 

56760CPR THIS ROUTINE DOES THE PROCESSING OF THE NEUTRAL 
5G770CPR NODE AND ELEMENT FILES INTO THE INTERNAL DATA BASE FILES 

56780CPR FOR THIS VERSION OF 3D TITS. 

5S790CPR IT CALLS RWNODS TO PROCESS THE NODAL INFO. 

56800CPR IT CALLS RWELMS TO PROCESS THE ELM INFO 

56810CPR IT CALLS PROFAC TO PROCESS THE FACE INFORMATION 

56820CPR 
56830CPR 

56840CPR CALLING SEQUENCE— CALL MKRGF 
56850CPR 

56860CPR CALLING ARGUMENTS-- NONE 

56870CPR NAME ATTRIBUTES DEFINITION 

56880CPR 

56890CPR 

56900CPR FILES USED— 

5G910CPR IFL3 — RANDOM FILE TO STORE THE NODAL COORDINATES 

56920CPR I FL4 — RANDOM FILE TO STORE THE ELM ATTRIBUTES 

56930CPR INTITS — HEAT TRANSFER OUTPUT FILE 

56940CPR IN0T6 — OUTPUT FILE FOR MESSAGES 

56950CPR IERF — ERROR REMARK FILE 

56960CPR 

56970CPR 

56980CPR COMMONS USED — 

56990CPR AFIL -- FILECODE COMMON 

57000CPR ELMDAT -- ELM ATTRIBUTE COMMON 

57010CPR SIZE -- SIZE VARIABLES COMMON 

57020CPR CNTLFL — INFO WHICH DESCRIBES ALL ELEMENTS 

57030CPR 

57040CPR 

57050CPR 

57060CPR FUNCTIONS/ROUTINES CALLED— 

57070CPR RWNODS 

57080CPR HEADEL 

57090CPR RWELMS 

57 100CPR READPS 

57110CPR PROFAC 

57120CPR EMINMX 

57 130CPR 
57 140CPR 

57 150CPR LIBRARIES ACCESSED— 

57160CPR NONE 

57170CPR 
57 180CPR 

57 190CPR LOCAL VARIABLES— 

57200CPR NNODMX — MAXIMUM NODE NAME 

57210CPR NNODMN — MINIMUM NODE NAME 

57220CPR NUMNOD — NUMBER OF NODES 

57230CPR NELMX -- MAXIMUM ELM NAME 

57240CPR NELMN — MINIMUM ELM NAME 

57250CPR NUME -- NUMBER OF ELEMENTS 

57260CPR CRD(45) — ARRARY OF COORDINATES FOR THE 15 POINTS IN A ELM 

57270CPR IFACPC -- PACKED THTD FACE NUMBERS 

57280CPR IE(8) — ARRAY OF CONNECTIVITY FOR ELM CORENER POINTS 

57290CPR AED( 56 ) — ELM ATTRIBUTE ARRAY 
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TITLE= 


57300CPR 

57310CPR 

57320CPR 

57330CPR 

57340CPR 

57350CPR 

57360CPR 

57370CPR 

57380CPR 

57390CPR 

57400CPR 

57410CPR 

57420CPR 

57430CPR 

57440CPR 

57450CPR 

58690 


THIS ARRAY CONTAINS 
NE , CRD , IE, I FAC PC , VOL 
NE — ELM NAME 

VOL — VOLUME OF THE ELM 

XY2 ( 3 ) -- COORDINATE ARRAY FOR EACH NODE 


SPECIAL COMMENTS ABOUT THIS ROUTINE 
NOTE THE 15 POINTS FOR EACH ELM ARE 
8 CORNERS . 6 MIDSIDES AND 1 CENTROID 


SPECIAL REMARKS/ INSTRUCT I ONS““ 

HOST * H6000 

*********************************** 

END 
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58700 SUBROUTINE PROFAC 

587 20C = * = = = = = = * = = = = = = = a = = = * = = = = = = = 

58730CPR 

58740CPR 

58750CPR ROUTINE-- PROFAC 

58760CPR 

58770CPR DATE/WRITTEN BY-- 8/04/82 RJ MAFFEO 

58780CPR DATE/REVISED BY— 

58790CPR 

58800CPR FUNCTION/PURPOSE — 

58810CPR THIS ROUTINE READS THE PACKED THTD FACE NUMBER FROM THE 

58820CPR END OF THE NEUTRAL ELEMENT FILE AND STORES THEM ON 

58830CPR THE INTERNAL DATA BASE FILE 

58840CPR THIS DATA BASE FILE IS RANDOM FILE IFL4, 

58850CPR THE ELEMENTAL ATTRIBUTE FILE 

58860CPR NOTE THE PACKED FACE NUMBERS EXIST ON THE HEAT 

58870CPR TRANSFER OUTPUT FILE 

58880CPR IN RECORDS OF THE FORM- <ELMT *> <FACE INFO> 

58890CPR 

58900CPR 

58910CPR CALLING SEQUENCE— CALL PROFAC 

58920CPR 

58930CPR CALLING ARGUMENTS-- NONE 

58940CPR NAME ATTRIBUTES DEFINITION 

58950CPR 

58960CPR 

58970CPR FILES USED— 

58980CPR IFL4 -- RANDOM ELEMENT ATTRIBUTE FILE 

58990CPR INTITS — HEAT TRANSFER OUTPUT FILE (INPUT TO TITAN) 

59000CPR IERF -- ERROR REMARK FILE 

59010CPR 

59020CPR 

59030CPR COMMONS USED — 

59040CPR AFIL — FILECODE COMMON 

59050CPR ELMDAT — ELEMENT ATTRIBUTE COMMON 

59060CPR SIZE -- SIZE OF PROBLEM VARIABLES IN THIS COMMON 

59070CPR 

59080CPR 

59090CPR 

59100CPR FUNCTIONS/ROUTINES CALLED— 

59110CPR RWELMS 

59120CPR 

59130CPR 

59140CPR LIBRARIES ACCESSED— 

59150CPR NONE 

59160CPR 

59 170CPR LOCAL VARIABLES — 

59180CPR NRED NUMMBER OF NUMBERS READ FROM NEF FOR ONE SET OF 

59190CPR ELEMENT FACE DEFINITIONS (LIMIT IS 200) 

59200CPR IFACPC PACKED FACE NUMBER 

59210CPR NCNT — TOTAL NUMBER OF ELMS THAT HAVE FACES DEFINED 
59220CPR 

59230CPR SPECIAL COMMENTS ABOUT THIS ROUTINE 

59240CPR 

59250CPR 

59260CPR SPECIAL REMARKS/INSTRUCTIONS — 

59270CPR HOST 3 H6000 

59280CPR 

59290CPR ********************************* 

59890 END 
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59930 SUBROUTINE SKIP (IFC.ISETS) 

59960CPR 

59980CPR==================== ======= ===================3*===== 

59990CPR 

60000CPR 

6OO10CPR ROUTINE-- SKIP 

60020CPR 

60030CPR DATE/WRITTEN BY— 08/31/83 SR MARTIN 

60040CPR DATE/REVISED BY-- 

60050CPR 

60060CPR FUNCT I ON/PURPOSE -- 

60070CPR THIS ROUTINE ALLOWS THE CALLING ROUTINE TO SKIP PAST 

60080CPR UNNEEDED TIMESTEP AND TEMPERATURE DATA 

60090CPR 

60100CPR CALLING ARGUMENTS — 

60110CPR NAME ATTRIBUTES DEFINITION 

60120CPR 

60130CPR IFC INPUT, INTEGER FILE CODE 

60140CPR I SETS INPUT , INTEGER NUMBER OF SETS TO SKIP OVER 

60150CPR 

60160CPR COMMONS USED -- 
60170CPR CNTLFL 

60180CPR SIZE 

60190CPR 
60200CPR 

60210CPR FUNCTIONS/ROUTINES CALLED -- 

60220CPR 

60230CPR 

60240CPR FILES USED -- 
60250CPR 

60260CPR FC S 

60270CPR 

60280CPR LIBRARIES ACCESSED -- 

60290CPR NONE 

60300CPR 

60310CPR 

60320CPR LOCAL VARIABLES — 

60330CPR 

60340CPR DUMMY (CHARACTER) USED TO DUMMY READ LINES IN THE FILE 

60350CPR 

60360CPR SPECIAL COMMENTS ABOUT THIS ROUTINE 

60370CPR 

60380CPR 

60390CPR SPECIAL REMARKS/INSTRUCTIONS — 

60400CPR HOST = H6000 

60410CPR 

60420CPR 

604 30C P R ================ = ======== E ===== 3 = C = K *=*==■ = = X S 3 S = B = S 

60630 END 
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6 06 50C = = = = = : * = = = = = = **« = «*e = = *s* = * = = ** = = = = = = ss = = ece*ss = = = =sn* = 

60670 SUBROUTINE HTFIO ( I FC , DATA . ICNT . I STAT ) 

60690C a3=SSBBaacs3=E *= s **=«s=aa3*a==aaBE*8*SBBa3aBES=Bssf b 

607 1 OC BBS 

60720CPR 

60730CPR 

60740CPR ROUTINE-- HTFIO 

60750CPR 

60760CPR DATE/WRITTEN BY-- 09/13/83 SR MARTIN 

60770CPR DATE/REVISED BY — 

60780CPR 

60790CPR FUNCTION/PURPOSE — 

60800CPR THIS ROUTINE READS INFO FROM THE HEAT TRANSFER INPUT FILE 
60810CPR 

60820CPR CALLING ARGUMENTS — 

60830CPR I FC (INPUT/INTEGER) FILE CODE OF FILE TO READ FROM 

60840CPR DATA (OUTPUT/REAL) ARRAY OF DATA PASSED BACK 

60850CPR ICNT ( INPUT/ INTEGER ) AMOUNT OF DATA IN ' DATA ' ARRAY 
60860CPR ISTAT (OUTPUT/INTEGER) RETURN STATUS, VALUES ARE — 
60870CPR O-OK. 1-END OF FILE, 2-READ ERROR 

60880CPR 

60890CPR COMMONS USED -- 

60900CPR 

60910CPR 

60920CPR FUNCTIONS/ROUTINES CALLED -- 

60930CPR 

60940CPR 

60950CPR FILES USED -- 
60960CPR 

60970CPR I FC S 

60980CPR 

60990CPR LIBRARIES ACCESSED -- 
6 1000CPR NONE 
6 1010CPR 
6 1020CPR 

6 1030CPR LOCAL VARIABLES — 

61040CPR 

6 1050CPR SPECIAL COMMENTS ABOUT THIS ROUTINE 
6 1060CPR 
6 1070CPR 

6 1080CPR SPECIAL REMARKS/INSTRUCTIONS — 

6 1090CPR HOST - H6000 

61 100CPR 

61110CPR 

6 1 1 20CPR -* *««***««* * amummwuummmmsmmmmm muummm 
61270 END 
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61280 SUBROUTINE RDCNTL 

6 1 290CPR 

6 1 300CPR = = sss3s32s*a = stss = essas3ssssssasss*rssssss = t5ssss*ss = 

61310CPR 

61320CPR 

61330CPR ROUTINE-- RDCNTL 

61340CPR 

61350CPR DATE/WRITTEN BY-- 09/13/83 SR MARTIN 

6 1360CPR DATE/REVISED BY-- 

6 1370CPR 

61380CPR FUNCTION/PURPOSE-- 

61390CPR READ USER HEADER INFORMATION AND FILE COUNTERS 

6 1400CPR 

61410CPR 

61420CPR 

61430CPR 

6 1440CPR CALLING ARGUMENTS — 

6 1450CPR NAME ATTRIBUTES DEFINITION 

61460CPR 
6 1470CPR 

6 1480CPR COMMONS USED — 

6 1 490CPR 

61500CPR AFIL 

6 1 5 IOC PR CNTLFL 

6 1 520CPR CNTLTM 

61530CPR SIZE 

61540CPR 

6 1 550CPR FUNCTIONS/ROUTINES CALLED — 

6 1 560CPR 

6 1570CPR HTFIO 

6 1 580CPR 

61590CPR FILES USED -- 
6 1600CPR 

61610CPR INTITS S 

6 1620CPR IN0T6 S 

6 1630CPR 

6 1640CPR LIBRARIES ACCESSED — 

6 1650CPR NONE 
6 1660CPR 
6 1670CPR 

6 1680CPR LOCAL VARIABLES — 

6 1690CPR 

61700CPR SPECIAL COMMENTS ABOUT THIS ROUTINE 

61710CPR 

61720CPR 

6 1730CPR SPECIAL REMARKS/INSTRUCTIONS — 

6 1 740CPR HOST * H6000 

61750CPR 

61760CPR 

6 1770CPR eiB8S3Stie8a = SB **B>:BE 

62640 END 
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62650 SUBROUTINE TFLSET 

62660CPR 

62670CPR* 

62680CPR 

62690CPR 

62700CPR ROUTINE-- TFLSET 

62710CPR 

62720CPR DATE/WRITTEN BY-- 09/14/83 SR MARTIN 

62730CPR DATE/REVISED BY-- 

62740CPR 

62750CPR FUNCTION/PURPOSE-- 

62760CPR THIS ROUTINE REWINDS THE HEAT TRANSFER FILE AND REPOSITIONS 

62770CPR THE FILE POINTER TO THE FIRST RECORD OF THE TIME/TEMP 

62780CPR INFORMATION 

62790CPR 

62800CPR 

62810CPR 

62820CPR 

62830CPR CALLING ARGUMENTS — 

62840CPR NAME ATTRIBUTES DEFINITION 

62850CPR 

62860CPR 

62870CPR COMMONS USED -- 
62880CPR 

62890CPR AFIL 

62900CPR CNTLTM 

62910CPR 

62920CPR FUNCTIONS/ROUTINES CALLED -- 

62930CPR 

62940CPR 

62950CPR FILES USED — 

62960CPR 

62970CPR INTITS S 

62980CPR 

62990CPR LIBRARIES ACCESSED — 

63000CPR NONE 

63010CPR 

63020CPR 

63030CPR LOCAL VARIABLES — 

63040CPR 

63050CPR SPECIAL COMMENTS ABOUT THIS ROUTINE 

63060CPR 

63070CPR 

63080CPR SPECIAL REMARKS/INSTRUCTIONS — 

63090CPR HOST * H6000 

63 100CPR 
631 10CPR 

63120CPR = CSI ** X 3 B 3 | S | BS»IS»«S»»IB 

63420 END 
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63440C**** 

63460 

63480C==== 

63500CPR 

63510CPR»= 

63520CPR 

63530CPR 

63540CPR 

63550CPR 

63560CPR 

63570CPR 

63580CPR 

63590CPR 

63600CPR 

63610CPR 

63620CPR 

63630CPR 

63640CPR 

63650CPR 

63660CPR 

63670CPR 

63680CPR 

63690CPR 

63700CPR 

63710CPR 

63720CPR 

63730CPR 

63740CPR 

63750CPR 

63760CPR 

63770CPR 

63780CPR 

63790CPR 

63800CPR 

63810CPR 

63820CPR 

63830CPR 

63840CPR 

63850CPR 

63860CPR 

63870CPR 

63880CPR 

63890CPR 

63900CPR 

63910CPR 

63920CPR 

63930CPR 

63940CPR 

63950CPR 

63960CPR 

63970CPR 

63980CPR 

63990CPR 

64000CPR 

64010CPR 

64020CPR 

64030CPR 

64040CPR 

64050CPR 

64060CPR 

64070CPR 

64080C°R 

64090CPR 

64100CPR 

64 1 10CPR== : 

65280 


SUBROUTINE GENTMP 




ROUTINE-- GENTMP 

DATE/WRITTEN BY-- 09/26/83 RU MAFFEO 

DATE/REVISED BY-- 

FUNCTION/PURPOSE-- 

THIS ROUTINE IS THE MINI-EXEC USED TO CALL ALL 

ROUTINES ASSOCIATED WITH COMPUTING THE FINAL VALUES 

OF TEMPERATURE FOR EACH STRESS NODE 

IT CAN DEAL WITH BOTH FINITE DIFFERENCE AND 

FINITE ELEMNT HEAT TRANSFER CODES 

IF THE HEAT TRANSFER CODE IS FINITE DIFFERENCE 

THE TEMPERATURES OF THE CORNERS OF THE HEAT TRANSFER 

ELEMENTS ARE FOUND VIA CALLS FROM THIS RUOTINE 

CALLING ARGUMENTS — NONE 


COMMONS USED -- 

MAFCOM 

AFIL 

MAFCOM 

SIZE 

MAFCOM 

WORK 

MAFCOM 

CNTLTM 

CNTLFL 

FUNCTIONS/ROUTINES CALLED — 

GETMPE 

SETFLP 

RDTEMP 

CWCRNT 

ETCORE 

NTCORE 

STSTMP 

FILES USED -- 


LIBRARIES ACCESSED -- 
NONE 


LOCAL VARIABLES — 

SPECIAL COMMENTS ABOUT THIS ROUTINE 

SPECIAL REMARKS/INSTRUCTIONS — 

HOST « H6000 


B=CSSBS3SXSBSB==aSSSSSS=SS3SS=3SSS==&BSSSS=X3«SaS 

END 
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65300C BBaEBBBBBBBBBBBaBBBBBBBBBBBBBaBaBBBBBBBBBEBBBBB * BBBB 
65320 SUBROUTINE STSTMP(TEMP ,TIMVAL ) 

65340C BBBaBBBBBBBBBBBBBB = BBBBBaBBBBaaaBBaBBBBBBBaaBBaBaBBB 

65360CPR 

65370CPR BBBBaBBBBBBBBBBBBBBaBBBBBBBBBBBBBBBBBBaBBBaBBBBBBBBBB 

65380CPR 

65390CPR 

65400CPR ROUTINE— STSTMP 

65410CPR 

65420CPR DATE/WRITTEN BY— 09/28/83 Rd MAFFEO 

65430CPR DATE/REVISED BY — 

65440CPR 

65450CPR FUNCTION/PURPOSE — 

65460CPR 

65470CPR THIS ROUTINE READS THE COEFF INFO FOR EACH STRESS POINT 

65480CPR (VIA CALL RWSPCF) GETS THE CORNER TEMPS FROM THE TEMP ARRAY 

65490CPR AND COMPUTES THE TEMPERATURE OF THE STRESS POINT BY 
65500CPR SUMMING THE COEFFS TIMES THE CORNER TEMPS 
65510CPR 

65520CPR CALLING ARGUMENTS — 

65530CPR NAME ATTRIBUTES DEFINITION 

65540CPR 

65550CPR TEMP (I) ARRAY CONTAINING TEMPS OF HEAT 

65560CPR TRANSFER CORNERS 

65570CPR TIMVAL (I) TIME STEP VALUE 

65580CPR 

65590CPR COMMONS USED -- 
65600CPR 

65610CPR MAFCOM 

65620CPR AFIL 

65630CPR MAFCOM 

65640CPR SIZE 

65650CPR SPDAT 

65660CPR 

65670CPR FUNCTIONS/ROUTINES CALLED — 

65680CPR 

65690CPR RWSPCF 

65700CPR WRSTMP 

65710CPR 

65720CPR FILES USED — 

65730CPR 

65740CPR 

65750CPR LIBRARIES ACCESSED — 

65760CPR NONE 

65770CPR 

65780CPR 

65790CPR LOCAL VARIABLES — 

65800CPR 

65810CPR SPECIAL COMMENTS ABOUT THIS ROUTINE 

65820CPR 

65830CPR 

65840CPR SPECIAL REMARKS/INSTRUCTIONS — 

65850CPR HOST « H6000 

65860CPR 

65870CPR 

65880CPR BB S>>>B>l*B>SB>SBSBB*asBB8SSBSBBS>aSBS8SS8ECB8SSSCZZ 
66460 END 
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66480C ** = = = * s = *a = = = = = = : = = = = = = = = e = 3s = = s = s = = s3s = aa = = sea = s = = = = = 

66500 SUBROUTINE WRSTMP ( NAME , TEMP . TIMVAL , NDONE ) 

66520C * = * — = = = = = = -- = = = = = = = s«s = = = = = = = s = = s = a = ses = s = = s = a = a = = = = 

66540CPR 

66550CPR =3 * a = ss = 3S =*==*=====*s===e=s=S3s===ass==ris=ss====srs 

66560CPR 

66570CPR 

66580CPR ROUTINE-- WRSTMP 

66590CPR 

66600CPR DATE/WRITTEN BY — 09/28/83 RJ MAFFEO 

66610CPR DATE/REVISED BY-- 

66620CPR 

66630CPR FUNCTION/PURPOSE-- 
66640CPR 

66650CPR THIS ROUTINE WRITES THE STRESS POINT TEMPERATURES 
66660CPR TO THE OUTPUT FILE (IN NEUTRAL TEMPERATURE FILE FORM) 

66670CPR 

66680CPR CALLING ARGUMENTS — 

66690CPR NAME ATTRIBUTES DEFINITION 

66700CPR 

66710CPR NAME (I) STRESS POINT LABEL 

66720CPR TEMP (I) TEMPERATURE OF STRESS POINT 

66730CPR TIMVAL (I) TIME VALUE OF STEP 

66740CPR NDONE (I) H OF NODES WRITTEN 

66750CPR 

66760CPR COMMONS USED -- 
66770CPR 

66780CPR MAFCOM 

66790CPR AFIL 

66800CPR 

66810CPR FUNCTIONS/ROUTINES CALLED -- NONE 

66820CPR 

66830CPR 

66840CPR FILES USED -- 
66850CPR 

66860CPR IOUTF S 

66870CPR 

66880CPR LIBRARIES ACCESSED — 

66890CPR NONE 

66900CPR 

66910CPR 

66920CPR LOCAL VARIABLES — 

66930CPR 

66940CPR SPECIAL COMMENTS ABOUT THIS ROUTINE 

66950CPR 

66960CPR 

66970CPR SPECIAL REMARKS/INSTRUCTIONS — 

66980CPR HOST * H6000 

66990CPR 

67000CPR 

6701 OCPR* *«»»*S3S*«3»»383S*S*BES*SrSB*383«a«83rB3»«3e»Ba*SSS 

67280 END 
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67300C e3 * SS3 = BS= * s * B = aa *ase*=s=s*s#r=aas=a=s=BS3=e=s=*sss= 

67320 SUBROUTINE ETCORE ( ETMP ) 

67340C S3SS63*3*S*SSSa«BSSSaS2=82SSSSSSCE31BS3B3S3SS:*SeSSS 

67360CPR 

67370CPR = = = a:a * ac * = = 3 — s ssssesx»a*3*sssM*sr = aaa*aaes! = ss*s** = as*s 

67380CPR 

67390CPR 

67400CPR ROUTINE-- ETCORE 

67410CPR 

67420CPR DATE/WRITTEN BY-- 09/29/83 RJ MAFFEO 

67430CPR DATE/REVISED BY-- 

67440CPR 

67450CPR FUNCTION/PURPOSE — 

67460CPR 

67470CPR THIS ROUTINE TRANSFERS THE FINITE DIFFERENCE ELEMENTAL 
67480CPR CORNER TEMPERATURES FROM SCRATCH FILE TO CORE 
67490CPR 

67500CPR CALLING ARGUMENTS — 

67510CPR NAME ATTRIBUTES DEFINITION 

67520CPR 

67530CPR ETMP (0) ARRAY CONTAINING ELM CORNER TEMPS 

67540CPR 

67550CPR COMMONS USED — 

67560CPR 

67570CPR AFIL 

67580CPR SIZE 

67590CPR 

67600CPR FUNCTIONS/ROUTINES CALLED — 

67610CPR 

67620CPR FILREW 

67630CPR 

67640CPR FILES USED — 

67650CPR 

67660CPR ICRNTF S 

67670CPR IERF S 

67680CPR 

67690CPR LIBRARIES ACCESSED -- 

67700CPR NONE 

67710CPR 

67720CPR 

67730CPR LOCAL VARIABLES — 

67740CPR 

67750CPR SPECIAL COMMENTS ABOUT THIS ROUTINE 

67760CPR 

67770CPR 

67780CPR SPECIAL REMARKS/INSTRUCTIONS — 

67790CPR HOST - H6000 

67800CPR 

67810CPR 

67820CPR SBaSSSIBaSs:EIBaaBBBBB * B *C B =*«S 

68230 END 
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682500**== 

68270 

682900==== 

S8310CPR 

68320CPR= : 

68330CPR 

68340CPR 

683500PR 

68360CPR 

68370CPR 

68380CPR 

68390CPR 

68400CPR 

68410CPR 

68420CPR 

68430CPR 

68440CPR 

68450CPR 

68460CPR 

68470CPR 

68480CPR 

68490CPR 

68500CPR 

685100PR 

68520CPR 

68530CPR 

68540CPR 

68550CPR 

68560CPR 

68570CPR 

68580CPR 

68590CPR 

68600CPR 

68610CPR 

68620CPR 

68630CPR 

68640CPR 

68650CPR 

68660CPR 

68670CPR 

68680CPR 

686900PR 

68700CPR 

68710CPR 

68720CPR 

68730CPR 

68740CPR 

68750CPR 

68760CPR 

68770CPR 

68780CPR 

68790CPR 

68800CPR 

68810CPR 

68820CPR 

68830CPR 

68840CPR 

688500PR 

68860CPR 

68870CPR 

68880CPR 

68890CPR 

68900CPR 

68910CPR 

68920CPR* 

69200 


SUBROUTINE NTCORE(TEMP) 

exeeeeeeeeeeessseeeexbe: 


ssseassa 


ROUTINE— MAIN 

DATE/WRITTEN BY-- 09/19/83 SR MARTIN 

DATE/REVISED BY-- 

FUNCTION/PURPOSE — 

THIS ROUTINE ASSUMES THE FILE POINTER IS POINTING AT THE 
FIRST NODE NUMBER, TEMPERATURE PAIR OF A 
TIMESTEP IT READS IN EACH PAIR AND STORES THE 
TEMPERATURE (FROM DATA( 2 ) ) INTO THE 'TEMP' ARRAY 
IN POSITION <NODE NUMBER> ( I E--DATA ( 1 ) ) 

THE RESULT BEING THAT, FOR SAY NODE 13, THE TEMPERATURE 
WOULD BE IN RECORD 13 OF THE 'TEMP' ARRAY 
UPON RETURN, THE FILE POINTER WILL POINT 
TO THE NEXT TIMESTEP VALUE 


CALLING ARGUMENTS -- 

NAME ATTRIBUTES DEFINITION 

TEMP REAL, OUTPUT ARRAY WHICH WILL CONTAIN 

TEMPERATURE VALUES NOTE — 
TEMP(I) WILL CONTAIN THE 
TEMPERATURE VALUE FOR NODE 'I' 


COMMONS USED — 

AFIL 

SIZE 

FUNCTIONS/ROUTINES CALLED -- 
HTFIO 

FILES USED — 


LIBRARIES ACCESSED — 
NONE 


LOCAL VARIABLES — 

DATA (REAL , ARRAY ) USED TO RETRIEVE THE NODE NUMBER, 
TEMPERATURE PAIR FORM THE INPUT FILE 
NODNUM (INTEGER) USED TO STORE THE NODE NUMBER 
AS AS INTEGER AFTER RETRIEVING IT VIA THE 
'DATA' ARRAY. 

SPECIAL COMMENTS ABOUT THIS ROUTINE 


SPECIAL REMARKS/INSTRUCTIONS — 
HOST « H6000 


sBsxsBESBBBXSEBESEEssBSBBsssssESBssssEassssEBscss 

END 
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69220C =*5===“2=====*=c=E=rs=E==aaaBc==s3=r==3==z====a=a2r= 

69240 SUBROUTINE WRCRNT ( TBAR , NODE ) 

69260C = 2 - E = caa23E:: = ::Er3r = = l: * = * = = = = = Bn = = = EZ = E = = = = = a*zi = 2 = = 2 
69280CPR 

69290CPR Basas = ==s * =ssa=B s=s=**«=3s==sa=ss»*«easeHssBsaesrss*s 

69300CPR 

69310CPR 

69320CPR ROUTINE— WRCRNT 

69330CPR 

69340CPR DATE/WRITTEN BY-- 09/29/83 RJ MAFFEO 

69350CPR DATE/REVISED BY — 

69360CPR 

69370CPR FUNCTION/PURPOSE — 

69380CPR 

69390CPR THIS ROUTINE WRITES THE FINITE DIFFERENCE CORNER 
69400CPR TEMPS TO THE SCRATCH FILE. 

69410CPR 

69420CPR CALLING ARGUMENTS — 

69430CPR NAME ATTRIBUTES DEFINITION 

69440CPR 

69450CPR TBAR (I) VALUE OF TEMPERATURE 

69460CPR NODE (I) H T ELM CORNER NAME 

69470CPR 

69480CPR COMMONS USED -- 
69490CPR 

69500CPR AFIL 

69510CPR 

69520CPR FUNCTIONS/ROUTINES CALLED -- 

69530CPR 

69540CPR 

69550CPR FILES USED — 

69560CPR 

69570CPR ICRNTF S 

69580CPR 

69590CPR LIBRARIES ACCESSED -- 

69600CPR NONE 

69610CPR 

69620CPR 

69630CPR LOCAL VARIABLES -- 
69640CPR 

69650CPR SPECIAL COMMENTS ABOUT THIS ROUTINE 

69660CPR 

69670CPR 

69680CPR SPECIAL REMARKS/INSTRUCTIONS — 

69690CPR HOST « H60O0 

69700CPR 

69710CPR 

69720CPR** ar*ie3i««*xas»*ri3**«a«*r**!:iiciii = : *i!*Ei*te*j«*it 
69930 END 
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699 50C 8CBS3SaSBaS8SSB«8S3BaSSS8eS8S88S*tSSSS3SS8SSaaSSS8S8 

69970 SUBROUTINE HTFPAS( I FC , NFP , I FI LP ) 

6999QC 

70010CPR 


70020CPR-***-*-- *** 


s s a s = s a 


70030CPR 

70040CPR 

70050CPR 

70060CPR 

70070CPR 

70080CPR 

70090CPR 

70100CPR 

701 10CPR 

70120CPR 

70130CPR 

70140CPR 

70150CPR 

70160CPR 

70170CPR 

70180CPR 

70190CPR 

70200CPR 

70210CPR 

70220CPR 

70230CPR 

70240CPR 

70250CPR 

70260CPR 

70270CPR 

70280CPR 

70290CPR 

70300CPR 

70310CPR 

70320CPR 

70330CPR 

70340CPR 

70350CPR 

70360CPR 

70370CPR 

70380CPR 

70390CPR 

70400CPR 

70410CPR 

70420CPR 

70430CPR 

70440CPR 

70450CPR 

70460CPR 

70470CPR 

70480CPR 

70490CPR 

70500CPR- =* 

70840 


ROUTINE-- HTFPAS 

DATE/WRITTEN BY-- 09/29/83 RJ MAFFEO 

DATE/REVISED BY-- 

FUNCT I ON/ PURPOSE -- 

THIS ROUTINE SKIPS PAST PORTIONS OF THE 
HEAT TRANSFER INPUT FILE 

IT CURRENTLY ONLY CAN SKIP PAST THE TIME-TEMPERATURE 
PORT I ON ( I FI LP*6 ) 

CALLING ARGUMENTS — 


NAME 

ATTRIBUTES 

DEFINITION 


I FC 

(i) 

FILE 

CODE OF FILE TO 

BE SKIPED 

NFP 

(i) 

# OF 

PARTITIONS TO SKIP 

IFILP 

(i) 

PORTION OF FILE (DATA 

TYPE )T0 BE SKIPPED 


I FI LP- 1 THRU 5 NOT CURRENTLY USED 
IFILP-6 SKIP TIME-TEMPERATURE RECORDS 


COMMONS USED -- 
SIZE 

FUNCTIONS/ROUTINES CALLED -- 


FILES USED -- 
I FC S 

LIBRARIES ACCESSED -- 
NONE 


LOCAL VARIABLES — 

SPECIAL COMMENTS ABOUT THIS ROUTINE 


SPECIAL REMARKS/INSTRUCTIONS ~ 
HOST « H6000 


sitt3:ziEis:iifxzcsxz:i:cs2i:33E = :c3iz;EZzcB2:2c::£: 

END 
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7 08 60C = = = = = « = = = = ::isn: = i:£c s t = ;cc3ii:: = :c 2 ri: = iatc := i = xcr:s 

70880 SUBROUTINE SORT( SEEDS # LAX . NDUM , IKEY ) 

70900C* = 

70910CPR 

70920CPR = = = * B “** = = = = *« = = = = a*:**a* = «***a = = sts**** = = ss*es***ssca = 

70930CPR 

70940CPR 

70950CPR ROUTINE-- SORT 

70960CPR 

70970CPR DATE/WRITTEN BY — ???????'> GE SORT ROUTINE 

70980CPR DATE/REVISED BY — 

70990CPR 

7 1000CPR FUNCTION/PURPOSE — 

7 10 IOC PR 

7 1020CPR THIS ROUTINE SORTS 1 DIMENSIONAL ARRAYS 
7 1030CPR 

71040CPR CALLING ARGUMENTS — 


7 1050CPR 

NAME 

ATTRIBUTES 

DEFINITION 

7 1060CPR 
7 1070CPR 

SEEDS 

(I/O) 

ARRAY 

TO BE SORTED 

7 1080CPR 

LAX 

(I) 

NUMBER 

OF ITEMS IN ARRAY 

7 1090CPR 

NDUM 

DUMMY 

DUMMY 


71 100CPR 

IKEY 

DUMMY 

DUMMY 


71110CPR 





71 120CPR 

COMMONS USED 

i — 

NONE 



71130CPR 

7 1 140CPR FUNCTIONS/ROUTINES CALLED — NONE 
71 150CPR 

71160CPR FILES USED -- NONE 
7 1 170CPR 

7 1 180CPR LIBRARIES ACCESSED -- 
7 1 190CPR NONE 
7 1 200CPR 

71210CPR LOCAL VARIABLES — 

7 1220CPR 

7 1230CPR SPECIAL COMMENTS ABOUT THIS ROUTINE 
7 1240CPR 
7 1250CPR 

7 1260CPR SPECIAL REMARKS/INSTRUCTIONS -- 

7 1 270CPR HOST * H6000 

7 1 280CPR 
7 1290CPR 
7 1 300CPR 
71550 END 
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715700=== 
71590 
716100=** 
7 1630CPR 
7 1640CPR 
7 1650CPR 
7 1660CPR 
7 1 6700 PR 
7 1680CPR 
7 1690CPR 
7 1700CPR 
71710CPR 
7 1 720CPR 
7 1 730CPR 
71740CPR 
7 1 750CPR 
7 1 760CPR 
7 1770CPR 
7 1780CPR 
7 1790CPR 
7 1800CPR 
71810CPR 
7 1820CPR 
7 1830CPR 
7 1840CPR 
7 1850CPR 
7 1860CPR 
71870CPR 
7 1880CPR 
7 1890CPR 
7 1900CPR 
71910CPR 
7 1 920CPR 
7 1930CPR 
71940CPR 
7 1950CPR 
7 1960CPR 
7 1970CPR 
7 1 980CPR 
7 1 990CPR 
72000CPR 
72010CPR 
72020CPR 
72030CPR 
72040CPR 
72050CPR 
72060CPR 
72070CPR 
72080CPR 
72090CPR 
72 1000 PR 
721 10CPR 
72120CPR 
72130CPR 
72140CPR 
72150CPR* 
72700 


SUBROUTINE HTICON( IHTIC ) 


S==S3S3B=SS==X 


ROUTINE-- HTICON 

DATE/WRITTEN BY-- 01/14/84 RJ MAFFEO 

DATE/REVISED BY-- 

FUNCT ION/ PURPOSE -- 

THIS ROUTINE CONVERTS FROM THE HEAT TRANSFER OUTPUT FORM 
TO THE INTERNAL 3D TITAN HEAT TRANSFER INPUT FORM 
VALID HEAT TRANSFER INPUT CODES ARE 

0 INTERNAL 3D TITAN HEAT TRANSFER INPUT FORM 

1 SINDA HEAT TRANSFER INPUT FORM 

2 MARC HEAT TRANSFER INPUT FORM (FORMATTED) 

3 — MARC HEAT TRANSFER INPUT FORM (SEQUENTIAL BINARY) 

4 THTD HEAT TRANSFER INPUT FORM 


CALLING ARGUMENTS -- 

NAME ATTRIBUTES DEFINITION 

IHTIC (I) HEAT TRANSFER INPUT CODE 

COMMONS USED — 

AFIL 

FUNCTIONS/ROUTINES CALLED -- 

SINTIT 

MARCRD 

FILES USED — 

IERF S 

LIBRARIES ACCESSED -- 
NONE 


LOCAL VARIABLES -- 

SPECIAL COMMENTS ABOUT THIS ROUTINE 

SPECIAL REMARKS/INSTRUCTIONS -- 
HOST * H6000 


*CeCXB«SaatB&*SXSS=:EtBBIBS8ISIBS:SB3SSSBB;B3:8ttB 

END 
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72710 SUBROUTINE MARCRD ( I FTYP , NSTEPS ) 

72720CPR 

7 27 30CPR* * = = = = = = = 2 = 2* = E = I = C3r = = = 22 = Z = = = CS = == = = 2 = = a = -S = = = = r = = = 

72740CPR 

72750CPR 

72760CPR ROUTINE-- MARCRD 

72770CPR i 

72780CPR DATE/WRITTEN BY-- 01/11/84 RD MCCLAIN 

72790CPR DATE/REVISED BY-- 

72800CPR 

72810CPR FUNCTION/PURPOSE — 

72820CPR 

72830CPR READ MARC POST TAPE AND WRITE TRANSFER MODULE 

72840CPR INPUT FILE 

72850CPR 

72860CPR CALLING ARGUMENTS -- 

72870CPR NAME ATTRIBUTES DEFINITION 

72880CPR 

72890CPR I FTYP INPUT FILE TYPE OF MARC POST TAPE 

72900CPR 0*F0RMATTED , 1*SEQ BINARY 

72910CPR , NSTEPS INPUT NUMBER OF STEPS TO READ FROM 

72920CPR MARC POST TAPE 

72930CPR 

72940CPR COMMONS USED -- 
72950CPR AFIL -- FILE CODES 

72960CPR 

72970CPR FUNCTIONS/ROUTINES CALLED — 

72980CPR 

72990CPR READR READS DATA FROM POST TAPE 

73000CPR SKPREC SKIPS RECORDS ON POST TAPE 

73010CPR 

73020CPR FILES USED — 

73030CPR 

73040CPR INTITS S TRANSFER MODULE INPUT FILE 

73050CPR IHTIN $ MARC INPUT FILE (SEQ BIN OR FORMATTED) 

73060CPR 

73070CPR LIBRARIES ACCESSED — 

73080CPR NONE 
73090CPR 
73 100CPR 

73110CPR LOCAL VARIABLES — 

73120CPR 

73130CPR SPECIAL COMMENTS ABOUT THIS ROUTINE 

73 140CPR 

73150CPR 

73 160CPR SPECIAL REMARKS/INSTRUCTIONS — 

73 170CPR HOST * H6000 

73 180CPR 
73 190CPR 

73200CPR 3 >3 >>3ii3ei:i:i:ser:e338iE«siK8Bi«iEts33&i=Etiiiii3t 
74770 END 
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74810 

74840CPR 

74850CPR 

74860CPR 

74870CPR 

74880CPR 

74890CPR 

74900CPR 

74910CPR 

74920CPR 

74930CPR 

74940CPR 

74950CPR 

74960CPR 

74970CPR 

74980CPR 


SUBROUTINE READR ( LFN, DATA , NDATA . CHAR .NWRDS , TTYPE, IFTYP , I END) 


sasssxssssss 


ROUTINE-- 


S2=SSSSS«=SB=a 


75160CPR 
75170CPR 
75 180CPR 
75190CPR 
75200CPR 
75210CPR 
75220CPR 
75230CPR 
75240CPR 


READR 


RD MCCLAIN 


DATE/WRITTEN BY— 01/11/84 
OATE/REVISED BY — 

FUNCT I ON/PURPOS E — 


READ FLOATING POINT , INTEGER . OR CHARACTER DATA 
FROM MARC POST TAPE 


CALLING ARGUMENTS — 


74990CPR 

75000CPR 

NAME 

ATTRIBUTES 

75010CPR 

LFN 

INPUT 

75020CPR 

DATA 

OUTPUT 

75030CPR 

NDATA 

OUTPUT 

75040CPR 

CHAR 

OUTPUT 

75050CPR 

NWRDS 

INPUT 

75060CPR 

75070CPR 

75080CPR 

75090CPR 

ITYPE 

INPUT 

75100CPR 

75110CPR 

75120CPR 

75130CPR 

IFTYP 

INPUT 

75140CPR 

COMMONS USED 

-- 

75150CPR 

AFIL — 

FILE CODES 


DEFINITION 

FILE CODE TO READ FROM 
FLOATING POINT ARRAY 
INTEGER ARRAY 
CHARACTER ARRAY 
NUMBER OF ITEMS TO READ 
TYPE OF DATA TO READ 

1- LOAD FP DATA INTO ARRAY DATA 

2- LOAD INTEGER DATA INTO ARRAY NDATA 

3- LOAD CHAR DATA INTO ARRAY CHAR 
FILE TYPE 

0- FORMATTED 

1 - SEQUENTIAL BINARY 


FUNCTIONS/ROUTINES CALLED -- 


FILES USED — 

LFN S 

IN0T6 S 


FILE CODE TO BE READ 
STANDARD OUTPUT FILE 


75250CPR 

75260CPR 

75270CPR 

75280CPR 

75290CPR 

75300CPR 

75310CPR 

75320CPR 

75330CPR 

75340CPR 

75350CPR 

75360CPR 

75370CPR 

75380CPR- 

75870 


LIBRARIES ACCESSED — 
NONE 


LOCAL VARIABLES — 

SPECIAL COMMENTS ABOUT THIS ROUTINE 


SPECIAL REMARKS/INSTRUCTIONS — 
HOST « H6000 


ENO 
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75900 SUBROUTINE SKPREC( IFI LE .NUMREC , IFTYP ) 

75930CPR 

75940CPR- - i:i==::===i::s:triBaE======r=;3 = i == a ::== 3- ===rtt - 3= 

75950CPR 

75960CPR 

75970CPR ROUTINE-- SKPREC 

75980CPR 

75990CPR DATE/WRITTEN BY-- 01/11/84 RD MCCLAIN 

76000CPR DATE/REVISED BY — 

76010CPR 

76020CPR FUNCTION/PURPOSE-- 
76030CPR 

7S040CPR SKIP RECORDS ON A SEQUENTIAL FILE 

76050CPR 

76060CPR 

76070CPR CALLING ARGUMENTS — 

76080CPR NAME ATTRIBUTES DEFINITION 

76090CPR 

76100CPR IFILE INPUT FILE CODE TO SKIP RECORDS ON 

76110CPR NUMREC INPUT NUMBER OF RECORDS TO SKIP 

76 120CPR IFTYP INPUT FILE TYPE OF MARC POST TAPE 

76130CPR O'FORMATTED, 1 “BINARY 

76140CPR 

76150CPR COMMONS USED — NONE 

76160CPR 

76170CPR 

76 180CPR FUNCTIONS/ROUTINES CALLED — 

76 190CPR 
76200CPR 

76210CPR FILES USED — 

76220CPR 

76230CPR 

76240CPR LIBRARIES ACCESSED — 

76250CPR NONE 

76260CPR 

76270CPR 

76280CPR LOCAL VARIABLES — 

76290CPR 

76300CPR SPECIAL COMMENTS ABOUT THIS ROUTINE 

76310CPR 

76320CPR 

76330CPR SPECIAL REMARKS/INSTRUCTIONS — 

76340CPR HOST * H6000 

76350CPR 

76360CPR 

76370CPR* «***»«*»*««*»*««b»*«»**o»= 

76690 END 
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767 10C 33 

76730 

76750C* 3 

76770CPR 

76780CPR 

76790CPR 

76800CPR 

76810CPR 

76820CPR 

76830CPR 

76840CPR 

76850CPR 

76860CPR 

76870CPR 

76880CPR 

76890CPR 

76900CPR 

76910CPR 

76920CPR 

76930CPR 

76940CPR 

76950CPR 

76960CPR 

76970CPR 

76980CPR 

76990CPR 

77000CPR 

77010CPR 

77020CPR 

77030CPR 

77040CPR 

77050CPR 

77060CPR 

77070CPR 

77080CPR 

77090CPR 

77100CPR 

77110CPR 

77120CPR 

77130CPR 

77 140CPR 

77150CPR 

77 160CPR 

77170CPR 

77 180CPR 

77 190CPR 

77200CPR 

77210CPR 

77220CPR 

77230CPR 

77240CPR 

77250CPR 

77260CPR 

77270CPR 

77280CPR 

77290CPR 

77300CPR 

77320CPR 

77340CPR 

77350CPR 

77360CPR 

77370CPR 

77380CPR 

77390CPR 

7740OCPR 

77410CPR 

77420CPR 


SUBROUTINE SINTIT 




ROUTINE— SINTIT 

DATE/WRITTEN BY— 01/09/84 RN PITTMAN 

DATE/REVISED BY- 

FUNCTION/ PURPOSE — 

MAIN EXECUTIVE ROUTINE FOR CREATION OF SINDA-T0-3DTITAN 
INTERFACE FILE 

CALLING ARGUMENTS — 

NAME ATTRIBUTES DEFINITION 

NONE 

COMMONS USED — 

WORK - STORES INTERNAL ARRAY DATA 


TEMP ARRAY OF ELEMENTAL TEMPERATURES 

CORRESPONDING TO ' ITELM ' 

IPKFAC ARRAY OF PACKED FACE NUMBERS 

CORRESPONDING TO ' IGELM' 

IGELM ARRAY OF ELEMENTS WITH GEOMETRY 

ITELM ARRAY OF ELEMENTS WITH TEMPERATURES 

AFIL - FILE CODES 


ELEMENTS DEFINITION 


IHTIN SINDA TEMPERATURE OUTPUT FILE (INPUT) 

ISINGF SINDA GEOMETRY FILE (INPUT) 

INTITS INTERFACE FILE (OUTPUT) 

I ERF HARDCOPY FILE (OUTPUT) 

KNTDAT - COUNTER DATA 


ELEMENTS DEFINITION 


NUMNOD NUMBER OF NODES (CORNER POINTS) 

IN MODEL 

NUMEWG NUMBER OF ELEMENTS WITH GEOMETRY 

NUMTIM NUMBER OF SOLUTION TIMES 

NNPE NUMBER OF NODES PER ELEMENT 

ECODE ELEMENT CODE 

*1 * 2D MODEL ( NNPE*4 ) 

*11 -3D MODEL ( NNPE*8 ) 

ACODE ANALYSIS CODE 

=1 “ FINITE DIFFERENCE ANALYSIS 

NUMEWT NUMBER OF ELEMENTS WITH TEMPERATURES 

NUMSKP NUMBER OF RECORDS TO SKIP ON 'IHTIN ' 

PRIOR TO BEGINNING OF TEMPERATURE DATA 
NUMFUL NUMBER OF FULL RECORDS OF DATA 

NUMREM NUMBER OF WORDS ON FINAL RECORD 

FUNCTIONS/ROUTINES CALLED — 

DATINT 

PRCSOF 

PRCGEO 

WRTNTF 

FILES USED — 

NONE 
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77430CPR 

77440CPR LIBRARIES ACCESSED — 

77450CPR NONE 

77460CPR 

77470CPR LOCAL VARIABLES — 

77480CPR NONE 

77490CPR 

77500CPR SPECIAL COMMENTS ABOUT THIS ROUTINE ~ 

77510CPR PROGRAM IS CURRENTLY DIMENSIONED TO PERMIT A MAXIMUM OF 

77520CPR 5000 ELEMENTS WITH GEOMETRY AND 5000 ELEMENTS WITH 

77530CPR TEMPERATURES 

77540CPR 

77550CPR SPECIAL REMARKS/INSTRUCTIONS — 

77560CPR HOST * HIS 

77570CPR 
77580CPR 
7 7590C 

78060 END 
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78100 SUBROUTINE DATINT 

78140CPR 

781 50CPR = S = *XSSSSSSES53 = S3 = = 3S = = = SS = SS = = 32S = = SSB = = XX = SS= = S = = = 

78160CPR 

78170CPR 

78 180CPR ROUTINE-- DATINT 

7S190CPR 

78200CPR DATE/WRITTEN BY-- 01/11/84 RN PITTMAN 

78210CPR DATE/REVISED BY-- 

78220CPR 

78230CPR FUNCTI ON/PURPOSE- - 

78240CPR THIS ROUTINE INITIALIZES ALL OF THE GLOBAL COMMON VARIABLES 

78250CPR 

78260CPR CALLING ARGUMENTS — 

78270CPR NAME ATTRIBUTES DEFINITION 

78280CPR NONE 

78290CPR 

78300CPR COMMONS USED — 

78310CPR WORK -- STORES INTERNAL ARRAY DATA 

78320CPR AFIL — FILE CODES (MAIN) 

78330CPR KNTDAT -- COUNTER DATA (MAIN) 

78340CPR 

78350CPR FUNCTIONS/ROUTINES CALLED — 

78360CPR NONE 

78370CPR 

78380CPR FILES USED — 

78390CPR NONE 

78400CPR 

78410CPR LIBRARIES ACCESSED ~ 

78420CPR NONE 

78430CPR 

78440CPR LOCAL VARIABLES — 

78450CPR NONE 

784G0CPR 

78470CPR SPECIAL COMMENTS ABOUT THIS ROUTINE — 

78480CPR NONE 

78490CPR 

78500CPR SPECIAL REMARKS/INSTRUCTIONS — 

78510CPR HOST = H6000 

78520CPR 

78530CPR 

78540CPR “ *>Bzir3ni:sancs3i = XKCiiacst«assrE = s:asi:3cii2«i;E = = acasE 
79020C . END OF SUBROUTINE 'DATINT ' 

79050 END 
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79090 SUBROUTINE PRCSOF 

79130CPR 

791 40CPR==== **=SS*eSS338B3SS3S3SBSS3SSS3S333B3SS83S3383S3SS383 

79150CPR 

79160CPR 

79170CPR ROUTINE-- PRCSOF 

79180CPR 

79190CPR DATE/WRITTEN BY-- 01/11/84 RN PITTMAN 

79200CPR DATE/REVISED BY — 

79210CPR 

79220CPR FUNCTION/PURPOSE-- 

79230CPR THIS ROUTINE CONTROLS THE INITIAL PROCESSING OF THE 

79240CPR 'SINDA OUTPUT FILE' (THE FILE WHICH CONTAINS THE 

79250CPR RESULTS OF THE SINDA ANALYSIS) 

79260CPR 

7 9270CPR CALLING ARGUMENTS — 

79280CPR NAME ATTRIBUTES DEFINITION 

79290CPR NONE 

79300CPR 

79310CPR COMMONS USED -- 

79320CPR AFIL — FILE CODES (MAIN) 

79330CPR 

79340CPR FUNCTIONS/ROUTINES CALLED — 

79350CPR GETEWT 

79360CPR GTNUMT 

79370CPR 

79380CPR FILES USED — 

79390CPR NONE 

79400CPR 

79410CPR LIBRARIES ACCESSED — 

79420CPR NONE 

79430CPR 

79440CPR LOCAL VARIABLES -- 
79450CPR NONE 

79460CPR 

79470CPR SPECIAL COMMENTS ABOUT THIS ROUTINE -- 

79480CPR NONE 

79490CPR 

79500CPR SPECIAL REMARKS/INSTRUCTIONS — 

79510CPR HOST » H6000 

79520CPR 

79530CPR 

79540CPR" 

79770C . END OF SUBROUTINE 'PRCSOF' 

79800 END 
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79840 SUBROUTINE GETEWT 

79880CPR 

79890CPR===*==— =*==*====«— =*======*=*===============*=====- 

79900CPR 

79910CPR 

79920CPR ROUTINE— GETEWT 

79930CPR 

79940CPR DATE/WRITTEN BY— 01/11/84 RN PITTMAN 

79950CPR DATE/REVISED BY — 

79960CPR 

79970CPR FUNCTION/PURPOSE — 

79980CPR 1 READS THE SINDA ELEMENT NUMBERS FROM 'IHTIN' INTO 'ITELM' 

79990CPR 2 READS FROM ' IHTIN' THE TOTAL NUMBER OF SINDA ELEMENTS 

80000CPR WITH TEMPERATURES INTO ' NUMEWT ' 

80O10CPR 3 COUNTS THE NUMBER OF RECORDS PRECEEDING THE FIRST TEMP- 

80020CPR ERATURE SOLUTION ON 'IHTIN' AND STORES THIS VALUE IN 

80030CPR ' NUMSKP' . 

8004 OC PR 

80050CPR CALLING ARGUMENTS — 

80060CPR NAME ATTRIBUTES DEFINITION 

80070CPR NONE 

80080CPR 

80090CPR COMMONS USED — 

80100CPR WORK —STORES INTERNAL ARRAY DATA 

80110CPR AFIL — FILE CODES (MAIN) 

80120CPR KNTDAT — COUNTER DATA (MAIN) 

80130CPR 

80140CPR FUNCTIONS/ROUTINES CALLED — 

80150CPR ERRPRT 

80160CPR 

80170CPR FILES USED — 

80180CPR IHTIN (S) SINDA OUTPUT FILE (INPUT TO THIS ROUTINE) 

80190CPR 

80200CPR LIBRARIES ACCESSED — 

80210CPR NONE 

80220CPR 

80230CPR LOCAL VARIABLES — 

80240CPR IPOS - POSITION OF THE BEGINNING OF EACH GROUP OF 10 

80250CPR ELEMENT NUMBERS IN THE 'ITELM' ARRAY 

80260CPR ***NOTE*** 

80270CPR IN THIS ROUTINE, THE GLOBAL VARIABLES 'NUMFUL' AND NUMREM' 

80280CPR ARE REDEFINED AS FOLLOWS- 

80290CPR NUMFUL - NUMBER OF 'FULL' 10-WORD RECORDS OF SINDA 

80300CPR ELEMENT NUMBERS CONTAINED ON 'IHTIN' 

80310CPR NUMREM - NUMBER OF ELEMENT NUMBERS REMAINING ON FINAL 

80320CPR RECORD OF ELEMENT NUMBERS 

80330CPR 

80340CPR SPECIAL COMMENTS ABOUT THIS ROUTINE -- 

B0350CPR NONE 

80360CPR 

80370CPR SPECIAL REMARKS/INSTRUCTIONS — 

80380CPR HOST * H60O0 

80390CPR 

80400CPR 

804 1 OCPR* «*»**at*** = *«»*:scc«c«!*3* = n = cB = = iij:it = 23t = 3»tnan*a3zi: = = 

8121 OC . END OF SUBROUTINE 'GETEWT' 

81240 END 
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81280 
8 1 320CPR 
8 1 330CPR 
8 1340CPR 
8 1350CPR 
8 1 360CPR 
8 1 370CPR 
8 1380CPR 
8 1390CPR 
8 1400CPR 
8 1410CPR 
8 1420CPR 
8 1430CPR 
8 1440CPR 
8 1450CPR 
81460CPR 
81470CPR 
8 1480CPR 
8 1490CPR 
81500CPR 
81510CPR 
8 1520CPR 
8 1 530CPR 
8 1540CPR 
8 1550CPR 
8 1 560CPR 
81570CPR 
8 1 580CPR 
8 1590CPR 
8 1600CPR 
81610CPR 
8 1620CPR 
8 1630CPR 
81640CPR 
8 1650CPR 
8 1660CPR 
8 1670CPR 
81680CPR 
8 1690CPR 
8 1700CPR 
8 1710CPR 
8 1720CPR 
8 1730CPR 
8 1740CPR 
8 1750CPR 
81760CPR 
81770CPR 
8 1780CPR 
8 1790CPR 
8 1800CPR 
81810CPR 
8 1820CPR 
8 1830CPR 
8 1840CPR 
8 1850CPR 
8 1860CPR 
8 1870CPR 
8 1880CPR 1 
82550C . 
82580 


SUBROUTINE GTNUMT 


esssssasssa 


ROUTINE— GTNUMT 

DATE/WRITTEN BY— 01/11/84 RN PITTMAN 

DATE/REVISED BY- 

FUNCTION/ PURPOSE — 

THIS ROUTINE COUNTS THE TOTAL NUMBER OF SOLUTION BLOCKS 
CONTAINED IN THE 'SINDA OUTPUT FILE' AND STORES THIS VALUE 
IN THE GLOBAL VARIABLE 'NUMTIM' 

CALLING ARGUMENTS — 

NAME ATTRIBUTES DEFINITION 

NONE 


COMMONS USED — 

WORK —STORES INTERNAL ARRAY DATA 
AFIL — FILE CODES (MAIN) 

KNTDAT — COUNTER DATA (MAIN) 

FUNCTIONS/ROUTINES CALLED — 

ERRPRT 

FILES USED — 

IHTIN (S) SINDA OUTPUT FILE (INPUT TO THIS PROGRAM) 

LIBRARIES ACCESSED — 

NONE 

LOCAL VARIABLES — 

DUMMY - VARIABLE USED IN DUMMY READS TO PASS THROUGH EACH 
SOLUTION BLOCK TO THE NEXT 

KNT - CURRENT COUNT OF THE NUMBER OF SOLUTION BLOCKS 
ENCOUNTERED ON ' IHTIN' 

***NOTE*** 

IN THIS ROUTINE, THE GLOBAL VARIABLES ' NUMFUL ' AND 'NUMREM' 
ARE REDEFINED AS FOLLOWS- 

NUMFUL - NUMBER OF ' FULL ' 6-WORD RECORDS OF TEMPERATURE 
DATA IN EACH SOLUTION BLOCK ON 'IHTIN' 

NUMREM - NUMBER OF TEMPERATURES ON FINAL RECORD OF EACH 
SOLUTION BLOCK 

THESE DEFINITIONS ARE NOT ALTERED FURTHER IN THIS PROGRAM 

SPECIAL COMMENTS ABOUT THIS ROUTINE — 

EACH TEMPERATURE SOLUTION BLOCK ON 'IHTIN' IS PASSED THROUGH 
WITH DUMMY READS UNTIL THE END OF FILE IS ENCOUNTERED 
THE SOLUTION COUNTER, 'KNT', IS INCREMENTED BY ONE AT THE 
BEGINNING OF EACH BLOCK 

SPECIAL REMARKS/INSTRUCTIONS — 

HOST - H6000 


END OF SUBROUTINE 'GTNUMT' 
END 
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82620 SUBROUTINE PRCGEO 

82660CPR 

8 2 670C P R ==**»*===== sbss:3:::s:bsss=scs:css a s a s * = s SSS= s = 3 s === = 


82680CPR 

82690CPR 

82700CPR 

82710CPR 

82720CPR 

82730CPR 

82740CPR 

82750CPR 

82760CPR 

82770CPR 

82780CPR 

82790CPR 

82800CPR 

82810CPR 

82820CPR 

82830CPR 

82840CPR 

82850CPR 

82860CPR 

82870CPR 

82880CPR 

82890CPR 

82900CPR 

82910CPR 

82920CPR 

82930CPR 

82940CPR 

82950CPR 

82960CPR 

82970CPR 

82980CPR 

82990CPR 

83000CPR 

83010CPR 

83020CPR 

83030CPR 

83040CPR 

83050CPR 

83060CPR 

83070CPR 

83080CPR 

83090CPR** 

83560C . . . 

83590 


ROUTINE-- PRCGEO 

DATE/WRITTEN BY-- 01/11/84 RN PITTMAN 

DATE/REVISED BY — 

FUNCTION/PURPOSE — 

1 READS THE VALUES CONTAINED IN THE 'DIRECTORY' RECORD OF 
THE 'SINDA GEOMETRY FILE' INTO THE GLOBAL COMMON VALUES 
'NUMNOD', 'NUMEWG' , AND 'NNPE ' 

2 SETS THE VALUE FOR 'ECODE' BASED UPON THE VALUE OF 'NNPE' 

CALLING ARGUMENTS — 

NAME ATTRIBUTES DEFINITION 

NONE 


COMMONS USED -- 

AFIL — FILE CODES (MAIN) 

KNTDAT — CCOUNTER DATA (MAIN) 

FUNCTIONS/ROUTINES CALLED -- 
ERRPRT 

FILES USED — 

ISINGF (S) SINDA GEOMETRY FILE (INPUT) 

LIBRARIES ACCESSED — 

NONE 

LOCAL VARIABLES — 

TITLE - CHARACTER ARRAY INTO WHICH THE 'ISINGF' TITLE RECORD 
IS READ 

SPECIAL COMMENTS ABOUT THIS ROUTINE -- 
NONE 

SPECIAL REMARKS/INSTRUCTIONS — 

HOST ■ H6000 


*tnES*«nci«i*»iaiaci*5**tiiiaioii3i* = 32 *u:;*ritS 2 E 

END OF SUBROUTINE 'PRCGEO' 

END 
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83630 SUBROUTINE WRTNTF 

83670CPR 

83680CPR BSSSSSE5XSSKS& SSSSSS3SEISSSSS3S&SSZ3SS8XSBSSaS8=SSSI= 

83690CPR 

83700CPR 

83710CPR ROUTINE-- WRTNTF 

83720CPR 

83730CPR DATE/WRITTEN BY— 01/11/84 RN PITTMAN 

83740CPR DATE/REVISED BY — 

83750CPR 

83760CPR FUNCTION/PURPOSE — 

83770CPR THIS ROUTINE CONTROLS THE WRITING OF THE 'INTERFACE FILE'. 

83780CPR 

83790CPR CALLING ARGUMENTS — 

83800CPR NAME ATTRIBUTES DEFINITION 

83810CPR NONE 

83820CPR 

83830CPR COMMONS USED -- 

83840CPR WORKS —STORES INTERNAL ARRAY DATA 

83850CPR AFIL — FILE CODES (MAIN) 

83860CPR , KNTDAT — COUNTER DATA (MAIN) 

83870CPR 

83880CPR FUNCTIONS/ROUTINES CALLED — 

83890CPR WRTDIR 

83900CPR WRTGEO 

83910CPR WRTEMP 

83920CPR ERRPRT 

83930CPR 

83940CPR FILES USED — 

83950CPR INTITS (S) INTERFACE FILE (OUTPUT) 

83960CPR ISINGF (S) SINDA GEOMETRY FILE (INPUT) 

83970CPR IHTIN (S) SINDA OUTPUT FILE (INPUT) 

83980CPR 

83990CPR LIBRARIES ACCESSED — 

84000CPR NONE 

84010CPR 

84020CPR LOCAL VARIABLES — 

84030CPR DUMMY - VARIABLE USED IN 'DUMMY' READS TO SKIP OVER ALL 

84040CPR DATA ON 'IHTIN' PRECEEDING THE TEMPERATURE 

84050CPR SOLUTION BLOCKS 

84060CPR 

84070CPR SPECIAL COMMENTS ABOUT THIS ROUTINE -- 

84080CPR NONE 

84090CPR 

84 100CPR SPECIAL REMARKS/INSTRUCTIONS -- 

84110CPR HOST * H6000 

84120CPR 

84130CPR 

84140CPR> : »»<» 

84850C . END OF SUBROUTINE 'WRTNTF' 

84830 END 
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84920 SUBROUTINE WRTDIR 

849S0CPR 

84970CPR* = s = sss = sss3sssc = = sssssssssssaessssssssssss = B = = sesssa = 

84980CPR 

84990CPR 

85000CPR ROUTINE— WRTDIR 

85010CPR 

85020CPR DATE/WRITTEN BY— 01/11/84 RN PITTMAN 

85030CPR DATE/REVISED BY — 

85040CPR 

85050CPR FUNCTION/PURPOSE— 

85060CPR THIS ROUTINE WRITES THE TITLE (RECORDS 1-4) AND DIRECTORY 

85070CPR (RECORD 5) RECORDS OF THE 'INTERFACE FILE' 

85080CPR 

85090CPR CALLING ARGUMENTS — 

85100CPR NAME ATTRIBUTES DEFINITION 

85 1 10CPR NONE 

85120CPR 

85130CPR COMMONS USED -- 

85140CPR WORK —STORES INTERNAL ARRAY DATA 

85150CPR AFIL — FILE CODES (MAIN) 

85160CPR KNTDAT — COUNTER DATA (MAIN) 

85170CPR 

85180CPR FUNCTIONS/ROUTINES CALLED — 

85190CPR ERRPRT 

85200CPR 

85210CPR FILES USED — 

85220CPR INTITS (S) INTERFACE FILE (OUTPUT) 

85230CPR ISINGF (S) SINDA GEOMETRY FILE (INPUT) 

85240CPR 

85250CPR LIBRARIES ACCESSED — 

85260CPR NONE 

85270CPR 

85280CPR LOCAL VARIABLES -- 

85290CPR TITLE - 80-CHARACTER TITLE CONTAINED ON RECORD 1 OF 'ISINGF' 

85300CPR AND WRITTEN ON RECORD 4 OF 'INTITS' 

85310CPR USRTTL - 80-CHARACTER ARRAY USED TO 'DUMMY-OUT' RECORDS 

85320CPR 1-3 OF 'INTITS' 

85330CPR 

85340CPR SPECIAL COMMENTS ABOUT THIS ROUTINE — 

85350CPR NONE 

85360CPR 

85370CPR SPECIAL REMARKS/INSTRUCTIONS -- 

85380CPR HOST » H6000 

85390CPR 

85400CPR 

854 10CPR° >S 

86090C . END OF SUBROUTINE 'WRTDIR' 

86120 END 
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86160 SUBROUTINE WRTGEO 

86200CPR 

86 2 1 OC PR5SS3s:sss:ssess3setS3SSSs: = s3rccKsss::zss = : = r = B;; = :3t 

86220CPR 

86230CPR 

86240CPR ROUTINE-- WRTGEO 

86250CPR 

86260CPR DATE/WRITTEN BY— 01/11/84 RN PITTMAN 

86270CPR DATE/REVISED BY — 

86280CPR 

86290CPR FUNCTION/PURPOSE— 

86300CPR THIS ROUTINE WRITES THE GEOMETRICAL DATA (COORDINATE, 

86310CPR CONNECTIVITY, AND FACES) ON THE INTERFACE FILE 

86320CPR 

86330CPR CALLING ARGUMENTS — 

86340CPR NAME ATTRIBUTES DEFINITION 

86350CPR NONE 

86360CPR 

86370CPR COMMONS USED — 

86380CPR WORK —STORES INTERNAL ARRAY DATA 

86390CPR AFIL — FILE CODES (MAIN) 

86400CPR KNTDAT — COUNTER DATA (MAIN) 

86410CPR 

86420CPR FUNCTIONS/ROUTINES CALLED — 

86430CPR ERRPRT 

86440CPR 

86450CPR FILES USED — 

86460CPR ISINGF (S) SINDA GEOMETRY FILE (INPUT) 

86470CPR INTITS (S) INTERFACE FILE (OUTPUT) 

86480CPR 

86490CPR LIBRARIES ACCESSED — 

86500CPR NONE 

86510CPR 

86520CPR LOCAL VARIABLES — 

86530CPR DUMMY - USED IN 'DUMMY READ' TO POSITION 'ISINGF' AT 

86540CPR BEGINNING OF COORDINATE DATA 

86550CPR N - ARRAY OF NODE NUMBERS DEFINING CONNECTIVITY 

86560CPR OF 'IGELM(I)' 

86570CPR NODE - CURRENT NODE (CORNER POINT) NUMBER 

86580CPR X - X-COORDINATE OF 'NODE' 

86590CPR Y - Y-COORDINATE OF 'NODE' 

86600CPR 2 - Z-COORDINATE OF 'NODE' 

86610CPR 

86620CPR SPECIAL COMMENTS ABOUT THIS ROUTINE — 

86630CPR NONE 

86640CPR 

86650CPR SPECIAL REMARKS/INSTRUCTIONS — 

86660CPR HOST * H6000 

86670CPR 

86680CPR 

86690CPR* BB>**8a>CSI3tSCa*aiSBSCBB«SBatBC8X£B«asaBSCS>XCBXKi:K 

87490C .. END OF SUBROUTINE 'WRTGEO' 

87520 END 


GE/TRANCITS 


WRTGEO 


PAGE* 82 


159 


TITLE* 


DATE *07/15/85 


TIME * 23 42 


PAGE=83 


87560 

87600CPR 

87610CPR*** 

87620CPR 

87630CPR 

87640CPR 

87650CPR 

87660CPR 

87670CPR 

87680CPR 

87690CPR 

87700CPR 

87710CPR 

87720CPR 

87730CPR 

87740CPR 

87750CPR 

87760CPR 

87770CPR 

87780CPR 

87790CPR 

87800CPR 

87810CPR 

87820CPR 

87830CPR 

87840CPR 

87850CPR 

87860CPR 

87870CPR 

87880CPR 

87890CPR 

87900CPR 

87910CPR 

87920CPR 

87930CPR 

87940CPR 

87950CPR 

87960CPR 

87970CPR 

87980CPR 

87990CPR 

88000CPR 

88010CPR 

88020CPR 

88030CPR 

88040CPR 

88050CPR*** 

88840C 

88870 


SUBROUTINE WRTEMP 


SSBSaSBSBBBSBBSaSSSBSSSBBSSSaSBSSBBSSBSBBBSBBSS 


ROUTINE— WRTEMP 

DATE/WRITTEN BY— 01/11/84 RN PITTMAN 

DATE/REVISED BY — 

FUNCTION/PURPOSE — 

THIS ROUTINE WRITES THE TIME AND TEMPERATURE DATA ON THE 
'INTERFACE FILE' 

CALLING ARGUMENTS — 

NAME ATTRIBUTES DEFINITION 

NONE 


COMMONS USED — 

WORK —STORES INTERNAL ARRAY DATA 
AFIL — FILE CODES (MAIN) 

KNTDAT -- COUNTER DATA (MAIN) 

FUNCTIONS/ROUTINES CALLED — 

ERRPRT 

FILES USED — 

IHTIN (S) SINDA OUTPUT FILE (INPUT) 

INTITS (S) INTERFACE FILE (OUTPUT) 

LIBRARIES ACCESSED — 

NONE 

LOCAL VARIABLES -- 

KNT - POSITION IN THE 'TEMP' ARRAY OF THE NEXT GROUP OF 
SIX TEMPERATURES TO BE WRITTEN TO 'INTITS' 

TIME - CURRENT VALUE OF TIME READ FROM 'IHTIN' FOR THIS 
TEMPERATURE BLOCK 

SPECIAL COMMENTS ABOUT THIS ROUTINE — 

NONE 

SPECIAL REMARKS/INSTRUCTIONS — 

HOST « H6000 


*»tcBCEEisiEn3=icB*c=*ctsiisEEaE3tairiat***i:cE:Er 

END OF SUBROUTINE ' WRTEMP 9 
END 
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88910 SUBROUTINE ERRPRT (IROUT> 

88950CPR 

88960CPR CS3S = EBeeescss SBSX3ssrsss3sizss333=e:sssersss2S£:ss£= 

88970CPR 

88980CPR 

88990CPR ROUTINE — ERRPRT 

89000CPR 

89010CPR DATE/WRITTEN BY-- 01/12/84 RN PITTMAN 

89020CPR DATE/REVISED BY-- 

89030CPR 

89040CPR FUNCTION/PURPOSE-- 

89050CPR THIS ROUTINE PROCESSES ALL FATAL ERROR CONDITIONS BY 

89060CPR PRINTING AN EXPLANATION OF THE ERROR ON THE HARDCOPY 

89070CPR OUTPUT FILE AND TERMINATING PROGRAM EXECUTION 

89080CPR 

89090CPR CALLING ARGUMENTS — 

89100CPR NAME ATTRIBUTES DEFINITION 

89110CPR - 

89 120CPR IROUT INTEGER/INPUT ROUTING VARIABLE FOR APPROPRIATE 

89130CPR ERROR MESSAGE 

89140CPR 

89150CPR COMMONS USED -- 

89160CPR AFIL -- FILE CODES (MAIN) 

89170CPR 

89 180CPR FUNCTIONS/ROUTINES CALLED ~ 

89 190CPR NONE 

89200CPR 

89210CPR FILES USED -- 

89220CPR IERF (S) HARDCOPY FILE (OUTPUT) 

89230CPR 

89240CPR LIBRARIES ACCESSED -- 

89250CPR NONE 

89260CPR 

89270CPR LOCAL VARIABLES — 

89280CPR NONE 

89290CPR 

89300CPR SPECIAL COMMENTS ABOUT THIS ROUTINE -- 

89310CPR NONE 

89320CPR 

89330CPR SPECIAL REMARKS/INSTRUCTIONS — 

89340CPR HOST * H6000 

89350CPR 

89360CPR 

89370CPR** »*i=s««s«a=*«**«n««i 883 **i*ar**Bi«sEs*Bi 3 i«ie»**** 

91950 END 
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9 1 970C sssarssssscasaass s * * = = * = = = = * = s = s = = 

91990 SUBROUTINE TEMFOR( ISPOC ) 

92010C = * S5S ** ses = ES *s»*=ss===s=a===a=5==a*=ssssc=sr=sssa== 

92030CPR 

92040CPR* a38S3SSSSSS8SSSSB38S3:tSSeii:S3S8B8SSSSS888SS8:SSSI8 

92050CPR 

92060CPR 

92070CPR ROUTINE-- TEMFOR 

92080CPR 

92090CPR DATE/WRITTEN BY-- 01/14/84 RJ MAFFEO 

92100CPR DATE/REVISED BY-- 

92110CPR 

92 120CPR FUNCTION/PURPOSE-- 

92130CPR THIS ROUTINE CALLS ALL THE REQUIRED ROUTINES TO FORMAT THE 
92 140CPR OUTPUT TEMPERATURES INTO FORMED NEEDED FOR THE STRESS ANALYSIS 
92150CPR PROGRAM 
92 160CPR VALID CODES ARE 

92170CPR 0 --- NEUTRAL TEMPERATURE OUTPUT FILE ONLY 

92 180CPR 1 NASTRAN FORMATTED TEMPERATURES 

92190CPR 2,3 — MARC FORMATTED TEMPERATURES 

92200CPR 

92210CPR CALLING ARGUMENTS — 

92220CPR NAME ATTRIBUTES DEFINITION 

92230CPR 

92240CPR ISPOC (I) STRESS PROGRAM OUTPUT CODE 

92250CPR 

92260CPR COMMONS USED -- 
92270CPR 

92280CPR AFIL 

92290CPR SIZE 

92300CPR 

92310CPR FUNCTIONS/ROUTINES CALLED — 

92320CPR 

92330CPR TPNAST 

92340CPR MARCTO 

92350CPR 

92360CPR FILES USED -- 
92370CPR 

92380CPR IERF S 

92390CPR 

92400CPR LIBRARIES ACCESSED -- 

92410CPR NONE 

92420CPR 

92430CPR 

92440CPR LOCAL VARIABLES — 

92450CPR 

92460CPR SPECIAL COMMENTS ABOUT THIS ROUTINE 

92470CPR 

92480CPR 

92490CPR SPECIAL REMARKS/INSTRUCTIONS — 

92500CPR HOST * H6000 

92510CPR 

92520CPR 

92530CPR* *3IBMIBiaS3ailIS«S3aaai3R8B«||SSStlsSXSS«SB8SHCC=: 

93110 END 
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93140 SUBROUTINE TPNAST ( NUMNOD , NUMTMP ) 

93 170CPR 

93 180CPR BBE3 * s =* s * ! s»**se==**«ssss=t*c==sssrsssBsss2snc*=ss3ss 

93 190CPR 
93200CPR 

93210CPR ROUTINE-- TPNAST 

93220CPR 

93230CPR DATE/WRITTEN BY-- 01/14/84 RD MCCLAIN 

93240CPR DATE/REVISED BY — 

93250CPR 

93260CPR FUNCT I ON/PURPOSE -- 

93270CPR 

93280CPR WRITE NASTRAN BULK DATA TEMPERATURE CARDS 

93290CPR 

93300CPR 

93310CPR CALLING ARGUMENTS — 

93320CPR NAME ATTRIBUTES DEFINITION 

93330CPR 

93340CPR NUMNOD INPUT NUMBER OF NODES ON NEUT TEMP FILE 

93350CPR NUMTMP INPUT NUMBER OF TEMPERATURE SETS 

93360CPR 

93370CPR COMMONS USED -- 
93380CPR 

93390CPR AFIL — FILE CODES 

93400CPR 

93410CPR FUNCTIONS/ROUTINES CALLED -- 

93420CPR 

93430CPR 

93440CPR FILES USED — 

93450CPR 

93460CPR IOUTF S NEUTRAL OUTPUT FILE 

93470CPR ITEMP S NASTRAN BULK DATA FILE 

93480CPR 

93490CPR LIBRARIES ACCESSED -- 

93500CPR NONE 

93510CPR 

93520CPR 

93530CPR LOCAL VARIABLES — 

93540CPR 

93550CPR SPECIAL COMMENTS ABOUT THIS ROUTINE 

93560CPR 

93570CPR 

93580CPR SPECIAL REMARKS/INSTRUCTIONS -- 

93590CPR HOST * H6000 

93600CPR 

93610CPR 

93620CPR»»«»>»»««»«»*»»««».««.....«...». s ,..,« B « e 
93910 END 
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93930C CS = S1::S::SSBS3C3 = SS33SS = S&SS3n 3BB3SSS838S338888SE3B3338838 

93950 SUBROUTINE I3DSF( CMNMX , CRD , XP , YP , ZP , NSP , 

93970CPR 

93980CPR==== =*=33833=333838=33333333338333333833333=333333333 

93990CPR 

94000CPR 

94010CPR ROUTINE— I3DSF 

94020CPR 

94030CPR DATE/WRITTEN BY— 11/28/83 R J MAFFEO 

94040CPR DATE/REVISED BY— 

94050CPR 

94060CPR FUNCTION/PURPOSE— 

94070CPR THIS ROUTINE COMPUTES THE VALUES OF THE LOCAL COORDINATES 

94080CPR OF A STRESS POINT WRT THE VERTICES OF THE 8-NODED HEAT 

94090CPR TRANSFER ELEMENT . I F THE STRESS POINT LIES IN THE ELM 
94100CPR IT COMPUTES THE 8 WEIGHTING COEFFICIENTS FOR THIS POINT 

94110CPR IF THE POINT LIES OUTSIDE THE ELM IT COMPUTES THE SURFACE 

94120CPR COORDINATES OF THE POINT PROJECTED IN THE LOCAL COORD SYSTEM 
94 130CPR AND COMPUTES THE WEIGHTING COEFFS BASED ON THESE SURFACE 
94 140CPR COORDINATES IT ALSO COMPUTES THE DISTANCE THE POINT IS AWAY 
94 150CPR FROM THE SURFACE 
94160CPR 

94170CPR CALLING ARGUMENTS — 


94180CPR 

NAME 

ATTRIBUTES 

DEFINITION 


94 190CPR 






94200CPR 

CMNMX 

a) 


ARRAY WITH MIN AND MAX WINDOW COORDS 

94210CPR 

CRD 

(i) 


ARRAY WITH ELM VERTICE COORDS 


94220CPR 

XP 

(i) 

» 



94230CPR 

YP 

(i) 

>>» 

STRESS POINT COORDS 


94240CPR 

ZP 

a) 

>> 



94250CPR 

NSP 

(i) 


STRESS POINT NAME 


94260CPR 

XLOC 

(0) 

>> 



94270CPR 

YLOC 

(0) 

>>» 

LOCAL COORDS OF STRESS POINT 


94280CPR 

ZLOC 

(0) 

>> 



94290CPR 

SF 

(0) 


WEIGHTING COEFFICIENTS 


94300CPR 

ICONV 

(0) 


CONVERGENCE PARAMETER 


94310CPR 




0 — DID NOT CONVERGE 


94320CPR 




1 CONVERGED TO POINT IN THE 

ELM 

94330CPR 




2 — - CONVERGED TO POINT OUT OF 

ELM 

94340CPR 

IPAS 

(I) 


SEARCH TYPE PARAMETER 


94350CPR 




-1 — - NO REARRANGE.NO ROTATION 


94360CPR 




NO SURFACING 


94370CPR 




0 NO REARRANGE. ROTATE 


94380CPR 




NO SURFACING 


94390CPR 




1 REAARANGE. ROTATE 


94400CPR 




PERFORM SURFACING 


94410CPR 

DIST 

(0) 


DISTANCE FORM OUTSIDE POINT TO 

SURFACE 

94420CPR 

XSUR 

(0) 

» 



94430CPR 

YSUR 

(0) 

»>> 

SURFACE COORDINATES 


94440CPR 

ZSUR 

(0) 

>> 




94450CPR 

94460CPR COMMONS USED — 

94470CPR 

94480CPR AFIL 

94490CPR 

94500CPR FUNCTIONS/ROUTINES CALLED — 
94510CPR 

94520CPR REAARY 

94530CPR SHAPFV 

94540CPR CHKPT 

94550CPR JACBCK 

94560CPR SFDIST 

94570CPR 

94580CPR FILES USED — 

94590CPR 

94600CPR 

94610CPR LIBRARIES ACCESSED — 
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PAGE=88 


94620CPR 

94630CPR 

94640CPR 

94650CPR 

94660CPR 

94670CPR 

94680CPR 

94690CPR 

94700CPR 

94710CPR 

94720CPR 

94730CPR 

94740CPR 

94750CPR 

94760CPR 

94770CPR 

94780CPR 

94790CPR 

94800CPR ! 

97960 


NONE 


LOCAL VARIABLES — 

I POINT ACCURACY INDICATOR 

0 -- STRESS POINT NOT RECALCULATED ACCURATELY 

1 -- STRESS POINT RECALCULATED ACCURATELY 
IJACB — JACOBIAN CHECK PARAMETER 

0 — JACOBIAN NOT POSITIVE 

1 — JACOBIAN IS POSITIVE 

SPECIAL COMMENTS ABOUT THIS ROUTINE 


SPECIAL REMARKS/INSTRUCTIONS 
HOST * H6000 


END 
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TITLE= 
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PAGE=89 


97980C=== 

98000 

98020CPR 

98030CPR* 

98040CPR 

98050CPR 

98060CPR 

98070CPR 

98080CPR 

98090CPR 

98 100CPR 

981 10CPR 

98120CPR 

98130CPR 

98140CPR 

98150CPR 

98160CPR 

98170CPR 

98 180CPR 

98 190CPR 

98200CPR 

98210CPR 

98220CPR 

98230CPR 

98240CPR 

98250CPR 

98260CPR 

98270CPR 

98280CPR 

98290CPR 

98300CPR 

98310CPR 

98320CPR 

98330CPR 

98340CPR 

98350CPR 

98360CPR 

98370CPR 

98380CPR 

98390CPR 

98400CPR 

98410CPR 

98420CPR 

98430CPR 

98440CPR 

98450CPR 

98460CPR 

98470CPR 

98480CPR 

98490CPR 

98500CPR 

98510CPR 

98520CPR 

98530CPR 

98540CPR 

98550CPR 

98560CPR 

98570CPR 

98580CPR 

98590CPR 

98600CPR 

98610CPR 

99530 


SUBROUTINE REAARY ( X , Y , Z , XP , YP , ZP , SF , ICLOS . ITRYC , DTOCN, XL B YL , 2L 


ROUTINE-- REAARY 

DATE/WRITTEN BY-- 11/28/83 R U MAFFEO 

DATE/REVISED BY-- 

FUNCT I ON/PURPOSE -- 

THIS ROUTINE REARRANGES THE CONNECTIVITY COORDS AND 
THE SHAPE FUNCTION ARRAYS IN AN ATTEMPT TO GET AN 
ORDER THAT WILL ENHANCE CONVERGENCE OF THE INVERSE 
SHAPE FUNCTION ROUTINE ( I 3DSF ) IT ALSO COMPUTES THE 
DISTANCE FROM THE STRESS POINT TO THE NEAREST ELM 
CORNER OPTIONALY IT WILL M UN-REARRANGE " THE 
SHAPE FUNCTION ARRAY AND THE LOCAL COORDS 

CALLING ARGUMENTS — 


NAME 

ATTRIBUTES 

DEFINITION 

X 

(I/O) 

>> 


Y 

(I/O) 

>>>> 

CONNECTIVITY COORD ARRAYS 

Z 

(I/O) 

>> 


XP 

(I) 

>> 


YP 

(I) 

>>>> 

COORDS OF STRESS POINTS 

ZP 

(I) 

>> 


SF 

(I/O) 


SHAPE FUNCTION ARRAY 

ICLOS 

(I/O) 


VALUE OF NEW I JOINT 

ITRYC 

(I) 


CONTROL ARGUMENT 

-1 -- REARRANGE SHAPE FUNCTION ARRAY 
0 -- PICK I JOINT AS JOINT NEAREST TO 
STRESS POINT 

N>0 -- USE VALUE OF N AS I JOINT 

DTOCN 

(0) 


DISTANCE TO NEAREST CORNER 

XL 

(I/O) 

»> 


YL 

(I/O) 

>»»> 

VALUES OF LOCAL COORDINATES 

ZL 

COMMONS USED 

(I/O) 

>>> 



FUNCTIONS/ROUTINES CALLED — 
S0RT2 

FILES USED -- 


LIBRARIES ACCESSED -- 
NONE 


LOCAL VARIABLES — 

SPECIAL COMMENTS ABOUT THIS ROUTINE 

NOTE THE IJ ARRAY IS USED TO RELATE THE REARRANGED 
CONNECTIVITY ARRAY TO THE JOINT PICKED AS THE I JOINT 

END 
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99S50C atSBSSB * 8BS88Bsa 8*»>*8»eSS«ISCSSSES8SS3«ISS8: 

99570 SUBROUTINE S0RT2 ( SEEDS . FOLLO . LAX ) 

99590 C* se = 3 Sa ====**=»== 3 isss=ss=s= 5 a*sa=ass 3 rsc 5 sssssEi;ssse 

99610CPR ROUTINE — S0RT2 

99620CPR 

99630CPR DATE/WRITTEN BY-- T?**?????? GE DOUBLE SORT 

99640CPR DATE/REVISED BY — 

99650CPR 

99660CPR FUNCTION/PURPOSE— 

99670CPR 

99680CPR THIS ROUTINE SORTS TWO 1 DIMENSIONAL ARRAYS 
99690CPR 

99700CPR CALLING ARGUMENTS — 

99710CPR NAME ATTRIBUTES DEFINITION 

99720CPR 

99730CPR SEEDS (I/O) PRIMARY ARRAY TO BE SORTED 

99740CPR FOLLO (I/O) ARRAY TO BE SORTED LIKE THE PRIMARY 

99750CPR LAX (I) NUMBER OF ITEMS IN ARRAYS 

99760CPR 

99770CPR COMMONS USED -- NONE 
99780CPR 

99790CPR FUNCTIONS/ROUTINES CALLED -- NONE 
99800CPR 

99810CPR FILES USED -- NONE 
99820CPR 

99830CPR LIBRARIES ACCESSED -- 

99840CPR NONE 

99850CPR 

99860CPR LOCAL VARIABLES -- 
99870CPR 

99880CPR SPECIAL COMMENTS ABOUT THIS ROUTINE 

99890CPR 

99900CPR 

99910CPR SPECIAL REMARKS/INSTRUCTIONS — 

99920CPR HOST * H6000 

99930CPR 

99940CPR 

99950CPR> a * S33 > C33 >»>*lt«BBaaiBXS«iaC8E3assSE>C*KBSBS*ICBJIS 

100260 END 
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1 00280C S**B3SSSS3SSBSZS3XSSSSSSIZBSrSSHSBXB8B38S«BBaESS8SB3 

100300 SUBROUTINE SFDIST( SF , DIST , XNC , YNC , ZNC , XSUR , YSUR , ZSUR , XP . YP , ZP , 

1 003 30C 83SSS8S8 88 8&SSS«SSSSS8B&SBSSS3S3 8E3SSSaS28SSSC3S3 = BS 

100340CPR 

100350CPR”««”““””»”” = ””»”" = * 8 Sas . = 88». ssa3 = 

100360CPR 

070^01? 

100380CPR ROUTINE" SFDIST 

100390CPR 

100400CPR DATE/WRITTEN BY— 11/15/83 R 0 MAFFEO 

1004 10CPR DATE/REVISED BY-- 

100420CPR 

100430CPR FUNCTION/PURPOSE — 

100440CPR 

100450CPR THIS ROUTINE COMPUTES THE DISTANCE FROM STRESS 

100460CPR POINTS THAT ARE OUTSIDE A HEAT TRANSFER ELM TO THE 

100470CPR SURFACE OF THE HEAT TRAN ELM IT ALSO COMPUTES THE 
100480CPR SURFACE COORDIANTES THIS SURFACE POINT IS FOUND BY 
100490CPR 'TRUNCATING' THE LOCAL COORDINATES TO 1 OR -1 
100500CPR 
1005 10CPR 
100520CPR 

100530CPR CALLING ARGUMENTS -- 


1OO540CPR 

NAME 

ATTRIBUTES 

DEFINITION 

100550CPR 

100560CPR 

SF 

(i) 


SHAPE FUNCTION ARRAY 

100570CPR 

DIST 

(0) 


COMPUTED DISTANCE 

100580CPR 

XNC 

(1) 

>>> 


100590CPR 

YNC 

(I) 

>>>>>>>>> 

LOCAL VALUES OF STRESS POINT 

100600CPR 

ZNC 

(I) 

»> 


1006 10CPR 

XSUR 

(0) 

»> 


100620CPR 

YSUR 

(0) 

>>>>>>>>> 

COORDS OF SURFACE POINT 

100630CPR 

ZSUR 

(0) 

>» 


100640CPR 

XP 

(I) 

>>> 


100650CPR 

YP 

(I) 

A 

A 

A 

A 

A 

A 

A 

A 

A 

COORDS OF STRESS POINT 

100660CPR 

ZP 

(I) 

»> 


100670CPR 

NNPE 

(I) 


NUMBER OF NODES PER ELM 

100680CPR 

CRD 

(I) 


COORDS OF VERTICES OF ELM 

100690CPR 

100700CPR 

COMMONS USED 





1007 10CPR 
100720CPR 

100730CPR FUNCTIONS/ROUTINES CALLED -- 
100740CPR 

100750CPR SHAPFV 
100760CPR 

100770CPR FILES USED — 

100780CPR 

100790CPR 

100800CPR LIBRARIES ACCESSED -- 

1008 10CPR NONE 

100820CPR 

100830CPR 

100840CPR LOCAL VARIABLES — 

100850CPR 

100860CPR SPECIAL COMMENTS ABOUT THIS ROUTINE 

100870CPR 

100880CPR 

100890CPR SPECIAL REMARKS/INSTRUCTIONS -- 

100900CPR HOST * H6000 

100910CPR 

100920CPR 

100930CPR asc8liassaSBBBeS8a ** SIBBI * cscsa ** ss8SSB,,BB8B * ,KaEzt>|iB3 
101290 END 
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10131 0C * s* «BSSSBBses*BBBSssssssBsssBssss 

101330 SUBROUTINE CHKPT ( SF , X , Y , Z , XS , YS , ZS , IPOINT ) 

101340CPR 


1013 50C PR 3 s = c 5crsss«; = ssa*ci»«r = r = *xzrirsiBfi = 2 c = 5ii:rt = = = cssiis 


101 360CPR 

101370CPR 

101380CPR 

101390CPR 

101400CPR 

101410CPR 

101420CPR 

101430CPR 

101440CPR 

101450CPR 

101460CPR 

101470CPR 

101480CPR 

101490CPR 

101500CPR 

1015 10CPR 

101520CPR 

101530CPR 

101540CPR 

101550CPR 

101560CPR 

101570CPR 

101580CPR 

101590CPR 

101600CPR 

101610CPR 

101620CPR 

101630CPR 

101640CPR 

101650CPR 

101660CPR 

101670CPR 

101680CPR 

101690CPR 

101 700CPR 

1017 10CPR 

101720CPR 

101730CPR 

101740CPR 

101750CPR 

101760CPR 

101770CPR 

101780CPR 

101790CPR 

101800CPR 

101810CPR 

101820CPR 

101830CPR 

101840CPR* 

101870C*** 

102110 


ROUTINE — CHKPT 

DATE/WRITTEN BY— 11/28/83 R J. MAFFEO 
DATE/REVISED BY— 

FUNCTION/PURPOSE — 

THIS ROUTINE RECOMPUTES THE STRESS POINT COORDS BASED 
ON THE VALUES OF THE LOCAL COORDINATES TO CHECK FOR 
ACCURACY 

CALLING ARGUMENTS — 


NAME 

ATTRIBUTES DEFINITION 

SF 

(I) 


ARRAY OF 

SHAPE FUNCTIONS 

X 

(I) 

» 



Y 

(I) 

>>>> 

ARRAY OF 

ELM CORNER COORDS 

Z 

(I) 

» 



XS 

(I) 

>> 



YS 

(I) 

»» 

ORIGINAL 

STRESS POINT COORDS 

ZS 

(I) 

>> 



IPOINT 

(0) 


ACCURACY 

INDCATOR 


0 — DID NOT REPRODUCE COORD 

1 — DID REPRODUCE COORD 


COMMONS USED — 


FUNCTIONS/ROUTINES CALLED 


FILES USED 


LIBRARIES ACCESSED — 
NONE 


LOCAL VARIABLES — 

SPECIAL COMMENTS ABOUT THIS ROUTINE 


SPECIAL REMARKS/INSTRUCTIONS — 
HOST * H6000 


,, ** I * ae * , ****4***a***fBan«sa»«*2ii*rri3i*B**BBi 

END 
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1 02 1 30C= = 3 = = = = = = = « = = 3 = = = = = * = = = e = = = s = = = = = s = = = s = 

102150 SUBROUTINE PDOFSF ( TETA , YNETA , ZETA , IGP . IND , PD ) 

102 190CPR 

102200CPR=== ==*»==*= ========= === s ====»==== =ass==a==s ======= --- 

1022 10CPR 
102220CPR 

102230CPR ROUTINE— PDOFSF 

102240CPR 

102250CPR DATE/WRITTEN BY-- 09/14/82 RJ MAFFEO 

102260CPR DATE/REVISED BY — 

102270CPR 

102280CPR FUNCTION/PURPOSE — 

102290CPR 

102300CPR THIS SUB COMPUTES THE DERIVATIVES OF THE 
102310CPR SHAPE FUNCTIONS WRT THE LOCAL SYSTEM 
102320CPR 

102330CPR CALLING ARGUMENTS — 

102340CPR NAME ATTRIBUTES DEFINITION 

102350CPR 

102360CPR TETA (I) > 

102370CPR YNETA (I) » LOCAL COORDINATES 

102380CPR ZETA (I) > 

102390CPR IGP (I) GAUSS POINT COUNTER 

102400CPR IND (I) 2D/3D INDICATOR 

1024 10CPR IND-0 >> 3D 

102420CPR IND= 1 THRU 6 >> 2D IND-FACE NUMBER 

102430CPR PD (0) ARRAY OF DERIVATIVES 

102440CPR 

102450CPR COMMONS USED — 

102460CPR 

102470CPR CORNDEF 

102480CPR 

102490CPR FUNCTIONS/ROUTINES CALLED — 

10250CCPR 

102510CPR 

102520CPR FILES USED -- 

102530CPR 

102540CPR 

102550CPR LIBRARIES ACCESSED — 

102560CPR NONE 

102570CPR 

102580CPR 

102590CPR LOCAL VARIABLES -- 
102600CPR 

1026 10CPR SPECIAL COMMENTS ABOUT THIS ROUTINE 

102620CPR 

102630CPR 

102640CPR SPECIAL REMARKS/INSTRUCTIONS — 

102650CPR HOST * H60O0 

102660CPR 

102670CPR 

1 02680C PR**** BB * s * ss= *** s * SB ***s8srsssss8asssaa*sssa*«s3B*sss 
103040 END 
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1 03060C = 3 = = » = = = = * = = = = = = a = = * = = = = = = c = 

103080 SUBROUTINE FND JAC ( X , Y , Z , P3D , X JAC ) 

103090CPR 

103100CPR SS3 s**=e*»==a8B*sssr=sstsa*«sa*Bar5sias*Bseet*as*stBs 
1031 10CPR 
103 120CPR 

103 130CPR ROUTINE-- FNDJAC 

103 140CPR 

103 150CPR DATE/WRITTEN BY— 11/28/83 R.J MAFFEO 

103160CPR DATE/REVISED BY — 

103170CPR 

103180CPR FUNCTION/PURPOSE — 

103190CPR THIS ROUTINE COMPUTES THE VALUE OF THE DETERMINANT OF 
103200CPR THE JACOBIAN EVALUATED AT THE STRESS POINT 
1032 10CPR 


103220CPR 

CALLING 

ARGUMENTS — 


103230CPR 

NAME 


ATTRIBUTES 

DEFINITION 

103240CPR 





103250CPR 

X 

(I) 

>> 


103260CPR 

Y 

(I) 

»» ARRAY 

OF ELM CORNER COORDS 

103270CPR 

Z 

(I) 

>> 


103280CPR 

P3D 

(I) 

ARRAY 

WITH PARTIAL DERIVATIVES OF 

103290CPR 



LOCAL 

COORD WRT THE GLOBAL COORD 

103300CPR 

X JAC 

(0) 

VALUE 

OF DETERMINANT OF JACOBIAN 

1033 10CPR 






103320CPR COMMONS USED -- 

103330CPR 

103340CPR 

103350CPR FUNCTIONS/ROUTINES CALLED — 

103360CPR 

103370CPR 

103380CPR FILES USED -- 

103390CPR 

103400CPR 

103410CPR LIBRARIES ACCESSED — 

103420CPR NONE 

103430CPR 

103440CPR 

103450CPR LOCAL VARIABLES — 

103460CPR 

103470CPR SPECIAL COMMENTS ABOUT THIS ROUTINE 

103480CPR 

103490CPR 

103500CPR SPECIAL REMARKS/INSTRUCTIONS -- 


1035 10CPR HOST * H6000 

103520CPR 

103530CPR 

103540CPR B> * 3 >»*** 3 »»BS>aBiei8>iiBiaie««BBBic8nB8i 

103570C BSS8 *BB«s»>aBBc«>sassaess*Bsc«sc«esBscsn*s«ss8sxs*ss33a«3saca 
103810 END 
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TITLE* 


DATE *07/15/85 


TIME = 23 42 


PAGE=95 


1 03830C==*== ==*======================SK=;=======r================== 

103850 SUBROUTINE SHAPFV ( SF f XNC , YNC , 2NC ) 

103860CPR 

1 03870CPR== Sa>BS8S8SB33SSS:SBSe3SSSt8S33888SSBSSSSSS8SSSBSS&8S 

103880CPR 

103890CPR 

103900CPR ROUTINE-- SHAPFV 

1039 10CPR 

103920CPR DATE/WRITTEN BY-- 11/28/83 R J MAFFEO 

103930CPR DATE/REVISED BY-- 

103940CPR 

103950CPR FUNCTION/PURPOSE -- 

103960CPR THIS ROUTINE FINDS THE VALUES OF THE SHAPE FUNCTIONS 

103970CPR FOR THE 8-NODED ISOPARAMETRIC ELM BASED ON THE VALUES 

103980CPR OF THE LOCAL COORDINATES 
103990CPR 

104000CPR CALLING ARGUMENTS — 

104010CPR NAME ATTRIBUTES DEFINITION 

104020CPR 

104030CPR SF (0) ARRAY WITH SHAPE FUNCTIONS 

104040CPR XNC (I) » 

104050CPR YNC (I) >>>> LOCAL COORDS 

104060CPR ZNC (I) » 

104070CPR 

104080CPR COMMONS USED -- 

104090CPR 

104 100CPR 

104110CPR FUNCTIONS/ROUTINES CALLED -- 

104 120CPR 

104130CPR 

104 140CPR FILES USED -- 

104150CPR 

104 160CPR 

1041 70CPR LIBRARIES ACCESSED -- 
104 180CPR NONE 
104 190CPR 
104200CPR 

104210CPR LOCAL VARIABLES — 

104220CPR 

104230CPR SPECIAL COMMENTS ABOUT THIS ROUTINE 

104240CPR 

104250CPR 

104260CPR SPECIAL REMARKS/INSTRUCTIONS — 

104270CPR HOST * H6000 

104280CPR 

104290CPR 

104 300CP R C£BB * SX*BKCBSB«XaBXSaBXS*SeSBCSBK = B*B*»CS*KX«BatS:S*SC 
1 04330C** *s»*s*aseB»*3*aB»«*B*»»**8sa**a«maBaBasBBsa»**B*Ba«*BB*8 
104490 END 
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TITLE* 


DATE *07/1 5/85 


TIME = 23 42 


PAGE*96 


1045100**= 

104530 

104540CPR 

104550CPR* 

104560CPR 

104570CPR 

104580CPR 

104590CPR 

104 6000 PR 

1046 10CPR 

104620CPR 

104630CPR 

104640CPR 

104650CPR 

104660CPR 

104670CPR 

104680CPR 


:=s=ca:a«;3Sisa3sac===ra«35scrsrc*=Ecsasj=*=r=Es;=: 

SUBROUTINE JACBCK( XL * YL , 2L , X , Y , 2 * X JAC , I JACB ) 


ROUTINE — 


JACBCK 


DATE/WRITTEN BY— 11/28/83 
DATE/REVISED BY — 


R . J MAFFEO 


FUNCTION/PURPOSE-- 

THIS ROUTINE CALLS THE REQUIRED ROUTINES TO CHECK THE 
VALUE OF THE DETERMINANT OF THE JACOBIAN AT THE STRESS 
POINT 

CALLING ARGUMENTS -- 


104690CPR 

NAME 

ATTRIBUTES DEFINITION 

104700CPR 
1047 10CPR 

XL 

(I) 

>> 


104720CPR 

YL 

(I) 

»» 

LOCAL COORDS 

104730CPR 

2L 

(I) 

>> 


104740CPR 

X 

(I) 

» 


104750CPR 

Y 

(I) 

>>>> 

ARRAY OF ELM CORNER COORDS 

104760CPR 

2 

(I) 

>> 

104770CPR 

XJAC 

(0) 


VALUE OF DETERMINANT OF JACOBIAN 

104780CPR 

I JACB 

(0) 


INTERGER VARIABLE 

104790CPR 
104800CPR 
1048 10CPR 
104820CPR 

COMMONS USED 



0 — JACOBIAN NOT POSITIVE 

1 — JACOBIAN IS POSITIVE 


104830CPR 

104840CPR 

104850CPR 

104860CPR 

104870CPR 

104880CPR 

104890CPR 

104900CPR 

104910CPR 

104920CPR 

104930CPR 

104940CPR 

104950CPR 

104960CPR 

104970CPR 

104980CPR 

104990CPR 

105000CPR 

105010CPR 

105020CPR 

105030CPR 

105040CPR 

105050CPR 

105060CPR** 

1050900***= 

105250 


FUNCTIONS/ROUTINES CALLED — 

PDOFSF 

FNDJAC 

FILES USED — 


LIBRARIES ACCESSED 
NONE 


LOCAL VARIABLES — 

SPECIAL COMMENTS ABOUT THIS ROUTINE 


SPECIAL REMARKS/INSTRUCTIONS — 
HOST * H6000 


BKKKESSaca 


Esssssscsac 

iBS*BS8taBri 


END 
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TITLE* 


DATE *07/15/85 


TIME * 23 42 


PAGE=97 


1 05 2 7 0C sssssasssssasssssssssssssssssssssssssssssssssssssssss 

105290 SUBROUTINE RWSPCF( IREC , B , IRW) 

1 05 3 1 0C = = = = = = = = = = = = = = c = = = = = = = = = 3 = = = = = = = = = aa = = = = = = = = = = = = = = = = 
105330CPR 

105340CPR*“*“=*«“=«*******»«*****«* m*mb*==wmwm**m 

105350CPR 

105360CPR 

105370CPR ROUTINE— RWSPCF 

105380CPR 

105390CPR DATE/WRITTEN BY— 04/07/83 RJ MAFFEO 

105400CPR DATE/REVISED BY— 

105410CPR 

105420CPR FUNCTION/PURPOSE— 

105430CPR 

105440CPR THIS ROUTINE READS AND WRITES DATA TO THE 
105450CPR THE RANDOM STRESS POINT COEFF FILE 
105460CPR THIS FILE CONTAINS THE FOLLOWING INFO PER RECORD 
105470CPR - STRESS NODE NAME (1) 

105480CPR - HEAT TRAN ELM CONTAINING THIS STRESS POINT (1) 

105490CPR - CONNECTIVITY OF THE HEAT TRAN ELM (8) 

105500CPR - WEIGHTING COEFF FOR THIS STRESS NODE (8) 

105510CPR 

105520CPR 

105530CPR 

105540CPR CALLING ARGUMENTS — 

105550CPR NAME ATTRIBUTES DEFINITION 

105560CPR 

105570CPR IREC (I) FILE RECORD NUMBER 

105580CPR B (I/O) DATA ARRAY 

105590CPR IRW (I) READ/WRITE OPTION 

105600CPR (O-READ 1-WRITE) 

105610CPR 

105620CPR COMMONS USED — 

105630CPR 

105640CPR AFIL 

105650CPR 

105660CPR FUNCTIONS/ROUTINES CALLED — 

105670CPR 

105680CPR 

105690CPR FILES USED -- 
105700CPR 

1057 10CPR IRSPCF R 

105720CPR 

105730CPR LIBRARIES ACCESSED — 

105740CPR NONE 

105750CPR 

105760CPR 

105770CPR LOCAL VARIABLES — 

105780CPR 

105790CPR SPECIAL COMMENTS ABOUT THIS ROUTINE 

105800CPR 

105810CPR 

105820CPR SPECIAL REMARKS/INSTRUCTIONS — 

105830CPR HOST * H6000 

105840CPR 

105850CPR 

1 05860CPR* *88S883C88VaCltaSCSZIIZXXCCSSB3CBCB>tatseXZB>Xt3EaS 

106170 END 
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TITLE 


DATE *07/15/85 


TIME 


23 42 


PAGE*98 




106210 SUBROUTINE ELTRAN 

106230C 

106250CPR 

1 06260CPR* «*aasBaB5sssaaBssssBS«*sssss*a*saaas**ss«s*a**sss*s= 

106270CPR 

106280CPR 

106290CPR ROUTINE— ELTRAN 

106300CPR 

106310CPR DATE/WRITTEN BY— 01/27/84 RJ MAFFEO 

106320CPR DATE/REVISED BY — 

106330CPR 

106340CPR FUNCTION/PURPOSE — 

106350CPR 

106360CPR Validate the values input for IECODE and IACODE then assign 

106370CPR the values to NNPE and NTPE depending on the values of 

106380CPR IECODE and IACODE 

106390CPR 

106400CPR 

106410CPR 

106420CPR 

106430CPR CALLING ARGUMENTS — 

106440CPR NAME ATTRIBUTES DEFINITION 

106450CPR 

106460CPR 

106470CPR COMMONS USED — 

106480CPR 

106490CPR AFIL 

106500CPR CNTLFL 

1065 10CPR 

106520CPR FUNCTIONS/ROUTINES CALLED — 

106530CPR 

106540CPR 

106550CPR FILES USED — 

106560CPR 

106570CPR I ERF S 

106580CPR 

106590CPR LIBRARIES ACCESSED — 

106600CPR NONE 

106610CPR 

106620CPR 

106630CPR LOCAL VARIABLES — 

106640CPR 

106650CPR SPECIAL COMMENTS ABOUT THIS ROUTINE 

106660CPR 

106670CPR 

106680CPR SPECIAL REMARKS/INSTRUCTIONS -- 

106690CPR HOST * H6000 

106700CPR 
1067 10CPR 

1 06 7 20C PR* eBS ** ei **a**8*a*BBS*t*aaascaa8*SBi«Bt«*ra«»*«ai»f *# 

107330 END 
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TITLE* DATE *07/15/85 TIME * 23 42 PAGE=99 

107350C X3SSSSSSS3SSSSSSSXSSSSS=SSSSS3SSSSa txs*=sssss=sssss=s * 

107370 SUBROUTINE INITCM 

1 07 390C Sa - SS3333333 SBSB = SSBSBSSSSSSSSSSSSSSS = =S3SBS = 3BSSSSS 

1074 10CPR 

107420CPR ssas = sssssiess3 = ssssssasEssssssssaaa3 * aaaa s aaaaaaaa3 s3s 

107430CPR 

107440CPR 

107450CPR ROUTINE— INITCM 

107460CPR 

107470CPR DATE/WRITTEN BY-- 01/31/84 RJ MAFFEO 

107480CPR DATE/REVISED BY — 

107490CPR 

107500CPR FUNCTION/PURPOSE— 

107510CPR 

107520CPR THIS ROUTINE INITIALIZES THE CRNDEF COMMON 
107S30CPR THIS COMMON IS USED IN CALCRD AND I3DSF ROUTINES 
107540CPR 

107550CPR CALLING ARGUMENTS — 

107560CPR NAME ATTRIBUTES DEFINITION 

107570CPR 

107580CPR 

107590CPR COMMONS USED -- 
107600CPR 

1076 10CPR CRNDEF 

107620CPR 

107630CPR FUNCTIONS/ROUTINES CALLED — 

107640CPR 

107650CPR 

107660CPR FILES USED -- 

107670CPR 

107680CPR 

107690CPR LIBRARIES ACCESSED — 

1O770OCPR NONE 
1077 10CPR 
107720CPR 

107730CPR LOCAL VARIABLES — 

107740CPR 

107750CPR SPECIAL COMMENTS ABOUT THIS ROUTINE 

107760CPR 

107770CPR 

107780CPR SPECIAL REMARKS/INSTRUCTIONS — 

107790CPR HOST = H6000 

107800CPR 

107810CPR 

107820 CPR* B *** 8 « B *** B »**»=* 3 a *=«*==a«**=se**=a=ta«Esa=s»s*» 3 = 

108390 END 
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1 084 1 0C ^BB3S3rB8SISCZS8SS(SSS3SSIlStCSSSE8BtBS(IESBSSBISB8SSSICtISSES 

108430 SUBROUTINE MARCTO ( I SPOC , NUMSP , NUMSOL ) 

108470CPR 

l08480CPR li 8SSBS3S33ESrSS8SSBIBSBISStBStrSSSS3S3SS8aSSSe 

108490CPR 

108500CPR 

108510CPR ROUTINE— MARCTO 

108520CPR 

108530CPR DATE/WRITTEN BY-- 03/08/84 RJ MAFFEO 

108540CPR DATE/REVISED BY — 

108550CPR 

108560CPR FUNCTION/PURPOSE — 

108570CPR 

108580CPR THIS ROUTINE CALLS ALL REQUIRED ROUTINES TO COMPUTE 
108590CPR THE MARC 8-NQDED ELM GAUSS POINT TEMPERATURES (OR THE 
108600CPR ELM CENTROIDAL TEMPERATURE) 

1086 10CPR 


108620CPR 

CALLING ARGUMENTS 

— 


108630CPR 

NAME 

ATTRIBUTES 

DEFINITION 

108640CPR 

108650CPR 

ISPOC 

(i) 

DESIRED ELM LOCATION 

108660CPR 

108670CPR 

108680CPR 

NUMSP 

(i) 

I5P0C*2 - 
ISPOC-3 - 
NUMBER OF 

-- CENTROIDAL TEMPERATURE 
-- 2X2X2 GAUSS POINT TEMPERATURES 
STRESS MODEL NODAL POINTS 

108690CPR 

NUMSOL 

(i) 

NUMBER OF 

TRANSIENT SOLUTIONS 


108700CPR 

1087 10CPR COMMONS USED — 

108720CPR 

108730CPR AFIL 

108740CPR WORK 

108750CPR 

108760CPR FUNCTIONS/ROUTINES CALLED — 
108770CPR 

108780CPR GPCOEF 

108790CPR GPTEMP 

108800CPR 

1088 10CPR FILES USED — 


108820CPR 

108830CPR ISNEF S 
108840CPR IOUTF S 
108850CPR ITEMP S 
108860CPR I ERF S 
108870CPR 


108880CPR LIBRARIES ACCESSED -- 

108890CPR NONE 

108900CPR 

108910CPR 

108920CPR LOCAL VARIABLES — 

108930CPR 

108940CPR SPECIAL COMMENTS ABOUT THIS ROUTINE 

108950CPR 

108960CPR 

108970CPR SPECIAL REMARKS/INSTRUCTIONS — 

108980CPR HOST * H6000 

108990CPR 

109000CPR 

109010CPR** **■■■*■•■■****««***««*«■*“*«** *■***» 

109950 END 
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TITLE- 


DATE =07/15/85 


TIME = 23 42 


PAGE*** 


1 09970C *«===» »**e* = s = *ss = * = ses = = *8»ss = ss* = »=sss3 = = s«**ss*:se = = *ss*s85s* = = 

109990 subroutine gpcoef(gpf) 

1 1 00 1 OC = a*3* = = ** = a = = s* = = = = = ss = c = se = * = c*ssBasss3 = * = = s sssssskbssssssssssss 

1 1OO30CPR 

1 1 004 OC PI? aB 5S2S2BssaBzsss**ssess*sBs«ssBsrsBssssss3asssrsss*as 
1 10050CPR 
1 10060CPR 

1 10070CPR ROUTINE— GPCOEF 

1 10080CPR 

1 10090CPR DATE/WRITTEN BY— 03/07/84 RJ MAFFEO 

1 10100CPR DATE/REVISED BY — 

1 101 10CPR 

1 101 20CPR FUNCTION/PURPOSE — 

1 10130CPR 

1 10140CPR THIS ROUTINE COMPUTES AND STORES THE WEIGHTING 
1 10150CPR COEFFS ASSOCIATED WITH THE 2X2X2 GAUSS POINTS 
1 10160CPR USED FOR THE MARC SOLID ELM( ELM TYPE 43) 

1 10170CPR 
1 10180CPR 
1 10190CPR 

1 10200CPR CALLING ARGUMENTS — 

110210CPR NAME ATTRIBUTES DEFINITION 

1 10220CPR 

1 10230CPR GPF (0) ARRAY WITH GAUSS PT COEFF 

1 10240CPR 

1 10250CPR COMMONS USED — 

1 10260CPR 
1 10270CPR 

1 10280CPR FUNCTIONS/ROUTINES CALLED — 

1 10290CPR 
1 10300CPR 

1 103 lOCPR FILES USED — 

1 10320CPR 
1 10330CPR 

1 10340CPR LIBRARIES ACCESSED — 

1 10350CPR NONE 
1 10360CPR 
1 10370CPR 

1 10380CPR LOCAL VARIABLES — 

1 10390CPR 

1 10400CPR SPECIAL COMMENTS ABOUT THIS ROUTINE 
1 10410CPR 
1 10420CPR 

110430CPR SPECIAL REMARKS/INSTRUCTIONS — 

1 10440CPR HOST * H6000 

1 10450CPR 
1 10460CPR 

1 10470CPR CKSSSS *Ct8tBB3S33SBBBt;BSBBSa:>B«SSI«BStS8BtSCBSEIS= 

110710 END 
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1 1 0730C ssssssssasssssKssessss sc ssssbsssssxsssssssssssskscbssssss 

110750 SUBROUTINE GPTEMP(TA . IE , GPF , ISPOC , TV) 

1 1 07 70C SSSSXSSSXSSSSSaCSSXSSSXKSSSSSSSSSSSBEBZBSSSSSSSSSSSSKSttSSS 

1 10790CPR 

110810CPR 
1 10820CPR 

1 10830CPR ROUTINE-- GPTEMP 

1 10840CPR 

1 10850CPR DATE/WRITTEN BY-- 03/17/84 RJ MAFFEO 

1 10860CPR DATE/REVISED BY — 

1 10870CPR 

1 10880CPR FUNCTION/PURPOSE— 

1 10890CPR 

1 10900CPR THIS ROUTINE COMPUTES THE VALUE OF THE ELM TEMPERATURE 
110910CPR FOR THE 8-NODED 30 MARC STRESS ELM IT WILL COMPUTE 
1 10920CPR TEMPS AT THE 8 2X2X2 GAUSS POINTS OR AT THE ELM 
1 10930CPR CENTROID 
1 10940CPR 
1 10950CPR 
1 10960CPR 

1 10970CPR CALLING ARGUMENTS -- 


1 10980CPR 
1 10990CPR 

NAME 

ATTRIBUTES 

1 DEFINITION 

1 1 1000CPR 

TA 

(I) 

NODAL TEMPERATURE ARRAY 

1 11010CPR 

IE 

(I) 

ELM CONNECTIVITY ARRAY 

1 1 1020CPR 

GPF 

(I) 

GAUSS PT COEFF ARRAY 

1 1 1030CPR 

ISPOC 

(I) 

TEMPERATURE LOCATION PARAMETER 

1 1 1040CPR 

TV 

(0) 

ELEMENTAL TEMPERATURE ARRAY 


1 1 1050CPR NOTE IF ISPOC-2 

1 1 1060CPR TV( 1 ) “CENTROIDAL TEMPERATURE 

1 1 1070CPR 

1 1 1080CPR COMMONS USED -- 
1 1 1090CPR 
111 100CPR 

1111 10CPR FUNCTIONS/ROUTINES CALLED — 

1111 20CPR 
111 130CPR 

111 140CPR FILES USED -- 

111 150CPR 

111160CPR 

111 170CPR LIBRARIES ACCESSED -- 
111 180CPR NONE 
111 190CPR 
1 1 1200CPR 

1 1 12 10CPR LOCAL VARIABLES — 

1 11220CPR 

1 1 1230CPR SPECIAL COMMENTS ABOUT THIS ROUTINE 
1 1 1240CPR 
1 1 1250CPR 

1 1 1260CPR SPECIAL REMARKS/INSTRUCTIONS — 

1 1 1270CPR HOST » H6000 

1 1 1280CPR 
1 1 1290CPR 

111 300CPR** •CSS«BSSZEaES=SXC8ISS£3Se£SS«CSaBBSm>>BBBeszBB33EE 

111590 END 


GE/TRANCITS 


GPTEMP 


PAGE* 102 LAST 


179 


APPENDIX D 


* Transfer Module Flow Diagram 

This appendix contains a program flow diagram of the transfer module. 
This diagram shows the general flow of the code and which subroutines are 
subordinate to each other. 


180 


TITLE- 


DATE -07/15/85 TIME = 23 58 PAGE- 1 


* * 

* THE SAVAM STACKER * 

* * 


THIS IS CALLING STACK STARTING FROM $$$ MAI 
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CO 

3 


AFIL 

C 

NEWSRC 

100011 

CALCRD 

CO 

3 


MAFCOM 

M 

NEWSRC 

100011 

CALCRD 

CO 

3 


ELMDAT 

C 

NEWSRC 

100011 

CALCRD 

CO 

3 


MAFCOM 

M 

NEWSRC 

100011 

CALCRD 

CO 

3 


SIZE 

C 

NEWSRC 

10001 1 

CALCRD 

CO 

3 


MAFCOM 

M 

NEWSRC 

100011 

CALCRD 

CO 

3 


WORK 

C 

NEWSRC 

100011 

CALCRD 

CO 

3 


CRNDEF 

c 

NEWSRC 

10001 1 

CALCRD 

SU 

5 

101 

GETPD 


NEWSRC 

10001 1 

CALCRD 

SU 

5 

102 

RWELMS 


NEWSRC 

10001 1 

CALCRD 

SU 

5 

103 

VOLM 


NEWSRC 

10O011 

CALCRD 

SU 

5 

104 

AREAD 


NEWSRC 

100011 

CALCRD 

SU 

5 

105 

REACRD 


NEWSRC 

100011 

CALCRD 

WR 

8 


YES 


NEWSRC 

10001 1 

GETPD 

SU 

1 


GETPD 


NEWSRC 

100012 

GETPD 

CO 

3 


CRNDEF 

c 

NEWSRC 

100012 

GETPD 

SU 

5 

101 

PARD 


NEWSRC 

100012 

GETPD 

SU 

5 

102 

GTGAUS 


NEWSRC 

100012 

VOLM 

SU 

1 


VOLM 


NEWSRC 

100013 

VOLM 

SU 

5 

101 

SHAPEF 


NEWSRC 

100013 

SHAPEF 

SU 

1 


SHAPEF 


NEWSRC 

100014 

SHAPEF 

CO 

3 


CRNDEF 

c 

NEWSRC 

100014 

PARD 

SU 

1 


PARD 


NEWSRC 

100015 

PARD 

CO 

3 


CRNDEF 

c 

NEWSRC 

100015 

AREAD 

SU 

1 


AREAD 


NEWSRC 

100016 

AREAD 

CO 

3 


CRNDEF 

c 

NEWSRC 

100016 

AREAD 

SU 

5 

101 

GTGAUS 


NEWSRC 

100016 

AREAD 

SU 

5 

102 

SHAPEF 


NEWSRC 

100016 

GTGAUS 

SU 

1 


GTGAUS 


NEWSRC 

100017 

REACRD 

SU 

1 


REACRD 


NEWSRC 

100018 

CORNCC 

SU 

1 


CORNCC 


NEWSRC 

100019 

CORNCC 

CO 

3 


MAFCOM 

' M • 

NEWSRC 

100019 

CORNCC 

CO 

3 


AFIL 

c 

NEWSRC 

100019 

CORNCC 

CO 

3 


MAFCOM 

M 

NEWSRC 

100019 

CORNCC 

CO 

3 


ELMDAT 

C 

NEWSRC 

100019 

CORNCC 

CO 

3 


MAFCOM 

M 

NEWSRC 

100019 

CORNCC 

CO 

3 


SIZE 

C 

NEWSRC 

100019 

CORNCC 

CO 

3 


MAFCOM 

M 

NEWSRC 

100019 

CORNCC 

CO 

3 


WORK 

C 

NEWSRC 

100019 

CORNCC 

CO 

3 


MAFCOM 

M 

NEWSRC 

100019 
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CORNCC 

CO 

3 


CNTLGM 

C 

NEWSRC 

100019 

CORNCC 

SU 

5 

101 

RWELMS 


NEWSRC 

100019 

CORNCC 

SU 

5 

102 

SORT 


NEWSRC 

100019 

CORNCC 

SU 

5 

103 

CCNVEC 


NEWSRC 

100019 

CORNCC 

SU 

5 

105 

GETCRD 


NEWSRC 

100019 

CORNCC 

SU 

5 

106 

CINV 


NEWSRC 

100019 

CORNCC 

SU 

5 

108 

RWCCF 


NEWSRC 

100019 

CORNCC 

SU 

5 

109 

PUTWFI 


NEWSRC 

100019 

CORNCC 

WR 

8 


YES 


NEWSRC 

100019 

CCNVEC 

SU 

1 


CCNVEC 


NEWSRC 

100020 

CCNVEC 

CO 

3 


MAFCOM 

M 

NEWSRC 

100020 

CCNVEC 

CO 

3 


AFIL 

C 

NEWSRC 

100020 

CCNVEC 

WR 

8 


YES 


NEWSRC 

100020 

GETCRD 

SU 

1 


GETCRD 


NEWSRC 

100021 

GETCRD 

CO 

3 


MAFCOM 

M 

NEWSRC 

100021 

GETCRD 

CO 

3 


AFIL 

C 

NEWSRC 

100021 

GETCRD 

WR 

8 


YES 


NEWSRC 

100021 

PUTWFI 

SU 

1 


PUTWFI 


NEWSRC 

100022 

PUTWFI 

CO 

3 


MAFCOM 

M 

NEWSRC 

100022 

PUTWFI 

CO 

3 


AFIL 

C 

NEWSRC 

100022 

RWCCF 

SU 

1 


RWCCF 


NEWSRC 

100023 

RWCCF 

CO 

3 


MAFCOM 

M 

NEWSRC 

100023 

RWCCF 

CO 

3 


AFIL 

C 

NEWSRC 

100023 

RWCCF 

RD 

7 


YES 


NEWSRC 

100023 

RWCCF 

WR 

8 


YES 


NEWSRC 

100023 

GETMPE 

SU 

1 


GETMPE 


NEWSRC 

100024 

GETMPE 

CO 

3 


MAFCOM 

M 

NEWSRC 

100024 

GETMPE 

CO 

3 


AFIL 

C 

NEWSRC 

100024 

GETMPE 

CO 

3 


CNTLFL 

C 

NEWSRC 

100024 

GETMPE 

CO 

3 


MAFCOM 

M 

NEWSRC 

100024 

GETMPE 

CO 

3 


SIZE 

C 

NEWSRC 

100024 

GETMPE 

SU 

5 

101 

HTFIO 


NEWSRC 

100024 

GETMPE 

SU 

5 

103 

MPERW 


NEWSRC 

100024 

MPERW 

SU 

1 


MPERW 


NEWSRC 

100025 

MPERW 

WR 

8 


YES 


NEWSRC 

100025 

SETFLP 

SU 

1 


SETFLP 


NEWSRC 

100026 

SETFLP 

CO 

3 


MAFCOM 

M 

NEWSRC 

100026 

SETFLP 

CO 

3 


AFIL 

C 

NEWSRC 

100026 

SETFLP 

SU 

5 

ioi 

TFLSET 


NEWSRC 

100026 

SETFLP 

SU 

5 

102 

HTFPAS 


NEWSRC 

100026 

RDTEMP 

SU 

1 


RDTEMP 


NEWSRC 

100027 

RDTEMP 

CO 

3 


MAFCOM 

M 

NEWSRC 

100027 

RDTEMP 

CO 

3 


AFIL 

C 

NEWSRC 

100027 

RDTEMP 

CO 

3 


MAFCOM 

M 

NEWSRC 

100027 

RDTEMP 

CO 

3 


SIZE 

C 

NEWSRC 

100027 

RDTEMP 

CO 

3 


MAFCOM 

M 

NEWSRC 

100027 

RDTEMP 

CO 

3 


CNTLFL 

C 

NEWSRC 

100027 

RDTEMP 

SU 

5 

101 

HTFIO 


NEWSRC 

100027 

RDTEMP 

WR 

8 


YES 


NEWSRC 

100027 

CWCRNT 

SU 

1 


CWCRNT 


NEWSRC 

100028 

CWCRNT 

CO 

3 


MAFCOM 

M 

NEWSRC 

100028 

CWCRNT 

CO 

3 


AFIL 

C 

NEWSRC 

100028 

CWCRNT 

SU 

5 

101 

FILREW 


NEWSRC 

100028 

CWCRNT 

SU 

5 

102 

RWCCF 


NEWSRC 

100028 

CWCRNT 

SU 

5 

103 

UNPAKF 


NEWSRC 

100028 

CWCRNT 

SU 

5 

104 

MPERW 


NEWSRC 

100028 

CWCRNT 

SU 

5 

105 

FCRNTP 


NEWSRC 

100028 

CWCRNT 

SU 

5 

106 

MAXTDF 


NEWSRC 

100028 

CWCRNT 

SU 

5 

107 

-WRCRNT 


NEWSRC 

100028 

UNPAKF 

SU 

1 


UNPAKF 


NEWSRC 

100029 

UNPAKF 

CO 

3 


MAFCOM 

M 

NEWSRC 

100029 

UNPAKF 

CO 

3 


AFIL 

C 

NEWSRC 

100029 

UNPAKF 

WR 

8 


YES 


NEWSRC 

100029 

FCRNTP 

SU 

1 


FCRNTP 


NEWSRC 

100030 

FCRNTP 

CO 

3 


MAFCOM 

M 

NEWSRC 

100030 

FCRNTP 

CO 

3 


AFIL 

C 

NEWSRC 

100030 

MAXTDF 

SU 

1 


MAXTDF 


NEWSRC 

100031 
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MAXTDF 

CO 

3 


MAFCOM 

M 

NEWSRC 

100031 

MAXTDF 

CO 

3 


AFIL 

C 

NEWSRC 

100031 

MAXTDF 

SU 

5 

101 

SORT 


NEWSRC 

100031 

MAXTDF 

WR 

8 


YES 


NEWSRC 

100031 

WRITMP 

SU 

1 


WRITMP 


NEWSRC 

100032 

WRITMP 

CO 

3 


MAFCOM 

M 

NEWSRC 

100032 

WRITMP 

CO 

3 


AFIL 

C 

NEWSRC 

100032 

WRITMP 

WR 

8 


YES 


NEWSRC 

100032 

CINV 

SU 

1 


CINV 


NEWSRC 

100033 

CINV 

CO 

3 


MAFCOM 

M 

NEWSRC 

100033 

CINV 

CO 

3 


AFIL 

C 

NEWSRC 

100033 

CINV 

WR 

8 


YES 


NEWSRC 

100033 

PRNTIM 

SU 

1 


PRNTIM 


NEWSRC 

100034 

PRNTIM 

CO 

3 


MAFCOM 

M 

NEWSRC 

100034 

PRNTIM 

CO 

3 


AFIL 

C 

NEWSRC 

100034 

PRNTIM 

WR 

8 


YES 


NEWSRC 

100034 

EMINMX 

SU 

1 


EMINMX 


NEWSRC 

100035 

EMINMX 

SU 

5 

101 

SORT 


NEWSRC 

100035 

FTIME 

SU 

1 


FTIME 


NEWSRC 

100036 

SSURCH 

SU 

1 


SSURCH 


NEWSRC 

100037 

SSURCH 

CO 

3 


MAFCOM 

M 

NEWSRC 

100037 

SSURCH 

CO 

3 


AFIL 

C 

NEWSRC 

100037 

SSURCH 

CO 

3 


MAFCOM 

M 

NEWSRC 

100037 

SSURCH 

CO 

3 


SIZE 

C 

NEWSRC 

100037 

SSURCH 

CO 

3 


MAFCOM 

M 

NEWSRC 

100037 

SSURCH 

CO 

3 


SPDAT 

C 

NEWSRC 

100037 

SSURCH 

SU 

5 

101 

KNTSPF 


NEWSRC 

100037 

SSURCH 

SU 

5 

102 

GNLSUR 


NEWSRC 

100037 

KNTSPF 

SU 

1 


KNTSPF 


NEWSRC 

100038 

KNTSPF 

CO 

3 


MAFCOM 

M 

NEWSRC 

100038 

KNTSPF 

CO 

3 


AFIL 

C 

NEWSRC 

100038 

KNTSPF 

SU 

5 

101 

FILREW 


NEWSRC 

100038 

KNTSPF 

SU 

5 

102 

RDSPF 


NEWSRC 

100038 

KNTSPF 

WR 

8 


YES 


NEWSRC 

100038 

RFILCR 

SU 

1 


RFILCR 


NEWSRC 

100039 

RFILCR 

CO 

3 


MAFCOM 

M 

NEWSRC 

100039 

RFILCR 

CO 

3 


ELMDAT 

C 

NEWSRC 

100039 

RFILCR 

SU 

5 

101 

RANSIZ 


NEWSRC 

100039 

RFILCR 

SU 

5 

102 

ZERORF 


NEWSRC 

100039 

ZERORF 

SU 

1 


ZERORF 


NEWSRC 

100040 

ZERORF 

WR 

8 


YES 


NEWSRC 

100040 

GNLSUR 

SU 

1 


GNLSUR 


NEWSRC 

100041 

GNLSUR 

CO 

3 


MAFCOM 

M 

NEWSRC 

100041 

GNLSUR 

CO 

3 


AFIL 

C 

NEWSRC 

100041 

GNLSUR 

CO 

3 


MAFCOM 

M 

NEWSRC 

100041 

GNLSUR 

CO 

3 


SIZE 

C 

NEWSRC 

100041 

GNLSUR 

CO 

3 


MAFCOM 

M 

NEWSRC 

100041 

GNLSUR 

CO 

3 


WORK 

C 

NEWSRC 

100041 

GNLSUR 

SU 

5 

101 

EINCOR 


NEWSRC 

100041 

GNLSUR 

SU 

5 

102 

RDSPF 


NEWSRC 

100041 

GNLSUR 

SU 

5 

103 

FNDELM 


NEWSRC 

100041 

GNLSUR 

WR 

8 


YES 


NEWSRC 

100041 

RDSPF 

SU 

1 


RDSPF 


NEWSRC 

100042 

RDSPF 

CO 

3 


MAFCOM 

M 

NEWSRC 

100042 

RDSPF 

CO 

3 


AFIL 

C 

NEWSRC 

100042 

RDSPF 

RD 

7 


YES 


NEWSRC 

100042 

FILREW 

SU 

1 


FILREW 


NEWSRC 

100043 

EINCOR 

SU 

1 


EINCOR 


NEWSRC 

100044 

EINCOR 

CO 

3 


MAFCOM 

M 

NEWSRC 

100044 

EINCOR 

CO 

3 


AFIL 

C 

NEWSRC 

100044 

EINCOR 

CO 

3 


MAFCOM 

M 

NEWSRC 

100044 

EINCOR 

CO 

3 


ELMDAT 

C 

NEWSRC 

100044 

EINCOR 

SU 

5 

101 

RWELMS 


NEWSRC 

100044 

FNDELM 

SU 

1 


FNDELM 


NEWSRC 

100045 

FNDELM 

CO 

3 


MAFCOM 

M 

NEWSRC 

100045 

FNDELM 

CO 

3 


AFIL 

C 

NEWSRC 

100045 

FNDELM 

CO 

3 


MAFCOM 

M 

NEWSRC 

100045 
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FNDELM 

CO 

3 


ELMDAT 

C 

NEWSRC 

100045 

FNDELM 

CO 

3 


MAFCOM 

M 

NEWSRC 

100045 

FNDELM 

CO 

3 


SIZE 

C 

NEWSRC 

100045 

FNDELM 

CO 

3 


MAFCOM 

M 

NEWSRC 

100045 

FNDELM 

CO 

3 


SPDAT 

C 

NEWSRC 

100045 

FNDELM 

su 

5 

101 

GETWND 


NEWSRC 

100045 

FNDELM 

su 

5 

102 

CKWIND 


NEWSRC 

100045 

FNDELM 

su 

5 

104 

I3DSF 


NEWSRC 

100045 

FNDELM 

su 

5 

105 

RWELMS 


NEWSRC 

100045 

FNDELM 

su 

5 

109 

RWSPCF 


NEWSRC 

100045 

FNDELM 

WR 

8 


YES 


NEWSRC 

100045 

GETWND 

su 

1 


GETWND 


NEWSRC 

100046 

GETWND 

CO 

3 


MAFCOM 

M 

NEWSRC 

100046 

GETWND 

CO 

3 


AFIL 

C 

NEWSRC 

100046 

GETWND 

CO 

3 


MAFCOM 

M 

NEWSRC 

100046 

GETWND 

CO 

3 


ELMDAT 

C 

NEWSRC 

100046 

GETWND 

CO 

3 


MAFCOM 

M 

NEWSRC 

100046 

GETWND 

CO 

3 


SIZE 

C 

NEWSRC 

100046 

GETWND 

CO 

3 


MAFCOM 

M 

NEWSRC 

100046 

GETWND 

CO 

3 


WORK 

C 

NEWSRC 

100046 

GETWND 

su 

5 

lOI 

RWELMS 


NEWSRC 

100046 

CKWIND 

su 

1 


CKWIND 


NEWSRC 

100047 

MKRGF 

su 

1 


MKRGF 


NEWSRC 

100048 

MKRGF 

CO 

3 


MAFCOM 

M 

NEWSRC 

100048 

MKRGF 

CO 

3 


AFIL 

C 

NEWSRC 

100048 

MKRGF 

CO 

3 


MAFCOM 

M 

NEWSRC 

10Q048 

MKRGF 

CO 

3 


ELMDAT 

C 

NEWSRC 

100048 

MKRGF 

CO 

3 


MAFCOM 

M 

NEWSRC 

100048 

MKRGF 

CO 

3 


SIZE 

C 

NEWSRC 

100048 

MKRGF 

CO 

3 


MAFCOM 

M 

NEWSRC 

100048 

MKRGF 

CO 

3 


CNTLFL 

C 

NEWSRC 

100048 

MKRGF 

su 

5 

101 

RWNODS 


NEWSRC 

100048 

MKRGF 

su 

5 

103 

EMINMX 


NEWSRC 

100048 

MKRGF 

su 

5 

104 

RWELMS 


NEWSRC 

100048 

MKRGF 

su 

5 

105 

PROFAC 


NEWSRC 

100048 

MKRGF 

RD 

7 


YES 


NEWSRC 

100048 

MKRGF 

WR 

8 


YES 


NEWSRC 

100048 

PROFAC 

su 

1 


PROFAC 


NEWSRC 

100049 

PROFAC 

CO 

3 


MAFCOM 

M 

NEWSRC 

100049 

PROFAC 

CO 

3 


AFIL 

C 

NEWSRC 

100049 

PROFAC 

CO 

3 


MAFCOM 

M 

NEWSRC 

100049 

PROFAC 

CO 

3 


ELMDAT 

C 

NEWSRC 

100049 

PROFAC 

CO 

3 


MAFCOM 

M 

NEWSRC 

100049 

PROFAC 

CO 

3 


SIZE 

C 

NEWSRC 

100049 

PROFAC 

su 

5 

101 

RWELMS 


NEWSRC 

100049 

PROFAC 

RD 

7 


YES 


NEWSRC 

100049 

PROFAC 

WR 

8 


YES 


NEWSRC 

100049 

SKIP 

su 

1 


SKIP 


NEWSRC 

100050 

SKIP 

CO 

3 


MAFCOM 

M 

NEWSRC 

100050 

SKIP 

CO 

3 


SIZE 

C 

NEWSRC 

100050 

SKIP 

RD 

7 


YES 


NEWSRC 

100050 

HTFIO 

su 

1 


HTFIO 


NEWSRC 

100051 

HTFIO 

RD 

7 


YES 


NEWSRC 

100051 

RDCNTL 

SU 

1 


RDCNTL 


NEWSRC 

100052 

RDCNTL 

CO 

3 


MAFCOM 

M 

NEWSRC 

100052 

RDCNTL 

CO 

3 


AFIL 

C 

NEWSRC 

100052 

RDCNTL 

CO 

3 


MAFCOM 

M 

NEWSRC 

100052 

RDCNTL 

CO 

3 


CNTLFL 

C 

NEWSRC 

100052 

RDCNTL 

CO 

3 


MAFCOM 

M 

NEWSRC 

100052 

RDCNTL 

CO 

3 


CNTLTM 

C 

NEWSRC 

100052 

RDCNTL 

CO 

3 


MAFCOM 

M 

NEWSRC 

100052 

RDCNTL 

CO 

3 


SIZE 

C 

NEWSRC 

100052 

RDCNTL 

su 

5 

101 

HTFIO 


NEWSRC 

100052 

RDCNTL 

su 

5 

102 

ELTRAN 


NEWSRC 

100052 

RDCNTL 

RD 

7 


YES 


NEWSRC 

100052 

RDCNTL 

WR 

8 


YES 


NEWSRC 

100052 

TFLSET 

SU 

1 


TFLSET 


NEWSRC 

100053 
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TFLSET 

CO 

3 


MAFCOM 

M 

NEWSRC 

100053 

TFLSET 

CO 

3 


AFIL 

C 

NEWSRC 

100053 

TFLSET 

CO 

3 


MAFCOM 

M 

NEWSRC 

100053 

TFLSET 

CO 

3 


CNTLTM 

C 

NEWSRC 

100053 

TFLSET 

SU 

5 

101 

SKIP 


NEWSRC 

100053 

TFLSET 

RD 

7 


YES 


NEWSRC 

100053 

GENTMP 

SU 

1 


GENTMP 


NEWSRC 

100054 

GENTMP 

CO 

3 


MAFCOM 

M 

NEWSRC 

100054 

GENTMP 

CO 

3 


MAFCOM 

M 

NEWSRC 

100054 

GENTMP 

CO 

3 


AFIL 

C 

NEWSRC 

100054 

GENTMP . 

CO 

3 


MAFCOM 

M 

NEWSRC 

100054 

GENTMP 

CO 

3 


SIZE 

C 

NEWSRC 

100054 

GENTMP 

CO 

3 


MAFCOM 

M 

NEWSRC 

100054 

GENTMP 

CO 

3 


WORK 

C 

NEWSRC 

100054 

GENTMP 

CO 

3 


MAFCOM 

M 

NEWSRC 

100054 

GENTMP 

CO 

3 


CNTLTM 

C 

NEWSRC 

100054 

GENTMP 

CO 

3 


CNTLFL 

C 

NEWSRC 

100054 

GENTMP 

SU 

5 

101 

TFLSET 


NEWSRC 

100054 

GENTMP 

SU 

5 

102 

GETMPE 


NEWSRC 

100054 

GENTMP 

SU 

5 

103 

SETFLP 


NEWSRC 

100054 

GENTMP 

SU 

5 

104 

RDTEMP 


NEWSRC 

100054 

GENTMP 

SU 

5 

105 

CWCRNT 


NEWSRC 

100054 

GENTMP 

SU 

5 

106 

ETCORE 


NEWSRC 

100054 

GENTMP 

SU 

5 

107 

NTCORE 


NEWSRC 

100054 

GENTMP 

SU 

5 

108 

STSTMP 


NEWSRC 

100054 

STSTMP 

SU 

1 


STSTMP 


NEWSRC 

100055 

STSTMP 

CO 

3 


MAFCOM 

M 

NEWSRC 

100055 

STSTMP 

CO 

3 


MAFCOM 

M 

NEWSRC 

100055 

STSTMP 

CO 

3 


AFIL 

C 

NEWSRC 

100055 

STSTMP 

CO 

3 


MAFCOM 

M 

NEWSRC 

100055 

STSTMP 

CO 

3 


SIZE 

C 

NEWSRC 

100055 

STSTMP 

CO 

3 


SPDAT 

C 

NEWSRC 

100055 

STSTMP 

SU 

5 

101 

RWSPCF 


NEWSRC 

100055 

STSTMP 

SU 

5 

102 

WRSTMP 


NEWSRC 

100055 

WRSTMP 

SU 

1 


WRSTMP 


NEWSRC 

100056 

WRSTMP 

CO 

3 


MAFCOM 

M 

NEWSRC 

100056 

WRSTMP 

CO 

3 


MAFCOM 

M 

NEWSRC 

100056 

WRSTMP 

CO 

3 


AFIL 

C 

NEWSRC 

100056 

WRSTMP 

WR 

8 


YES 


NEWSRC 

100056 

ETCORE 

SU 

1 


ETCORE 


NEWSRC 

100057 

ETCORE 

CO 

3 


MAFCOM 

M 

NEWSRC 

100057 

ETCORE 

CO 

3 


MAFCOM 

M 

NEWSRC 

100057 

ETCORE 

CO 

3 


AFIL 

C 

NEWSRC 

100057 

ETCORE 

CO 

3 


MAFCOM 

M 

NEWSRC 

100057 

ETCORE 

CO 

3 


SIZE 

C 

NEWSRC 

100057 

ETCORE 

SU 

5 

101 

FILREW 


NEWSRC 

100057 

ETCORE 

RD 

7 


YES 


NEWSRC 

100057 

NTCORE 

SU 

1 


NTCORE 


NEWSRC 

100058 

NTCORE 

CO 

3 


MAFCOM 

M 

NEWSRC 

100058 

NTCORE 

CO 

3 


MAFCOM 

M 

NEWSRC 

100058 

NTCORE 

CO 

3 


AFIL 

C 

NEWSRC 

100058 

NTCORE 

CO 

3 


MAFCOM 

M 

NEWSRC 

100058 

NTCORE 

CO 

3 


SIZE 

C 

NEWSRC 

100058 

NTCORE 

SU 

5 

101 

HTFIO 


NEWSRC 

100058 

WRCRNT 

SU 

1 


WRCRNT 


NEWSRC 

100059 

WRCRNT 

CO 

3 


MAFCOM 

M 

NEWSRC 

100059 

WRCRNT 

CO 

3 


MAFCOM 

M 

NEWSRC 

100059 

WRCRNT 

CO 

3 


AFIL 

C 

NEWSRC 

100059 

WRCRNT 

WR 

8 


YES 


NEWSRC 

100059 

HTFPAS 

SU 

1 


HTFPAS 


NEWSRC 

100060 

HTFPAS 

CO 

3 


MAFCOM 

M 

NEWSRC 

100060 
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